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While I feel that the general nature of the picture shown on Figure 
1 and Figure 2 is correct, some of the details such as the wider mushroom 
fallout area at Seattle than Denver, for example, are open to suspicion 
since the data used to prepare this chart are not homogeneous. That is, 
at each elevation, we used the best possible average wind, but at, say, 
10,000 feet, we used a mean chart based on winds collected from about 10 
years of Observations,while at 80,000 feet average winds based on records 
for only 4 years were available. Further, we are always plagued with 
the fact that high level winds are obtained mainly during periods of fair 
weather and low wind speeds, while low-level winds are collected more 
frequently. 

Figure 2, the summer patterns, also shows predominant fallout to 
the east of the burst. The wind epeeds become much lighter and the west 
coast stations actually have a mean northward airflow. The southern tier 
of States is located in an area of very light winds. Thus, at New Orleans, 
for example, the mean fallout picture shows no preferential direction. 
Western Europe would be more like eastern U. S. rather than the U. S. 

West Coast. 

Having viewed the mean picture, we now pass on to the variability 
of individual cases about this mean. In the troposphere, or up to 30 - ' 
45,000 feet, where one is more familiar with wind patterns, it is well 
recognized that above the lowest levels (say above 10,000 feet), the 
prevailing westerly flow is the result of waves or undulations super¬ 
imposed on a mean zonal flow. At times, the daily picture in the mid¬ 
troposphere looks much like the mean with almost straight west-to-east 
flow, while at other times, we find large amplitude waves or closed 
centers. 


These waves or closed centers most frequently move across the 
country from west to east, but occasionally, at certain preferred places 
such as the southwest U. S., remain stationary for many days. 

In the winter, in the stratosphere (Figure 3), the winds normally- 
also blow west-to-east ov^r the temperate latitudes. Although this map 
shows the flow only at 53*000 feet, direction of the airflow is much the 
same up to 100,000 feet. However, on anomolous occasions, as in Figure 
4, which shows the airflow at 67,000 feet in February 1952, large areas 
of the temperate latitudes possess easterly winds. In this map, northern 
U. S. and southern Canada is such an area in which the normal wind direc¬ 
tion is reversed. Fallout in such a case would probably still be to the 
east of ground zero in Northern U. S., but not nearly as far away from 
ground zero as in the average case. . _ 

DOE ARCHIVES 

However, the typical and almost invariant summer picture above about 
60,000 feet (Figure 5 shows east-to-west flow). This will result in many 
cases of fallout from the stratosphere occurring to the west of ground zero 
because of light westerly tropospheric winds. 
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Figure U 

Anomolous Stratospheric Flow Pattern in Winter, with 
Easterly Winds over much of the United States and Canada 










Figure 6 shows the distribution of points at which fallout will 
occur for hypothetical bombs detonated at Washington, D. G. (where the 
lines intersect) for January and July and coming from two elevations; 
one in the middle or upper troposphere ( 30,000 feet), and one in the 
lower stratosphere (60,000 feet). The points actually indicate where 
a 100 micron particle of 2.5 density would land if started at the indi¬ 
cated height over ground zero. The arcs are at 200-mile intervals. In 
the winter, while there ere pretty wide departures from direct west-to- 
east fallout, there are no cases of deposition to the west of ground 
zero during the two Januarys used in this study from either level. During 
the summer, not only does the 100 micron particle fall much closer to 
ground zero, but there actually are cases of fallout to the left or the 
west of ground zero. The filled-in circles show the cases in which it 
was raining or snowing at the time of the wind observation. In the 
winter, it is evident that precipitation is fair more likely with fallout 
occurring northeast of the burst. The distribution of the points provides 
a measure of the scatter of directions of, say, the effective mean winds 
for the stem in the left side, and lower mushroom on the right side of 
the figure. 

The last bit of climatological statistics appears in Figure 7* Since 
the .distant fallout is felt to come from the mushroom, the size of the 
area in which potential danger exists depends, in part, on the amount of 
the directional shear of the winds in the mushroom. Although the clima¬ 
tological data for this figure was taken from Washington, D. C., the 
results are drawn for Pittsburgh, Pa. We have on this figure the cumula¬ 
tive frequencies of directional shears in the mushroom. Thus, in January, 
there is almost a 60 % likelihood that the directional shear will be no 
greater than the smallest angle of 5° and almost likelihood ttiat it 
wil-l he no greater than 10°, etc. The shear is seen to be much greater 
in summer. For each of the angles, one may get an idea of the fraction 
of the northeastern U. S. on which fallout would be observed. 

The conclusion to be drawn from this figure is the fact that there 
is only a small likelihood of large directional shear} especially in 
winter, a saving factor for some Civil Defense considerations. I have 
no statistics to confirm this state of affairs everywhere in the temperate 
latitudes, however. 

Figures 8-13 illustrate actual cases of fallout with typical 
weather situations. For example, in Figure 8, on the left is the weather 
situation near Washington, D. C., at 0630 (Greenwich Time) on 10 January 
1953 showing the isobars and fronts in the usual fashion and the precipit¬ 
ation area in shading. A storm is approaching Washington from the south 
with a wide precipitation area to the north and west. To the right of 
the weather picture, we find the fallout plot in the event- of a hypothet¬ 
ical burst at, say, the Pentagon. The radial lines are the lines along 






Figure 6 

f Fall-out free 30*000 and 60,000 feet in January and July, Washington, D 









which fallout from the Indicated height will he deposited. The heavy line 
shows the places on the ground at which the 100 micron particle of a 
specific gravity 2.5 will be deposited, while the dotted line shows where 
fallout will reach at 8 hours after burst. The numbers on the fallout 
diagram are the heights in tens of thousands of feet: 10,000 , 20,000, 

30,000 feet and so forth. The direction of motion of the stem below 
10,000 feet can be deduced from the surface isobars. No attempt is made 
to extrapolate beyond the height of observed winds. 

It is evident that fallout in Virginia, Maryland, and southern 
Pennsylvania will be in the rain. No further interpretation of the 
significance of this fact will be attempted at this time. This is an 
example of a cyclonic disturbance approaching the target from the south 
or southwest and shows the fallout to occur typically to the northwest 
from the lowest levels and the north, or northeast from the middle and 
upper levels. 

Other typical winter weather situations are shown in Figures 9 ** 11« 
Figure 9 shows the fallout pattern following the passage of a winter cold 
front, with fallout occurring to the southeast of Washington. Figure 10 
illustrates the conditions well within the cold air, with northwesterly 
winds in the lowest levels, shifting to southwesterly at the higher levels. 
Figure 11 typifies conditions after the center of the high has passed 
Washington and a new low is approaching from the west. Figures 12 and 13 
show two examples of summer fallout patterns. Figure 12 illustrates 
conditions following the passage of a weak summer cold front with north¬ 
westerly winds prevailing in the lower levels and very weak easterlies 
in the stratosphere, while Figure 13 shows a well developed Bermuda high 
with southerly winds near the surface becoming easterly inlfae stratosphere 
resulting in fallout to the northeast of Washington. 

1 should like to complete the presentation by returning to the role 
which meteorology plays in fallout and take a look into the future. The 
two contributions which the meteorologist can make, as mentioned before, 
are: Collection of climatological fallout data for planning, and opera¬ 
tional fallout forecasts for day-to-day use. I feel that the presentation 
I have made this morning answers only a very limited number of planning 
questions. The kind of problem it does not answer is: If the Pentagon 
■were the target for an x megaton bomb detonated at the ground, what is 
the probability that New York or Philadelphia or Boston might be subjected 
to a dangerous fallout? If this ^rmposium provides a method by which one 
may either delineate dosages from the wind data or even delineate, say, 
two kinds of areas showing different dangers, we feel that we can answer 
the kind of question which Civil Defense or military commanders are prone 
to ask. The simple plan to be described represents ideas contributed by 
many groups besides our own and machine methods may be required in its 
execution. 



Figure 8 

Surface Weather Map and Associated Fall-out Pattern, January 10, 1953 






Figure 9 

Surface Weather Map and Associated Fall-out Pattern, January 12, 1953 






Figure 10 

Weather Map and Associated Fall-out Pattern, January 12, 1954 






Pattern , January 14 > 1954 







4 ,. Figure 12 

Surface Weather Map and Associated Fall-out Pattern, July 2, 1954- 





Surface Weather Map and Associates Fall-out Pattern, July 14, 1954 







We would first divide the area of interest, for example, the United 
States, into regions with what we consider to be homogeneous areas of 
fallout climatology, an example of which is shown in Figure l4. A single 
upper air wind reporting station would be selected in each region and 
considered to be representative of the entire region. As we see in Figure 
15, fallout diagrams with the areas of, say, dangerous and very dangerous 
fallout computed from each weather station would then be prepared each day 
for a number of months of past weather records. For example, the heavy 
shading might illustrate the very dangerous areas and be surrounded by 
dangeroils areas. On each day, we overlay the same grid of points and then 
singly count the frequency of occurrence of the point in each area. These 
nay then be tabulated for each season and analyzed to provide probability 
distributions shown on the right. When the probability distribution for 
occurrences of very dangerous fallout is superimposed on a potential 
target, one can readily read off the percentage probability at, say, a 
nearby city and obtain the likelihood that it would lie within the very 
dangerous fallout region for the given target. 

This, I believe, would satisfy very many planning requirements. A 
different story must be told for the day-to-day operational capabilities 
of the meteorologist. 

First of all, even if we use persistency as a forecasting tool - 
that is, take wind observations, say, every 6 hours and use the last 
observation to transport the fallout debris, I am not even sure our 
present observational network is adequate. There are two main defects: 
first, the average height of the wind runs is only 55>000 feet and even 
at Silver Hill here in the Washington area where ideal operating condi¬ 
tions exist, the average height of the wind runs is only 75>OOQ feet. 
Frequently, we do get runs in excess of 100,000 feet, but these are 
us ual ly in warm climates and with wind speeds which are not excessive. 

My guess is that, in the winter, with current equipment, the average 
height of the upper wind observations over northern U. S. is well below 
55,000 feet. Second, if the balloon is too near the horizon, it either 
cannot be tracked at all or the wind data becomes uncertain. Thus, if 
the speeds are high, making low elevation angles, our wind data becomes 
much poorer. In addition, the present method of analyzing and reporting 
upper winds is not particularly suited for fallout computations and cer¬ 
tain unnecessary errors are introduced. In Figure 1 6 , we see part of an 
actual plot of wind speed with height. What is normally used in preparing 
fallout plots are the winds at 5>000-foot levels rather than vertically 
integrated winds. The latter would be more desirable from a fallout 
viewpoint, since winds at all levels transport the particles. Although 
chosen at random, you can see in Figure l6 that the speed at each 5,000- 
ft. level is lighter, by chance, than the mean through the layer in which 
the wind is centered so that the predicted fallout would fall short of 
where it would actually be. 
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Figure 14 

Examples of Areas with Similar Fall-cut Climatology 









It is highly unlikely, if a bomb is actually detonated or we know 
one will be detonated, that a persistency forecast will be acceptable. 
Weathermen will be asked to make best possible forecasts where so much 
is at stake. Unfortunately, we have practically no forecasting experi¬ 
ence at levels above about 30,000 feet, and, in fact, do not routinely 
even draw weather maps above 55,000 feet. Our experience, in fact, is 
so l imi ted to the lover atmosphere that, at the Nevada Test Site, the 
absolute .wind forecast error doubles from the 20,000 layer to the 40,000- 
foot layer. Further, time and space variability of the wind can be ade¬ 
quately included only by the meteorologist. Concerning the need to 
incorporate time and space changes in the wind, I have an interesting 
case showing a comparison of using point winds with time and space vari¬ 
able winds, in Figure 17• The straight dashed lines are the lines along 
which fallout occurs from a given height if the test site wind at burst 
time alone is used at Nevada while the curved solid lines incorporate 
time and space variability. The numbers are the observed fallout along 
the roads. It is clear that the time and space variability is superior, 
although it may not always provide as good a forecast as shown here. The 
need for time and space variability is much greater in the case of mega¬ 
ton weapons than in Nevada test operations since fallout is important 
for longer times and to greater distances. 

Finally, the meteorologist must be able to provide the Civil Defense 
official or military commander with a measure of the uncertainty in the 
forecast. We can all imagine cases in which, if the uncertainty in fall¬ 
out forecast is too great, it may do more harm than good. 

In conclusion, the meteorologist even today can provide certain 
valuable generalities about the probable areas of radioactive fallout, 
but much remains yet to be done to answer questions about dosages or 
dangerous areas and about the details of the fallout probabilities. 
Finally, the weatherman needsmuch more development research and training 
before he can really assist in operational fallout forecasting. 
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Figure 17 

Effect, of Time and Space Variability of the Wind on the Location of Fall-out 



PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS 
FUNCTIONS OF YIELD AND METEOROLOGICAL CONDITIONS, 
AIR WEATHER SERVICE METHOD 

Lt. Col. Clifford A. Spohn, USAF 
Air Weather Service 


My talk today is, in a sense, sailing under false colors. The 
problem under consideration during this portion of the symposium deals 
with the determination of ground levels of contamination resulting from 
fallout, and I am listed in the program as presenting the Air Weather 
Service method. Actually, the method of fallout prediction used by Air 
Weather Service does not provide an answer giving ground contamination 
levels. The technique currently used is concerned only with the meteor¬ 
ological portion of the problem, that is, specifying where and when fall¬ 
out will occur. 

The technique itself is a very simple one, and is based on the 
following assumptions: 

That the point winds near the detonation point (observed or 
forecast) are representative of the time-space wind field through which 
fallout takes place. 

That winds at successive 10,000-foot layers in the atmosphere 
are representative of the wind flow through their respective layers. 

That the various contaminated particles will have constant fall 
rates, although the fall rates are unspecified, and need not be related 
to particle size. 

That fallout takes place upon a level surface, i.e. terrain 
Variations are ignored. 

Given these assumptions, the fallout plot is easily p£2£® 
wind sounding, observed or forecast, for the place in question, is used 
to construct composite wind vectors. The 5000' wind vector (representing 
flow in the layer from sea level to 10,000') is laid off from the ground 
zero position. From the end of this vector, the 15,000' wind vector (for 
the layer 10,000’ - 20,000') is next laid off. The 25,000' wind vector 
is laid off from the end of the 15 , 000 * wind vector, and the process 
continued to as high a level as desired, or to the end of the wind data. 
Radial lines are now srawn from ground zero to the ends of the various 
wind vectors and labeled appropriately. Thus, the radial line to the 
end of the 15 , 000 ' wind vector is labeled 20,000'; to the end of the 
25,000' wind vector ,30,000', etc. These radial lines are the resultant 


wind vectors from the specified level to the ground, and indicate the 
line away from ground zero along which fallout from the given level will 
take place. The lengths of these resultant vectors is measured and 
divided by their labeled heights expressed in tens of thousands of feet 
(i.e., the length of the resultant vector labeled 20,000' will be divided 
by "2", the vector divided by "3", etc.). The numbers obtained 

represent the speed at which fallout will progress outward from ground 
zero along the various radial lines. 

The finished plot is prepared by laying off the resultant vectors, 
and marking off along them successive hourly positions of fallout. By 
connecting all of the points corresponding to a given time after zero, 
the line along which fallout will occur at that time is indicated. 

The completed plot thus shows three items. First, the general area 
in which fallout will take place; second, the lines along which fallout 
from a given level will occur; and, third, the lines along which fallout 
will occur at particular times. 

The validity of the underlying assumptions must certainly be 
considered. The latter three assumptions, that of a constant fall rate, 
fallout on a level surface, and using a wind as representative of a 
10,000' layer, will not, in general, introduce much error. The first 
assumption, however, that a spot wind sounding will be representative 
of the space-time wind field throughout which fallout takes place, is 
less dependable. On a comparatively small test sample from the Nevada 
Proving Ground, using wind soundings taken two hours prior to detona¬ 
tion, and for fallout computed out to not more than nine hours after 
detonation, 75$ of the resultant vectors deviated by less than 15° from 
a post-analyzed position. If the fallout computations are extended 
considerably beyond nine hours, or if old wind data (i.e., 24-hour 
forecasts or observations 6-12 hours old) are used, this assumption 
becomes increasingly less reliable, and may lead to quite serious errors 
in the calculated fallout plot. For this reason, current Air Weather 
Service practice is to use, if at all possible, observed wind data less 
than six hours old, and to compute fallout out to 6 - 9 hours from 
time zero. ^ 
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Note: Examples of problem solutions by the Air Weather Service method 

are given in the paper, "Comparison of Problem Solutions by Various 
Methods," by Cdr. R. W. Paine. 
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PREDICTION OF DOSE-FATE AND DOSAGE CONTOURS AS 
. FUNCTIONS OF YIELD AND METEOROLOGICAL CONDITIONS, 
NAVY WEATHER SERVICE METHOD 

Dr. F. W. van Straten 
Office of the Chief of Naval Operation* 


The Navy method for predicting the location of potential radio¬ 
active areas after a surface burst of an atomic weapon over land or 
shallow water is designed to permit routine computations in advance 
of an actual explosion. It is expected that transparent overlays 
showing potential radioactive areas will be prepared, when necessary, 
after each wind sounding. These overlays will then be available for 
superimposticn on a map should an atomic explosion occur during the 
valid time of the wind sounding. Although it is realized that this 
procedure gives only approximate results, it is felt that immediate 
operational decisions can be made on the basis of these approximations. 
Time and data permitting, more refined computations can be made 
subsequent to the event. 

To permit this computation in advance of an explosion - when the 
height of the explosion is unknown, the size and nature of the bomb 
is not known and the terrain characteristics are undetermined - a 
number of assumptions is required. The limitations resulting from 
these assumptions are acceptable in view of the availability of a 
rough diagram from which immediate tactical decisions can be made. 


It is assumed that the explosion penetrates the tropopause. This 
assumption is made so that the diagram will show the maximum area 
which my be affected. If it is determined, subsequently, that the 


radioactive cloud did not, in fact, penetrate the tropopause, the 


diagram can later be reduced. 
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Lack of knowledge concerning the nature of the bomb and terrain 
necessitates assumptions about the distribution of particles in the 
radioactive cloud and their rate of radioactive decay. It is assumed 
that any horizontal cross section through the stem is identical with 
any other horizontal cross section. In other words, it is assumed 
that all particles sizes are distributed uniformly throughout the 
cloud. Further, radioactive decay is postulated to occur at a rate 
proportional to the square of the diameter of the particle and in¬ 
versely proportional to the square of the time from the explosion. 


Although it is recognized that a complete spectrum of particle 
sizes is probably present, only a size range between 7f> and lf>0 
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microns is considered. It is assumed that particles larger than 150 
microns will fall out so close to ground aero as to produce no 
additional hazard. The smaller size particles, less than 75 microns, 
are of concern, principally, because of the great distances they can 
travel in the time it takes them to fall from high altitudes. For 
particles smaller than 75 microns, this tine is so great that consider¬ 
able modification in the wind field is possible, and a computation 
based on a current wind sounding may lead to completely erroneous 
results. Particles of 75-micron size are assumed to fall at a rate 
of 5,000 feet per hour, while those of 150-micron size fall at a rate 
of 5,000 feet in 3/10 of an hour. 

For the original construction of the potential radioactive area, 
it is assumed that validity is limited to 6 hours in time, and 250 
miles in radius. 

It should be noted that provision is made for subsequent extension 
in time and area. Further, while the diagram plotted in advance is 
designed to show fall-out on the surface, the method also lends it¬ 
self to subsequent computations to show radioactive areas at any 
altitude and time - air radox plots - following the explosion* 


In brief, the method used by the Navy consists of the following 
steps t 


Determination of the mean wind vectors between specified 

levels• 

From the mean wind vector and the rate of fall, the position 
of fall-out and the time taken to fall are computed for each size 
particle. 

The positions marked along the mean wind vectors permit 
definition of the area of fall-out. 

The time marks can be used to draw lines showing the 
approximate time at which radioactivity reaches the surface within 
the affected area. 

The degree of radioactive contamination in any given unit 
area is considered to bs proportional to the number of mean vectors 
traversing that area. From the vector concentration and the time 
lines, it is possible to compare relative radioactivity from unit area 
to unit area, A series of lines is drawn on the diagram connecting 


area of equal radioactivity. 
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In practice a conventional winds aloft plotting board, a trans¬ 
parent plastic sheet and a radex rule are used. Figure 1 is a dia¬ 
grammatic sketch of the radex rule. 
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RADEX RULE 
ML— 502 /UM 


Figure 1 



plot of surface and 5OO0‘ 
VECTOR ON PLOTTING BOARD 
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Figure 2 
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Across the base of the radex rule is a scale used for plotting 
-wind speeds* When a given mean wind vector is placed under the radex 
rule in such a way that it is parallel to the horizontal lines on the 
rule and extends from the vertical side to the dark line hypotenuse, 
the slot in the radex rule marks the position of fall-out of the 150 
micron particles and the outer hypotenuse - edge - marks the position 
of fall-out of the 75 micron particles. The line extending fVom 
zero on the base scale to the hypotenuse, together with the tick marks 
on either side of this line, is used to expand a plotted diagram by 10° 
to take care of dispersion and minor errors in wind information. 

The following specific steps are followed to produce the Navy 
diagram: 

On the plotting board, add the 5,000 foot wind vector to l/2 
the surface wind vector. On the end of the 5,000 foot vector add the 
10,000 foot vector. Continue this process with each 5,000 foot wind 
vector to the top of the sounding* The radex rule is used for marking 
the wind speeds. See figure 2* 

Place a transparent, plastic sheet over the plotting board. 
Mark the center of the plotting board and then mark the 100° position 
of the plotting board as north. Transcribe the mid point of each wind 
vector on the plastic sheet. Thus the mid point of the 5,000 foot 
vector is transcribed and labeled 5. The mid point of the 10,000 foot 
vector is labeled 10, etc. 

Connect the center indicated on the plastic sheet with the 
points marked 5, 10, 15, etc., extending the line well beyond the 
marked points. 

Place the radex rule on each vector in turn so that the center 
indicated on the plastic sheet is on the vertical edge of the rule, 
and the numbered point is on the dark line hypotenuse with the vector 
parallel to the horizontal lines. Mark through the slot on to the 
vector and mark the edge of the radex rule on the vector. See figure 

3. 

For the vector labeled 5 the slot marking is labeled 3 indi¬ 
cating 3/10 of an hour. The outer mark is labeled 10 indicating one 
hour. On the vector labeled 10 these corresponding points are 6 and 
20. The 15 vector shows 9 and 30, etc. See figures 4 and 5. 

Using the radex rule, expand the diagram by adding 10° on 
either side of the outer vectors. See figure 6. 

Estimate the position of "10 w - the one hour position - on 
each of the vectors, and draw a line through these points. Repeat 
for 2 hours, 3 hours, etc. 

Enclose an area on the overlay limited by the 10° expansion 
lines, the vector ends, the 6-hour line arid the 250 mile radius arc, 
whichever is smallest (the scale of the radex rule is such that 20 M 
is equal to 250 miles), 
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Radex Rule in position on Vector 


Figure 3 
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TABLE I 










Place a rectangular grid over the diagram and count up the 
nuster of vectors traversing each grid area. The grid is drawn in 
squares 8/10 of an inch frcm the side. On the grid draw a line 
connecting the 3> 6* 9, etc., points. Those portions of the vectors 
which fall closer to the center than this line represent large 
particle fall-out. The remainder of the vector represents small 
particle fall-out. Using Table I, the number of vectors and the tine 
lines calculate the relative radioactivity in each area. Connect areas 
of equal indicated activity. See figures 7 and 8, 

The diagram is now completed since it shows the area of probable 
fall-out, the time of fall-out and the relative activity. 

The steps for preparing an air radex plot will not be outlined 
in detail. Suffice it to say that appropriate mean wind vectors ar® 
selected. For example, if an air radex plot at 5,000 feet and H + 1 
hours is required, small particle fall-out will originate from the 
10,000 foot level and large particle fall-out from the 20,000 foot 
level. Thus, only vectors between 5,000 and 20,000 feet are used. 

To expand the diagram timewise, an appropriate air radex plot is 
first constructed and then a new wind field used to determine how 
particles fall from the selected level to the surface. 

One final note is required, it is assumed that the intensity of 
radiation is proportional to the ratio of the square of the diameter 
of the particle to the square of time. This ratio is given as 
relative activity. To convert relative activity to actual activity it 
is necessary to have a proportionality constant. Presumably a single 
measurement anywhere within the area can provide this constant and 
will permit conversion of all relative values to actual values. 
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Figure 





PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS 
FUNCTIONS CF TIELD AND METEOROLOGICAL CONDITIONS, 

ARMY SIGNAL CORPS METHOD 

' Donald M. Swingle 

Signal. Corps Engineering Laboratories 

The Army Problem 

During the latter half of 19$h the Signal Corps received an urgent 
requirement to provide fall-out computations for a study. The two 
Signal Corps organizations interested in meteorology, i.e«, the Signal 
Corps Engineering Laboratories and the Array Electronic Proving Ground, 
united to perform the required job. SC EL was to determine or develop 
a suitable fall-out computation method and AELCTPG was to apply the 
method to the computation of the patterns required. 

The requirements of the Army problem called for the quantitative 
computation of radioactive fall-out material and dosage under a number 
of different meteorological conditions . A survey of existing computa¬ 
tion methods revealed none that would obtain the required results. It 
was, therefore, decided to devise a suitable method subject to the 
following requirements: 

The method must be physically reasonable, starting with a 
specified initial state of the cloud, accounting for modifying influ¬ 
ences occurring during fall-out, and yielding quantitative plots of the 
amount and times of arrival of fall-out material on the ground. Modify¬ 
ing influences would include the effects of wind, diffusion, varying 
fall velocities, decay of radioactivity during fall-out, and special 
effects such as the occurrence of precipitation and cloud and details 
of the terrain on which the fall-out would occur. 

The method must yield quantitative plots of dose rate re¬ 
ferred to H plus 1, total dose, time of arrival of fall-out and time 
of ending of fall-out. 

The method must yield results in agreement with the observa¬ 
tions made at CASTLE BRAVO. ■ 

The method must account for real winds in the atmosphere. 


The method should be applicable to 1? MT detonations, any 


day, anywhere. 
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The method must use the best available background information. 


The method must be as simple and fast in its application as 
possiliLe, consistent with not further degrading the existing sparse 
data on fall-out from high yield weapons. 
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The method should be amenable to scaling, if possible. The 
method derived below in this report was based on a 15 MT detonation but 
may be scaled either by use of the AFSWP scaling laws (which are not 
based on our cloud model but are convenient) or by the computation of 
another model of the same nature and having similar treatment but dif¬ 
fering in the values utilized. 

Sources and Acknowledgements 

Since the personnel assigned to this task had not had previous 
contact with the fall-out problem, arrangements were made to visit all 
groups active in the field which were located in or near Washington, D» 

C. Contacts were made with AFSWP, JTF-7, Hq. Air Weather Service, Office 
of the Chief of Naval Operations, Hq. Air Research and Development 
Command, and the D.S. Weather Bureau. 

Personnel interviewed included the following: 

Captain Russell Maynard 
Commander R.W. Paine 
Lt. Col. C.D. Bonnot 
Lt. Col. House 
lit. Col* C.A. Spchn 
Dr* Florence Van Straten 

Lt. Col. N.M. Luleji&n 
Dr. Lester Machta 

Source documents utilized included the following: 

AFSWP Report TAR-507 and others issued by AFSWP 
RAND Report R-265-ABC 

Instruction books of the various services and offices. 

Special acknowledgement must be given to the assistance given 
us by Captain Maynard and for the information contained in the RAND re- 5 
port. Without his extensive knowledge of sources of information and 
without the full description of the cloud growth and structure contained 
in the RAND report, it would have been quite difficult, if not impossible 
to derive the cloud model utilized in this stuc^y. 

Derivation of the Signal Corps Method DOE ARCHIVES; 

The derivation of the computation method used by the Signal Corps 
is described below in terms of a) basic concepts, b) information on the 
initial conditions within the cloud structure and related factors, c) 


Hq. AFSWP 
Hq. AFSWP 
JTF-7 
JTF-7 

Hq. Air Weather Service 
Office of the Chief of 
Naval Operations 
Hq. ARDC 

D.S. Weather Bureau 
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analysis of the basic data, d) derivation of model to obtain the dis¬ 
tribution of fall-out material on the ground, and e) the abbreviated 
method of obtaining total dose* 

Concept : 

As a result of the understanding which evolved during our inquiry, 
it was decided to attempt a simple physical model of the cloud and of 
its fall-out processes. It was conceived that the initial conditions 
of distribution of radioactivity vs. fall velocity, cloud dimensions,, 
and distribution of cloud radioactivity in height and distance could be 
specified. This cloud would then be considered as consisting of a num¬ 
ber of different "slabs" of suitable thickness* Each slab would be 
divided into "wafers" representing particles of given size classes lo¬ 
cated within the slab and each wafer would be followed to the earth 
from the initial location of the center of the slab in accordance with 
the fall velocity of the particles and the wind field. The division 
into slabs and wafers would be so done that each wafer, upon landing 
on the earth would exert "unit radioactive pressure" URP, i.e», would 
yield a given unit number of roentgens per hour referred to H plus 1 
hour over the area which it covered. The UKPs would be tallied over a 
grid system and summed to indicate the total accumulated radioactivity 
at that point. This pattern would be analyzed to yield the fall-out 
pattern* By noting the times of arrival of each wafer, its contribu¬ 
tion to total dose could also be evaluated. 

Initial Conditions ; 

It was decided to do the basic analysis in terns of a 15 MT deton¬ 
ation in order to use the realtively good size, shape, and fall-out dat# 
obtained from CASTIE BRAVO, 

It was assumed that the cloud shape would be as follows: The 
washroom top has a radius of 25 mil.es and extends from 60,000 to 95,000 
feet* The stem has a radius of 2.5 miles, extends as a cylinder to 
20,000 feet, and than flares out to mushroom diameter. In addition a 
spike may develop on the top of the mushroom cloud, but is ignored in 
our analysis since it contains bub a small volume of the total cloud. 

It was assumed that the particle size distribution and fall 
velocity curves contained in the RAND report were the best available 
and sensibly correct. This implied a relationship between radioactivity 
and fall velocity. 


DOE ARCHIVES 

It was assumed that the RAND 9.1 ratio of fall-cut activity coming 
from the mushroom tc that coming from the stem was correct. 






Having specified the initial data and assumptions vre next attempted 
to obtain as simple a mathematical description of the cloud as poosible. 

It appeared unreasonable that the radioactivity of the mushroom 
would be evenly distributed in height but, similar to the distribution 
of water vapor in turbulent air, it seemed reasonable that the mixing 
ratio of radioactivity per unit mass of air should be nearly constant 
within each of the two major regions, mushroom and stem. Uniform height 
distribution in the mushroom would have placed half the activity in the 
upper third of the mushroom mass. Ary stratification of radioactivity 
would also tend in the direction of increasing the content of the lower 
portions of the cloud. It was, therefore, decided to represent the 
activity in the cloud by a "mixing ratio model." This means that equal 
pressure thicknesses of cloud contain equal initial radioactivity, sub¬ 
ject to the differential of activity between mushroom and stem. 

The above decision led to considering all aspects of the problem, 
including cloud shape and particle fall velocities, in terms erf pressure 
as the vertical coordinate. Pressure has uniformly been related to 
height through the U.S, Standard Atmosphere. 

Representation of the flare of the stem was next considered and an 
effort made to obtain a simple continuous function which would approxi¬ 
mate the reported shape* The flare is concave in height coordinates! 
in pressure coordinates it must therefore be more so. The simple shape 
chosen was that erf an ellipse in pressure-radius coordinates. It was 
assumed that for our purposes all detonations could be assumed to occur 
at low elevations. The nominal surface pressure was taken to be 1000 
mb. An ellipse having a vertical tangent at 1000 mb and ground zero, 
and a horizontal tangent at the lower edge of the mushroom was assumed. 
The equation of this ellipse in miles and millibars is 

(1000 - p) 2 + b 2 (25 - R) 2 * c 2 



Figure 1 illustrate© this in both pressure and height coordinates. It 

will be noted that the only pertinent parameters here are the height of * 

the mushroom base and its radius. This shape was somewhat large at 

20,000 feet but, in agreement with the Rand report, assigned^ jfet&StCHlVES 

volume, and hence little activity, below 10,000 feet, DOE ARCHIVES' 

Th© volumes of tha mushroom and of the stem in pressure-radius 
space were next computed taking the mushroom as a cylinder and the stem 
as a figure of revolution bounded by tha ellipse. Using the given 
ratio of total activity led to determining that the "radioactivity mix** 
ing ratio" of the mushroom was about 12 times that of the stem, 
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mushroom 



METHOD OF DETERMINING THE STEM PROFILE 







Using the URP concept Introduced above, the mushroom was divided 
into five "slabs ' 1 of about 12 mb thickness, with the consideration that 
this number of slabs would adequately represent the differences of wind 
and of particle fall time with height. Each slab was to contain nine 
"wafers" corresponding to various particle size classes so chosen bb to 
account for all the activity within the mushroom and to adequately rep¬ 
resent the differences in fall velocity expected between various sizes 
of particles. Finally, the largest of the size classes (230 microns) 
was divided into five classes ranging up to 500 microns in order to 
better represent the fastest falling particles. The 230 micron class 
was deleted from all slabs and five wafers having an equal total activi¬ 
ty and covering the new large sizes were assigned to the middle slab of 
the mushroom. Each slab then contained 6 wafers and the central slab 
was assigned the five additional wafers, thus retaining the total of 1*5 
wafers. 

Since the stem contains primarily large particles, four large 
particle size wafers were determined and a thickness of the stem having 
the previously calculated stem radioactivity mixing ratio was computed 
which would yield unit radioactive pressure at the center of stem slabs, 
i.e, the stem was divided into UEP wafers. The volume of each such stem 
pressure slab was determined and an effective cylindrical radius, de¬ 
fining an equal volume, was computed. The stem cloud is then represented 
by a series of cylindrical slabs having varying radii. These slabs were 
about 70 mb thick and were assigned from about 70 to 280 mb. Beyond 
280 mb, double thickness slabs and double weight wafers were assigned 
down to about 700 mb. 

To obtain the particle size classes mentioned above, the Rand 
particle size data were replotted and appropriate fractional portions 
were scaled off. For the mushroom each slab contains wafers represent¬ 
ing the following nominal sizes: 20 , 1 * 2 , 68 , 85 , 100 , 110 , 125 , and 
150 microns. In addition the center slab contains particles of 180, 

190, 230, 290, and 500 microns* The stem contains wafers representing 
particles of 180, 260, 375, and 61*0 microns. Figure 2 shows how these 
size classes account for the total cloud fall-out radioactivity. 

The result of this breakdown is to yield 69 URP or double URP 
wafers totaling 81 weights in the initial cloud. To obtain a fall-out 
pattern, each of the wafers is plotted on an overlay at the point where 
the winds would move it. This process may be rather burdensome but it 
can lead to quantitative results depending upon the wind. Calibration 
to CASTLE BRATO showed that each wafer "weighed" about 100 roentgens 
per hour at H+ 1 hour. 

doe archives 

It was noted that the information on which the dimensions of the ■ 
model cloud was based came entirely from tropiesl ocean areas. The 
problem arose of how to account for varying atmospheric conditions,, In 
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the present work no account Is being taken of cloud, precipitation, dif¬ 
fusion and terrain effects. The possibilities of major changes in cloud 
height and shape with changes in the hydrostatic structure of the atmos¬ 
phere, which might occur from the test condition to perhaps polar conti- 
nental conditions, appeared important enough for some thought. Based on 
consideration of the increased buoyancy of cold air at low levels, the 
effects of entrainment of this colder air, and the relatively warm and 
low tropopaus© characteristic of polar air masses, it was estimated that 
the following table of mushroom top and base height would be expected 
as the tropopause height varied: 

Height in 

Feet of: 

Tropopause 55,000 U5,000 35,000 

Mushroom base 60,000 55,000 50,000 

Mushroom top 95,000 80,000 65,000 

A complete model computation was run on each of these structures. The 
results of this work are shown in Figure 3, which gives the structures 
in pressure terms, and in Figure U, which gives the structures in height 
terms. It, will be noted that the assumed behavior of the mushroom re¬ 
sulted in maintaining the pressure thickness of the mushroom almost 
constant, resulting in similar mixing ratio ratios between mushroom and 
stem and in identical numbers of layers having similar radii. 

Fall Time and Trajectory Calculations : 

In order to determine the point on the earth'a surface on which 
any wafer would bs deposited, it was necessary to devise a method-which 
would integrate the effect of the winds through which the particles 
foil. In principle, (0 thig required 
o'- f V jP] A n 
S ~*{>W(p,d) dp 

where V(p) is the horizontal wind at pressure p, and w(p,d) is the fall 
velocity(in millibars per hour) g depending on d, the particle size, and 
on pressure. Precis© evaluation of this Would be difficult and un¬ 
warranted by the available data. The approximation used was 

_ I n "y. 

S * W(d) |j ~K~ A p i D0E ARCHIVES 

where i indicates the reported wind level representing a pressure thiek- 
nessAPi* Jc l(p) accounts for the deviations of fall velocity from con¬ 
stancy, and w(d) reflects the fall velocity for different particle sizes. 


As indicated above, the Rand fall time curves were assumed. Thess 
.wero replotted in pressure terms for th© particle sizes used in the mod* 
el described above* All curves were found to ba concave upward, re- 
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THE THREE MODELS PLOTTED IN RADIUS - HEIGHT COORDINATES 









fleeting a greater fall velocity, in millibars per hour, at the lower 
pressures. Since this effect was about the same amount percentage-wise 
for all the smaller particle sizes and since the real effect of such a 
change in fall time was small for the larger sizes, the effect was treated 
by determining the reduction in time which the UO micron particles would 
spend in any given pressure level as a result of the curvature compared 
to the time they would spend in an equal slice near the surface. This 
gave the K(p) factors required to weight the winds representing such 
pressure levels so that U.0 micron particles passing through these levels 
at a fall velocity assumed constant and equal to that found in the lower 
atmosphere would be appropriately moved by this modified wind. Fall 
velocities of all other particle si2e classes were scaled to this UO 
micron fall velocity so as to introduce no error in the center of the 
mushroom region, yielding the W(d) factors. Thus w(p,d) was treated 
as 

W(p,d)= Klp)W(d) 

A continuous estimate of the wind as a function of pressure V(p), 
is not available to the meteorologist under normal conditions. In 
practice, he is given a number of winds representing the air velocity 
at a number of given heights, commonly based upon averages over time 
intervals of the order of one minute (i.a. average over height intervals 
of the order of 1000 feet). Climatological data may not include winds 
except for the "mandatory levels,” i.e. 1000, 850, 700, 500, and 300 
millibars. In using wind data to compute trajectories, the winds are 
weighted in proportion to the pressure thickness for which they are 
(and must be) taken as representative. For example, the 850 mb wind 
may be used to represent the atmospheric flow between 925 and 775 milli¬ 
bars and is weighted by 1.5 since it represents I50 millibars in our 
vertical coordinate. 

Thus wind data are weighted to account for the deviation of parti¬ 
cle fall velocity from constancy and to account for the thickness of 
the atmosphere which they represent. A single weighting function table 
is used. True (i.e. Rand) fall times are used, however, for determin¬ 
ing the time of first and last arrival of wafers, since the fall time 
for each particle size from each of the nominal pressure levels is "" 
easily tabulated. 

Pose Rate and Total Dose: 

Since the study for which fall-out computations were required had 
a number of facets, it was necessary to determine both the dose rate 
referred to Htl hours, D, and the total dose, Dp, received up to an 
arbitrary time, T, later than the ending of fall-out. POE ARCHIVES 

The position of each wafer as it arrived at the surface was noted 
and the numoer of UR? wafers lying over any required point was noted. 
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For CASTLE BRAVO it was found that each wafer represented about 100 r/hr 
referred to H + l. This calibration was used for all other calculations. 


since it was the only good point of contact between wafers on the sur¬ 
face and observed fall-out. Following tbe apDroach of AFSWP Report $07 
we have assumed a decay rate given by t"^*^ for all times of interest, 


i 


although deviations from this which would be the same for all wafers 


could be accounted for without producing significant complication of 


the computational problem. 


The reauirement for some sort of integrated dose related to human 
effects at first appeared to be exceedingly difficult to meet. Inquiry 
relating to the present knowledge of biological recovery, however, re¬ 
vealed that (a) the facts are not too well known, particularly as thqy 
relate to human beings, and (b) the recovery within the first few days 
following beginning of exposure was estimated to be less than 30&> With¬ 
in the accuracy of our computation method and considering the uncertain¬ 
ties of the recovery problem, it was decided not to allow for recovery 
but rather to compute the total dose. Since each URP wafer is followed 
to the surface and its time of arrival determinable, one could sum the 
total dose which each delivered up to any given time. This, however, 
with some 69 wafers per problem, would be rather lengthy. It was felt 
that some mathematically simplifying assumption regarding the rate of 
accumulation of radioactivity on the ground would be desirable. The 
assumption made was that all the fall-out arrived at a constant rats 
from the time of first arrival, a, and the time of ending, e. Assuming 
total fall-out to be 1 r/hr at H+l, the total dose to time T, where 
T>e, is given by the following simple formula: 


Q = 6.25-— 

i e — a 


— 5 T 


The first term corresponds to the infinity dose delivered by unit dose 
rate accrued between a and e and the second to the dose which would b® 
delivered between T and infinity. In use, the dose rate V is multi¬ 
plied by the factor given by the formula above, or determined by 
entering a chart with a and e, to obtain the total dose to T. Dp so 
obtained was plotted against coordinates a and e and analyzed as a 
contour chart. 



of the Model and Method 


DOE ARCHIVES* 


In the practical aoDlication of the Signal Corps method devised 
above, ail effective wind trajectory was constructed by connecting the 
weighted winds corresponding to each interval head to tail. The wind 
for any given slab was picked off by interpolation along the length of 
the appropriate vector. These winds were actually plotted on a scale 
such that the 68 micron particles fell directly along the plotted 
trajectory. All other particle size destinations, as explained above, 
then scaled linearly with this wind. In order to simplify the practi- 







cal work, a triangle was constructed similar to that described earlier 
by Dr. Van Straten, except that it accounted for a,larger number of 
particle sizes. When the length of the 68 micron trajectory is laid 
parallel to the base of the scale, with the origin at the left-hand 
edge and the end of the vector intersecting the 68 micron line, the 
destination of all other sizes may be plotted at the intersection of the 
appropriate sloping lines with the 68 micron vector. Thus, all eight 
(or thirteen) wafers corresponding to a slab within the mushroom can be 
plotted with ease, and similarly, the four wafers representing each stem 
slab can be rapidly laid out. A "poker chip" having a diameter corres¬ 
ponding to the proper radius for the slab which it represents (i.e. on 
the same scale as the 68 micron trajectory plot) is located at each UR? 
wafer destination point. An overlay, ruled in a suitable grid# is 
positioned and those grid sections which are at least 50% covered by a 
poker chip are marked. This is done for each of the eleven slabs cor¬ 
responding to 11 effective winds, using the same overlay. When all the 
poker chips have been tallied, the number of URP's lying over each area 
are summed and the resulting pattern analyzed as a pressure field. Since 
each URP was determined to represent 100 r/hr at H+l, the fall-out field 
is.readily obtained. 

At the same time that the "poker chips" are tallied in the overlay, 
the leading and trailing edge of the chip are sketched and marked with 
the appropriate time of arrival. From the field of these values, con¬ 
tours of earliest time of arrival, designated above as a, and of latest 
time of arrival, designated above as the time of ending, e, are drawn. 

Using the assumption that fall-out (in units of r/hr at H-f 1) accumu¬ 
lates linearly between a and e, and the equation for total dose dis¬ 
cussed above, total dose to T = a + U3 is computed. In our later 
practice, a contour chart entered in a and e and having contours label¬ 
ed in terms of total dose to a + U3 per unit dose rate is used* The 
values so obtained are again analyzed to obtain the total dose contours. 

Since all three models have the same number of slabs and wafers, ■ 
there is no modification of the computation method necessary except for 
the location of the nominal initial pressure of each slab. Having 
selected the model, the computation procedures may be summarized as 
follows: (1) The weighted winds to the top of the model are plotted 
in an effective trajectory, (2) the 68 micron trajectory vectors are 
drawn for each slab, (3) the wafers are positioned using the scaling 
triangle, (it) URP contributions are tallied on the overlay along with 
a and e, (5) the sum of tallies and a and e are contour analyzed, (6) 
total dose is evaluated for a suitable number of points and (7) a total 
dose chart Drepared. 

__ ‘ DOE ARCHIVES 

«~Xn practical work a 15 mile sauare grid was used, while for Nr. Barnett's 
paper a 3.5 mile square grid was used. 



A paper of this nature would not be complete without reference to 
those assumptions which are considered weak and consideration of pos¬ 
sible means of improving the method. 

We aknesses ; * 

Weaknesses in the method stem from both inadequate basic data and 
from deficiencies in the development of the method. In the former class 
the major deficiencies lie in the fact that useful data on total fall¬ 
out, original distribution of radioactivity for successive fractions of 
this radioactivity (or, alternatively, distribution of activity by 
particle size and height and fall velocity), biological recovery effects, 
decay of radioactivity, diffusion, and terrain effects are either sparse 
or non-existent. The major deficiencies in the model lie in the length 
of time required to perform a computation, the linear arrival assump¬ 
tion, and the finite number of wafers considered. It will be noted that, 
without the use of machine computation methods, an increase in the num¬ 
ber of wafers (which would yield somewhat smoother patterns) would further 
aggravate the time problem. 

It is assumed that many of the informational deficiencies will be 
met by future tests and that more certain estimates of those which can¬ 
not be completely met will result from the analysis of the data obtained,, 

With regard to the model and method, a number of improvements appear 
possible: 

► 

Increase the number of slabs and wafers per slab—this 
would produce smoother patterns. 

Perfect an alternate method of obtaining a and e based on the 
determination of the maximum and minimum effective winds in each azimuth, 

or 

Replace the linear accretion assumption with immediate sum- ..." 
mation of total contribution from arrival time to a + lj.3 hours* „ 

Adapt the method to machine computation methods. 

Use time-varying trajectories. DOE ARCHIVER 

Use a better approximation of stem shape. 

Incorporate any new test results to refine the model and 

method. 
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Provide an improved method of scaling either through alter¬ 
nate models or otherwise. 

Modify model to take account cf diffusion, especially if 
smaller yield weapons are to be considered. 

Use an overlay grid ruled in coordinates suitable to the 

problem. 

Account for variations in particle size distributions of 
various soils and targets. 

Incorporate biological recovery corrections. 

Modify method so as to reduce computation time while retain¬ 
ing quantitative results of reasonable accuracy. 

&OE ARCHIVER 
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QUESTIONS AND ANSWERS FOLLOWING DR. SWINGLE’S PRESENTATION 


Question: LTCOL Spohn, AWS 

On the last paper, I am curious of the time interval it 
takes to actually work through a plot - how long does it 
take to run a plot? 

Answer: Dr. Swingle, OCSO 

After some experience in a well equipped office, three or 
four man hours. Most of the work we have done was not for 
operational use; we set it up on a production line basis 
and at one time we had as many as 20 men working on it. 

Question: FCDA 

What typical times between arrival of the smaller particles 
and larger particles; that is, in your computation to get 
dosage figure? 

Answer: Dr. Swingle, OCSO 

I think you will see that in-Mr. Barnett's paper tomorrow, 
but it runs for something like half an hour to a matter of 
ten hours, or more. 

Question: LTCOL Spohn, AWS 

In applying wind data to this, do you use the same assumption 
as we (AWS) do; that is, assinning point winds apply throughout 
the time of fall-out? 



Answer: 


Dr. Swingle, OCSO 


DOE ARCHlV Ef 


Yes, so far. In the study which we have done we have been 
concerned with getting what fall-out might be. A single 
sounding might be just as good as changing it. We have not 
had an occasion to want to. 
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PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS FUNCTIONS OF YIELD 
AND METEOROLOGICAL CONDITIONS, ARDC METHOD 

LtCol N. M. Lulejian, USAF 
Air Research and Development Command 


I have analyzed the fall-out from the first shot of CASTLE test 
operation. The method of analysis is based on Stokes Law of fall and 
upon the following assumptions of particle size distribution and 
activity with height within the atomic cloud: 
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vhere h - Height above target in thousands of feet (or height in 
cloud) 

A = Percentage activity in each 5000 ft slice,of the cloud. 

£A is the cumulative percentage activity. 

djnean ■ Mean particle diameter in microns. 

I want to now refer to my earlier talk on the percentage fall-out 
as a function of scaled height. During this discourse I showed that 
as the scaled height was reduced the percentage activity in the stem 
of the cloud increased as compared to the activity in the mushroom,. 
Armed with this fact and also with the analysis of particle size 
in the tower shots of TIMSLER-SNAPPER and UPSHOT-KNOTHOLE test opera¬ 
tions, X arrived at the assumptions made in paragraph 1 above. For 
greater details in this regard, refer to ARBC Report No. C4-23 676 . 


For the winds we have used the composite winds of Eniwetok, 

Rongerik, and Bikini for H - hour and H + 2:15 hours. For H + 8 and 
H * l4 hours we have used the composite winds of Eniwetok and Rongerik 
weather stations. Figure 2 shows the time variation of the winds, and 
figures 3 and 4 show the composite wind analysis for different sized 
particles. Figure 1 gives the names of the islands of interest in 
the Pacific. ’ 

DOE ARCHIVES 


The main point of difference between our method of analysis and 
that of RAND is that their assumptions of particle size and activity 
distribution within the cloud are radically different from ours. I 
believe this is o nly possible because the stratospheric winds had a 


270 













Figure 3 





Figure It 



directional shear of approximately l 80 ° as compared to tropospheric 
winds during Bravo shot. I say that the fall-out occurred from the 
stem; RAND says it occurred primarily from the mushroom; hut since 
the ground trajectory of the stem and mushroom actually coincide in 
many pointB (see figure 2), the problem cannot be resolved until at 
some future test when the wind shears between the troposphere and 
the stratosphere are at approximately 90 ° • When and if this happens 
at some future test, we will he able to resolve the problem uniquely, 
provided we are there to catch the fall-out. 

I want to make it clear that my assumptions (given in paragraph l) 
apply only to a surface burst on land of an atomic bomb In the range 
of 10 to 20 MI. It does not apply to a surface burst over water 
(barge shot). I don't know anything about particle size distribution 
or activity distribution in a r Vater ?1 cloud. All I know is that 
former speakers of this conference (ERDL, Scripps and New York Opera¬ 
tions Office of the AEG) have indicated that the total fall-out may be 
2 to 5$ (Mr. Eisenbud) or it may be 13 to 58 $ (NRDL-Scripps). All 
speakers agreed that they make no claim to have caught all of the fall¬ 
out. So, as far as I know, fall-out from water shots may he 2 to 58 $ 
or more. The activity within a water cloud may be anything. I don't 
even have a clue as to what it, may be. 


I have used Stokes Law of fall as compared to aerodynamic fall. 

The difference between the two rates of fall is relatively small in 
the size range of 60 to 100 microns, and it is in this very size range 
that we find most of the activity in the cloud. Therefore, for ARDC 
method, the use of Stokes law or aerodynamic fall makes relatively 
little difference. However, for those other units who assume that 
particle sizes between 100 and 1000 microns hold most of the cloud 
activity, the use of aerodynamic fall becomes mandatory. It is my 
claim that the primary difference in the methods presented is the 
particle size distribution with activity and height, not the different 

rates of fall. DOE ARCHIVES* 

Figures 5 and 6 show the measured dose rates on the islands in the 
Pacific for the first shot of Operation CASTLE. These two figures are 
taken from CASTLE Project 2.5a preliminary report by NRDL. They are 
shown here as the data that we must try to approximate with our plot 
of fall-out. Figure 7 is the calculated fall-out area using the 
assumptions made in paragraph 1. Figure 8 shows the ARDC fall-out 
plot in isodose contour lines of dosage integrated from time of fall¬ 
out to 48 hours after the bomb detonation, using the t"^*^ relation. 

It should be noted that all points on figure 7 are calculated points. 
They do not show any measured data. All the measured data is- in 
figures 5 and 6. 
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Now let me talk to you about what I call the "Volume-Effect" of 
fall-out which was found during TUMBLER-SNAPPER and UPSHOT-KNOTHOLE 
test operations in Nevada. In pages 65 to 67 I discuss this matter 
in greater detail in ARDC Report No. c4— 23676 . Figure 9 » which is 
taken from AFSWP Report WT-558, shows this "Volume-Effect" at Lincoln 
Mine (approximately 40 miles north of ground zero) during TUMBLER- 
SNAPPER Shot No. 5 . Similar "Volume-Effects" of fall-out were 
obtained during TUMBLER-SNAPPER and UPSHOT-KNOTHOLE tower shots. 

Figure 9 shows that the 48 hour integrated dose calculated by 
t~i- 2 relation must be multiplied by I .67 to account for the "Volume- 
Effect" of fall-out. In ARDC Report C4-23676, we have indicated what 
would happen if thiB same "Volume-Effect" is taken into account for 
the CASTLE Bravo fall-out. During the evening sessions of this 
symposium I have been advised by AFSWP and NRDL personnel that the 
"Volume-Effect" was not observed in and around the Bikini Lagoon. I 
have been told that there were six continuously recording stations 
on the islands of Bikini Atoll and that these traces did not show the 
"Volume-Effect" of the type I found in the desert. It was also con¬ 
jectured by AFSWP personnel (in night session) that if this "Volume- 
Effect" occurs at all in fall-out from multi-megaton weapons* it 
would be observed 500 to 1000 miles downwind (where the total fall¬ 
out would be small). It is my opinion that this "Volume-Effect" 
should have appeared on the northern portions of Rongelap and 
Rongerik quite strongly. It should also have appeared to a lesser 
extent on the southern sectors of these two atolls. The fact that 
they were not observed at Rongerik may be due to the fact that there 
was very little measured data downwind during the fall-out period. I 
believe that the distinct possibility exists that Rongelap natives 
may have received from 1.2 to 1.5 times the 48 hour integrated dose 
used by the investigators in this field of endeavor. I don't know 
whether this same "Volume-Effect"would occur from fall-out coming 
from bombs detonated not on the ground* but on water. As I said 
earlier, I know nothing about fall-out from "water clouds". How¬ 
ever, I believe that this "Volume-Effect" will be observed in future . 
tests for surface burst weapons because I believe these things scale 
from Nevada to the Pacific. -, ; 

t)OE ARCHIVED 

It is believed that the ARDC method of fall-out analysis ade¬ 
quately explains the fall-out picture shown in Figures 5 and 6. The 
important points that must be kept in mind when an attempt is made 
to fit assumed particle size and activity data to the actual fall-out 
case are the following: 


The time variation of the winds was so great that composite 
analysis of the winds is required. The winds at ground zero are not 
valid over the downwind path after the lapse of only one or two hours. 
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Aliy analysis of the fall-out on Bikini Lagoon must take into 
account the amazing fact that the intensity of fall-out on the lagoon 
islands to the south and southwest of ground zero is only 8 to 15 
r/hr at H+l hour, even for islands as close as approximately 5 nauti¬ 
cal mil es from ground zero. On the other islands of Bikini Lagoon, 
the contamination goes to 1000 to 2100 r/hr at H+l hour (see figures 5 * 

and 6 ). 

The atoll of Rongelap is contaminated excessively in a 
steep gradient as shown in figure 5 . 

The native village on Rongelap Atoll (southern portion of 
atoll) is in approximately the same radiation intensity field as the 
weather station on Rongerik. 

Ailinginae shows less contamination than Rongelap or Ronge¬ 
rik. 

Bikar has less contamination than Rongerik. 

Task Force ships 50 miles southeast of ground zero showed 
less contamination than Ailinginae. 

The islands of Taongi and Utirik were least contaminated. 

It was later found that the position of the Japanese fishing 
"boat (Fulcuru Maru) was not 10° north of east (80 miles) of ground * 

zero, hut thanks to Mr. Eisenbud of MYOQ, AEG, the exact position is 
shown in figures 1 and 8 as point X and the direction of motion of 
the ship is shown to he from point X to point A. Extrapolation of 
the Eisenbud data on the Japanese ship shows that it reached the 12° 

H. latitude line at approximately E +6 to E *8 hours. The fall-out was 
supposed to have stopped at h 48 hours according to Japanese accounts. 


We have kept in mind the measured data mentioned in the paragraph 
above, and our model fits the above data quite well (see figures 7 and 
8 ). The winds are those obtained from the shot island (Bikini), 
Rongerik and Ed wet ok. Looking over the papers presented by different 
members of this symposium I note the following: 

The RAM) Aureole report shows CASTLE Bravo analyst 
overestimates the fall -out around ground zero (especially on the 
southerly islands of the atoll). The downwind fall-out misses Ronge¬ 
rik and Rongelap and if memory serves me right, it also misses Bikar. 
This model then does not succeed in approximating the measured Bravo 
fall-out. At this symposium, RAND has presented new and as yet 
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unpublished data which adequately accounts for the fall-out in and 
around the islands of Bikini Atoll. It also appears to give the 
correct order of magnitude of the contamination on the northern por¬ 
tion of Rongelap Atoll. The analysis does not go any further than 
that as far as I can tell. 

The KRDL analysis of Bravo event was widened too far to get 
the fall-out to fall on Rongelap and Rongerik. It is believed that 
since KRDL used only the ground zero winds (not taking the time varia¬ 
tion of winds into account) then they were forced to widen their area 
of fall-out to account for the observed fall-out. 

The particle size and the activity distribution in the cloud 
shown in paragraph 1 above, was obtained after careful analysis of 
the TUMBLER-MAPPER and UPSHOT-KMOTHOLE tower shots and the JANGLE 
surface shot. Reference is made to ARDC Report 03-364-17* Cta- page 4-0 
of this report is shown that the following particle size distribution 
fits the UPSHOT-KNOTHOLE tower shots (17 to 50 KT from 300 foot 
towers): 

125 microns in lower third of stem 

100 microns in middle third of stem 
90 microns in upper third of stem 

80 to 70 microns in mushroom of cloud for a cloud reaching 
a maximum altitude of 40,000 feet msl 

It is assumed that as the cloud height increases the particle sice 
distribution at the higher levels in the cloud decreases. It is 'from 
this background of detailed analysis of the tower shots of UP SHOT- 
KNOTHOLE and TUMBLER-SNAPPER that the assumptions of paragraph 1 
were made. So you see, we didn't pick them out of thin air after 
alii It is unfortunate that most other units have neglected analyses 
of the Nevada shots as a guide to Bravo fall-out. 

The next point I want to discuss is my "reasoning" for assuming 
that the activity in the stem is greater than the activity in the 
mushroom. During TUMBLER-SNAPPER tower shots (12 to l4 KT from 
300 ft towers) the activity in the stem was approximately the same 
as the activity in the mushroom. However, during the UPSHOT-KNOTHOLE 
test operation (17 KT to 50 KT from 300 ft towers) the activity in 
the stem was 2 to 3 times the percentage activity in the mushroom. 

For the JANGLE surface shot the activity in the stem was approximately 
thrice times that found, in the mushroom. It is for this reason that 
we assumed in paragraph 1, that approximately Sofa of the activity 
is in the stem and only in the mushroom. This appears to be 
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contrary to common sense, but I can assure you I*m not doing this to 
be ornery. I am making this assumption because I believe that we can 
reasonably scale from tower shots to surface shots and we can also 
scale for increasing yield. 

We employ the following scaling law: 

Pi 5 Pg for a given scaled height, 
where 

P - EAR Equation 1 

g^t=1.2 

P = Percentage fall-out of residual activity 
A * Area in square miles 
R = Dose rate at time of fall-out in r/hr 
W = Bomb yield in kllotons 
t - Time in hours after bomb detonation 


k B Constant * 12 in our calculations 


This means that we assume that the percentage fall-out from a given 
scaled height is constant. This refers to residual activity deposited 
on soil debris which falls out within approximately one day due to the 
gravity of the particle. Equation 1 may be readily written as: 




Equation 2 


Then for any element of Area, A^, inclosed by a contour line R^, we 
may be permitted to write: DOE ARCHIVES 

t = t r^irAi ' 2 

2 *■ W.JL f«J Equation 3 


We can also say: 

PT o' f W 2l f f ' 1 1,2 Equation V 

R *~ R » L w[J L"£j 

We note above that the area scaling with yield depends upon the time. 
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In the normal formula for scaling it is assumed that 

* * T W z -) z /3 ' Equation 5 

V A, l-j^J 

/ ! 

The time variable from Equation 5/may be eliminated by setting 
t ~ Z ~ f where Z is distance in the vertical direction. 

If we do this we get 




Equation 6 


It may be argued that Equation 6 and Equation 5 are not too different; 
and this may be true, except that we have a problem of extrapolating 
from 1 KT to 15,000 KT, where the use of the 0.60 exponent produces 
a factor of two differences from the use of the O.667 exponent. The 
reason we have gone through the above analysis is to inpress the 
reader with the fact that the time of start of fall-out is very impor¬ 
tant in any scaling factor for integrated dose. Paradoxically there 
is least amount of experimental data on this. 


Once we have the Bravo model and ve try to extrapolate back to 
1 KT, we find that we overestimate the problem by a factor of 3 or 4 
using Equation 5 and by 6 to 8 when Equation 6 is used. We have 
applied a correction factor to compensate for this. We believe that 
by such a procedure we may be able to use the Bravo shot to scale 
down to 1 MT and up to 60 MT quite adequately. For details see 

ARDC Report C4-23676. DOE ARCHIVE! 

During the conference a number of personnel referred to the ARDC 
method as implying 90$ activity in the stem and 10$ activity in the 
mushroom. Actually in our report (c4-23676) and in the first para¬ 
graph above, we showed that our method assumes approximately 80$ 
activity in the stem and 20$ in the mushroom, whereas RAND assumed 
10$ in the stem and 90$ in the mushroom. We believe that large 
particles rise high into the mushroom, immediately as the cloud forms, 
but it is our opinion that all this falls out fast (within 10 to 30 
minutes) and it all lands on or near ground zero or the lagoon. Ho 
significant amount of activity falls out downwind in excess of 20$ 
from the mushroom except after 24 hours or so, at which time both 
dilution and decay would have reduced the radiation intensities to 
insignificant amounts. For lower yield bombs such as 5 MT, the 
activity in the mushroom my increase to 30 or 40$, with the stem 
activity still as high as 60 to 70$ • This is because we assume that 




the particle size distribution decreases with height, and the 5 
cloud would not reach as high as the 15 MT cloud of the Bravo shot. 

Figure 10 shows what would happen to this country if 111 bombs 
of 15 HE are detonated over 106 cities whose population is 100,000 
or more and over five other selected areas. For details refer to r 

ARDC Report CU-23676. The slide is here shown to illustrate that 
if 100 or more hydrogen bombs are used, most of the country will be 
contaminated almost independently of meteorology. For example, if 
the winds aloft were all easterly, it would not radically change the 
contamination pattern over the country. Similarly for northerly or 
southerly winds. The slide also shows that evacuation and dispersal 
may not be advisable, for it may clearly be a case of evacuating from 
the frying pan into the fire. 

X would like to make a few rather cbvious remarks about world¬ 
wide contamination that say be derived from our analysis of percentage 
fall-out for surface burst bombs. It will be assumed that 80$ of the 
residual activity of a bomb is deposited downwind within 15 to 20 
hours after bomb detonation. This means that if we need 10,000 bombs, 
airburst, of a certain yield to uniformly contaminate the world, then 
we will need only approximately 200 of the same yield bombs to ccnt&mi'- 
nate this country to the same level. 

This presentation of the AllDC method of prediction of dose rate 
and dosage contours was not presented in the manner indicated above. 

However, the ARDC method has not been altered in this presentation. ' 

This paper represents rather the many "conferences" within this con¬ 
ference that were held during these three days. Some of the items, 
notably the scaling method was not presented at the symposium due to 
shortage of time, but it is included here for sake of completeness. 

This paper is in effect a condensation of ARDC Report C3-36V17 and 

its supplement, and ARDC Report C4-23676. The illustrations shown 

here are from the above reports. This paper in substance indicates 

not only my talks at this symposium, but also my thoughts by the ./ 

close of the conference after X had seen the presentation by other ^ 

members of the symposium. My final comment is that it is unfortunate 

that RAM) and other organizations did not utilize the wealth of 

information available from analysis of TUMBLER-SHAPPER and UPSHOT- 

KEOTHOLE tower shots. Had they done so, they would not have placed 

bo much scavengable activity in the mushroom. v 
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SURVEY OF FALIDUT AREA AT SEA FOLLOWING 
YANKEE - NECTAR, OPERATION CASTLE 

Mr. Iferril Eisenbud 
U. S. Atomic Energy Commission 
New York Operations Office 
Health and Safety Laboratory 


The fallout following Yankee and Nectar was estimated by means of 
aerial surveys. Evaluation of the data indicated that in the region 
of from 25 to 150 miles the estimated Yankee fallout of fission pro¬ 
duct, adjusted to Hfl days, was 12ljQ megacuries, the estimated Nectar 
foil nut , at HU days, was 293 megac uries. 

SETCF'TtrTs^po ssilW to measure surface activity “by aircmftT"' 
It is necessary only to know the law of height attenuation for the 
characteristic gamma radiations measured. In the case of fission pro¬ 
duct the height attenuation is shown in figure 1. Figure 2 illustrates 
also Idle correlation between the surface gamma activity as estimated 
from the aircraft, and measurements taken on shipboard. The estimate 
of total activity could be made since there was also correlation be¬ 
tween surface activity levels and total fission product load in depth 
to the bottom of the thermocline, figure 3. 

Figure U is the trace of activity taken on H+72 hours following 
Yankee, figure 5 is the related diagram of wind vectors, figure 6 re¬ 
presents the estimated isodose traces. Figures 7, 8,9, 10, and 11 
bear the same relationship to Nectar. It is interesting to note that 
the wind vector charts effectively display the probable path of fall¬ 
out. In the case of Nectar the low level fallout, apparently stem 
debris, was found on H4-21+ hours but disappeared during the next 2U 
hour period. 

The method, beyond providing information on quantity of fallout, 
supplies data which will permit prediction of the direction of fallout 
paths, a means of determining safe areas for surface vessel oper¬ 
ations and a general source of information for military and civilian 
defense planning for personnel protection. 

Future planning should include the use of telemetering equipment 
plus automatic altitude compensation of readings to eliminate pilot 
error. The use of telemetering will permit -a central plotting arrange¬ 
ment and is necessary since more than one plane will be required to 
provide a traverse of a 360° pattern. Sea water surface measurements 
should be related to depth measurements from the time immediately 
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following the burst until ths main body of activity disappears since 
no early data has been taken. With better equipment and adequate 
preparation it may be desirable also to survey beyond the 1?0 mile 
radius. 
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PREDICTION OF DOSE-RATE AND DOSAGE OONTODftS AS 
FUNCTIONS OF HELD AND METEOROLOGICAL ■ CONDITION, 

LASL METHOD 

PART I 

FORECASTING OF THE 10 R ISODOSE LINE 
Dr. Gaelen Felt 

Los Alamos Scientific Laboratory- 


Introduction 


The problem of long-range fall-out from very large devices 
(megatons) was first examined prior to Operation IVY, At that time 
the methods of analysis were based on a simple theoretical model de¬ 
vised to evaluate the hazards from the Jangle shots and on the 
empirical results from low-yield tests In Nevada* Neither source was 
truly applicable In detail to Mlke,^^^&pell-known at the time, but 
the general result, as we know frcmJpHHp was correct; i.e., fall-out 
from Mike under adverse condi tions|lrouidnave been very severe at dis¬ 
tances of 200 to 300 miles* The actual conditions on Mike-Day were, of 
course, favorable but at the same time rather unusual in that the lo¬ 
cation of the main fall-out was well clear of all populated areas* 

The rare occurrence of ideal conditions, the length of the Castle 
(operation, and the very evident hazards from the devices made necessary 
a re-examination of the problem of fall-out* In the field such a re¬ 
examination could be very crude at best, geared as it was to the 
Immediately practical aim of operational forecasting. It was necessary 
to devise a system of forecasting simple enough to provide results 
based on the latest possible weather information, both observations and 
forecasts, conservat ive e^g gh to guarantee no repetition of the un¬ 
fortunate results ofHBljl|rand yet daring enough to enable one to 
take advantage of w«ftne^conditdons different from those known to be 
ideal (no northern components from surface to 100,000 ft and axis of 
fall-out between about 31?° and 0L5°). DOE ARCHIVES 

The treatment generally used in Nevada with moderate success was 
pretty clearly not well suited to fall-out forecasting at the Pacific 
Proving Grounds except possible under the condition of "deep easter¬ 
lies", when all wind vectors from the surface up are easterly (an 
alternate "ideal" situation, incidentally for Eniwetok Atoll but not 
for Bikini)* The basic difficulty with using the Nevada system for 




analysis of the effects of large devices in the Pacific is that the 
wind structure in the Pacific is not primarily unidirectional in the 
pertinent altitude range. The usual fall-out pattern in Nevada is a 
narrow band, broadened at times by the presence of light and variable 
winds and by the occasional presence of abnormally great directional 
shear but st ill confined to a snail sector. This general character is 
not often found in the Pacific, where the directional shear is almost 
always very great just above the east-northeast flow in the low al¬ 
titudes and again at the tropopausa (usually about 55,OCX) to 60,000 ft). 
In the Pacific the sector which includes all the distant fall-out is 
very frequently 180° and occasionally even greater* 

While the sector is very large in the Pacific it is not true that 
the fall-out is uniformly distributed within it, and it therefore be« 
comes necessary to look closely at the angular distribution of fall¬ 
out as well as the radial* It is this new factor in the analysis which 
makes the Castle system differ fundamentally from the Nevada system* 

It does not necessarily follow that tho Nevada system should not also 
have this feature since there are evident failings there which may 
perhaps be removed in this way* 

The Dose Index 

The first step in the new system of analysis was to calculate the 
relative dose at various points about the origin for a given wind 
structure. The mechanics of preparing the isodose plots will be de¬ 
scribed in the following sections. Here we shall discuss the assump¬ 
tions and general principles of the net hod* The assumptions used in 
the calculations are listed below* 

The activity is uniformly distributed in height. This assumption is 
obviously a poor one and was subsequently altered to emphasize the 
middle region of the cloud and to depress the stem and the top. In 
the first plot this modification was applied at the plotting stage 
rather than at the calculations! stage* 

The distribution of particle sizes is uniform. This 
is not true but is justified by several arguments, 
view—distant fallout—the number of particles two feet in diameter 
relative to the number 100 microns in diameter makes not the slightest 
difference since the former all fall in the region of the crater any¬ 
way. On the other hand the relative number of very small particles is 
equally uninteresting since these do not fall quickly enough to affect 
the results significantly. What we have done in effect is to assume a 
’’white” distribution for simplicity and to restrict the region of 
validity to distances between approximately fifteen miles and a couple 
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From our point of 



of hundred miles from the origin. Again, in the absence of any de¬ 
tailed information on the real distribution, the assumption we hare 
made is probably the safest. 

The amount of activity deposited by an active particle is proportional 
to its area. One may fairly argue that in some cases the activity 
should be proportional to the particle volume. The choice of area 
dependence is based on the idea of plating on molten solids and on the 
idea of scavenging by water droplets. 

Th| jjotal activity is proportional to the total yield and decays as 

The particles settle by Stoke's law. This assumption is always 
challenged but no suitable substitute has yet been suggested. In the 
range of particle sizes of interest it is probably as good a law as 
any modification of it. 

The area at the surface covered by particles falling from a given 
height is proportional to the square of the time of fall. This 
assumption expresses the fact of divergence which may arise from any 
number of reasons. Mo numerical value is assigned to the divergence. 

The net radial distance that a particle travels is proportional to the 
time of fall. 


With the above assumptions one obtains the dose index in the 
following inelegant way. At a given point "P" on the surface (Fig. 1) 
active particles can arrive from certain discrete altitudes of the 
cloud determined by the intersection with the he dograph of the radius 
vector from ground zerc through the point. In the example these are 
the altitudes 79.000 ft at "A" and 39,000 ft at "B". The hodograph is 
customarily drawn for a fall rate of 5,000 ft per hour so that particles 
with that fall rate which start from 39,000 ft at nominally zero time 
will land at "B" in 7.3 hours. These particles will have a definite 
diameter "d". We have said that the activity brought down by 


any particle is proportional to its area; i.e., 
contains a factor n d^ rt . 


the dose index "D n 
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The remainder of the argument concerns the time. Two other 
assumptions contain time terms, the natural decay and the area of 
deposition. Now the dose rate is expressible in activity per unit 
area. For example, one megacurie per square mile is roughly equiva¬ 
lent to 3 r/hr about 3 ft above a fission fragment deposition. The 
dose rate is therefore expressible as 
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rate 


a ^ t 
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and the integrated dose, apart from other factors, is expressible as 


dose 


l-2.2 


We chose to define the dose index' as 

a 2 


D 


( 2 ) 


(3) 


where "d" is the particle diameter in microns and "t" is the time at 
which it arrives in hours. In the event that particles arrive from 
two heights as indicated for point "P" in Fig. 1, one merely adds the 
two indices arithmetically. We have felt, in view of the crudeness of 
these arguments, that very little would be achieved by keeping the 
power of "t n at -2.2 and have used -2.0 instead. 

One will note that in this form the dose index does not yet con¬ 
tain Stoke’s law but states merely that particles of size n d” arrive 
at tii® tt t". The use of Stoke’s law and some of the other assumptions 
permits the dose index to be written in a wide variety of alternate 
forms. Stoke's law states that the terminal rate of fall is pro¬ 
portional to the area of the particle. In our case the terminal rate 
is reached so quickly that one may write 


h 


K d* 


U) 


where tt h H is the starting height, n t n is the time of fall, n d n is the 
particle diameter and "K” is a constant of proportionality containing 
the density, the viscosity, and numerical factors. The viscosity, in¬ 
cidentally, is temperature dependent, but the variation ie small com¬ 
pered to the other uncertainties in this system and we have chosen to 
keep the viscosity constant. Substitution of equation (U) permits one 
to writ© the dose index as 


h . 1 
K T3 


and further use of assumption (7) above permits the form 


D ** h .1 
K- "3 


(3a) 
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where "K‘ * is some new constant and "R n is the radial distance along 
a bearing from ground zero. 

Numerical values of the dose index can be computed from any of the 
forms of equation (3) with the proper selection of unite. 

We defined the numerical value of the dose index as the ratio of 
the square of the particle diameter in microns to the square of the 
time of fall in hours. In these units the dose index is of the same 
order as the dose in roentgens from a 10-megaton yield as determined 
by a rough match ^thflHHMfaata. The adjustment to other yields 
in the megaton range SgS&ae by direct proportionality as indicated 
in assumption (li). ” 


ration of the Plotr Static Case 


We proceed now to describe the method of constructing the isodose 
contours for the simpler case of the static hodograph. The static 
hodograph presents a plan view of the vertical structure of the fore¬ 
cast winds at zero time and at the origin. Fig. 1 is a sample hodo¬ 
graph constructed from data in Table I. The criticisms against the 
static hodograph are that it should not be expected to persist un¬ 
changed even at the origin throughout the time of fall-out, and that 
it does not in principle apply at all to points displaced from the 
origin. These criticisms become more important as one attempts to 
forecast the low-level isodose lines which are established at con¬ 
siderable distances from the origin and many hours after shot time. 
The static hodograph does, however, provide a useful guide and has 
the virtue that a plot can be prepared from it very quickly indeed. 


Using the sample hodograph of Fig. 1 let us compute the dose in¬ 
dex for point "P" at h0° and 60 nautical miles from ground zero. There 
will be two components for "P" and all other points along the radius 
from "Z n through "P*. from the hodograph the intercept heights and 
distances are 79,000 ft and i*5 miles at B A", and 39,000 ft and 73 
miles at B B". As mentioned above the total index for w P n is the sum 
of the indices computed from the two intercepts. 

Let us consider first the index from intercept n A". 
which fall at "A" on the surface start at 79,000 ft above Zero. As 
they fall to 65,000 ft they travel westward approximately 75 miles. 

In falling the next 5,000 ft they travel approximately 10 miles south- 
vest. Thereafter they swing around through the south, travel east¬ 
ward between 50,000 and U0,000 ft, gradually turn to the northeast, 
and at about 15,000 ft settle into the surface trades for the last 20 
miles of westward travel to .point n A", Since the hodograph is drawn 
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Figure 1 

Sample Holograph 








for the convenient rate of 5,000 ft per hour the particles arrive at 
"A" 15.8 hours after zero* Particles arriving at "P" from 79,000 ft 
must spend a slightly longer time in each altitude layer; their tra¬ 
jectory is similar to that for the fas ter-falling particles but is ex¬ 
panded. The relative time for particles to arrive at B P B compared to 
n A" is given by the ratio of the vector sums or by the factor 60/U5 ■ 
1*33* The time in hours is therefore 21,1 hours. 

Nov we have assumed a relationship between fall rate and particle 
diameter given by Stoke 1 s lav. Specifically, for land shots at the 
Pacific Proving Grounds, we assume that 75-micron particles settle at 
h,000 ft per hour. Particles arriving at "P" from 79,000 ft settle 
at 3,7liO ft per hour and are therefore slightly smaller than ?5 
microns, Stoke's law gives 72ek microns. Using equation (3) with the 
numerical values 72.1* microns and 21.1 hours one obtains a dose index 
at "P" from the intercept "A" of approximately 12, 

A similar computation for the dose index from intercept ”B a gives 
the much higher value 203. The total index at "P" is therefore 220. 
This means that for 10 megatons one would expect an integrated life¬ 
time dose of 220 roentgens in the vicinity of point n P". 


It is obviously time-consuming to compute the index for a large 
number of points in the field in the way just described. In order to 
speed up the computation of dose indices the nomogram in Fig. 2 was 
constructed from an alternate form of the dose index equation. The 
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where ” is the hodograph fall rate (5,000 ft per hour), "h 0 w is the 
intercept height in feet, K X S W is the Stoke's law constant (U,000/5625 
in our units), n R (> n is the intercept radius in miles, and ”R n is the 
distance in miles to a point on the bearing n Q ”, One should note 
that for constant 51 Q " the intercept height and distance are also 
constant and that the dose index along this bearing decreases as a R“^ w , 
Consequently if one can determine the dose index at one point along the 
bearing the indices for other points can be read very quickly from a 
straight line of slope -3, on log-leg graph paper. The obvious point 
to determine an index is at "R ”, For all bearings the dose index at 
R " R q is determined by the intercept height alone. It becomes very 
simple therefor© to make a suitable nomogram by drawing a line of 
slope *»2 which gives the dose index n D * for R 3 R 0 * The doses at 
any "R" are then obtained by placing a straightedge with slope -3 
through the point n D Q , R 0 K « For quick raadix^ several arbitrarily 
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spaced lines of slope -3 have been dr aim in Fig. 2, 

Using again the sample wind structure of Fig. 1 one obtains the 
index for point "P® from the intercept "B" by entering the nomogram 
at h “ 39,000 ft at the top of the page, coming down to the unweighted 
height line at index 110, moving left to a point above the distance 
R 0 ■ 73 miles, then up a line of slope -3 to index 210 above R “ 60 
miles* 

The nomogram in Fig. 2 has two height lines, the curved one 
corresponding to a weighted height distribution. This second curve 
is believed to be somewhat more realistic inasmuch as it reduces the 
weight of the very low and very high parts of the cloud and increases 
the weight of the middle cloud. This rather general weighting is based 
entirely on observation and experience and is not derived analytically. 

With a very little experience in the use of the nomogram one can 
from the hodogrsph prepare the plot shown in Fig. 3 in approximately 
one-half hour. The weighted height line has been used in this plot. 


For operational use a 10 R dose is a handy dividing point between 
“dangerous" and "harmless" levels, A dose index of the order of 10 is 
therefore perhaps the most important of the family of lines* An index 
of 10 is however established at rather late times — it is completed 
for the 90,000 ft level at approximately 23 hours, and its extreme 
point may be a considerable distance from the origin. Consequently 
it is in general desirable to make appropriate corrections for time 
and space changes in the wind structure. 

The method devised ie based on the same general assumptions as for 
the static case and is again simple enough to permit the preparation 
of a plot in a short time. The complication here was that a complete 
weather analysis of the entire area could not be made every few hours 
by the limited crew at the Eniwetok Weather Central. Consequently it 
was necessary to have a forecaster present to make appropriate 
corrections as the plot developed. 

Since in the dynamic case the trajectories of particles of differ® 
ent sizes falling from th© sme height are no longer similar, the 
terminal points do not fall along a straight line running out from the 
origin. The index falls off inversely as the cub© of the distance 
along the line of terminal points instead of along a line of constant 
bearing. 
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It Its possible but at present very time-consuming to compute the 
lines of terminal points for a set of altitudes complete enough for an 
adequate plot. The best that could be done in the time available was 
to locate the 10 R point for about a dozen altitudes and base the 10 R 
line on these points alone, A complication arises when the points 
fall at different distances on roughly the same angular bearing. This 
result corresponds to a contribution from two or more altitudes and 
means that the 10 R line must be drawn some distance beyond the outer¬ 
most calculated point* The adjustment is easy to make approximately by 
simple application of the inverse R^ relation to the inner points to 
find the contribution at the outermost and a second application of in¬ 
verse to the sum. To take into account the height weighting one 
locates the terminal point and then contracts cs* elands radially, 
again by inverse R^, These approximate methods are an unhappy feature 
of the system but are necessary if the plot is to be completed in a 
reasonable time. We were unable to discover any quick, reliable 
analytic or graphical method of making these adjustments more precise. 


The mechanics of locating the 10 R points (before any adjustments) 
are straightforward enough. The dose index equation (3) gives the dose 
in R for 10 megatons. For a predicted yield different from 10 megatons 
one may write a dose equation 



(?) 


where"I B is the dose in roentgens, "W is the yields in megatons, and 
"d” and n t n are defined as before. This equation follows, of course, 
from assumption (U)« The farm (?) is not convenient for calculations, 
since in the construction of a trajectory for a falling particle one is 
concerned only with the time the particle spends in a stratum character¬ 
ised by a fixed mean wind* It is convenient then to introduce Stoke 5 s 
law and to solve equation (?) for the fall rate * f v w . The result is 
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(Sa) 

where "K B is the constant 1*000/562?, n h H is the starting height in 
feet and®the unit of "v" is feet per hour. The time spent in a given 
stratum n Ah B is them 



At ■ Ah o (6) 

v 

Equations (?a) and (6) were used to prepare a worksheet listing the 
strata, times in the strata, and accumulated times from the start 
to the strata. Forecast mean winds were then entered on the work 
sheet as the trajectory plots were made and the times and space 



locations became known. In the actual cases computed the changes in 
wind structure were never very significant and the resulting plots 
differed little fWmn those made from the static hodograph. The 
dynamic system does, however, have the capability of handling a rapidly 
changing situation should it appear. 

Conclusions 

The fall-out forecasting systems described above have a large 
number of very obvious deficiencies which one would hope to remove by 
the next operation. For example one would hope to include a good 
representation of the particle size distribution for both land end 
water shots. The question whether area or volume of a particle Is 
more significant for the deposition of activity should be investigated. 
A careful estimate of the height distribution of the activity should 
be made. Some attention should also be given to the effects of finite 
lateral cloud dimensions and to the spreading of the cloud. These are 
points which will refine the system. 

Dr. Tom Whit© of LASL began to make such refinements during Castle 
in an attanpt to remove one outstanding defect of this system, its 
inherent inability to reveal ary detail near the origin. His method 
is described in a separate report. 

Two difficulties with the system cannot be removed at all. The 
case where the hodograph passes directly over the origin is not un¬ 
likely and gives embarrassing results, since R ■ 0 and the dose in¬ 
dex vanishes except at the origin. This result is at least circularly 
symmetric. The other awkward case is that in which one car more ele¬ 
mental wind vectors is radial, since one then has a continuous 
altitude range contributing to the dose index along a bearing in place 
of one or two discrete altitudes. This is in reality precisely the 
kind of wind structure which leads to very intense narrow bands, but 
the system as set up cannot estimate the magnitude of the dose index 
along such a bearing. One could perhaps handle this situation by 
revising the method such that each 5,000-ft altitude interval made a 
contribution as determined by the mean height of the interval along a 
bearing line through this mean height. Such a change would introduce 
difficulties in evaluating the contributions at one point from several 
heights since the bearing lines through the mean heights would not in 
general coincide. 

One difficulty will continue to plague us even 
refined system® At this time it is very difficult to obtain a weather 
forecast precise enough to justify the effort which must be put into a 
good fall-out forecast. Because of the inherent uncertainties in 


archives 

with any sort of 


315 





weather forecasting one is tempted to conclude that a refined system 
of fall-out forecasting will for the time being be most useful for 
post-shot analysis and that a detailed fall-out forecast based oh a 
normal weather forecast would be misleading. The weather forecasting 
presently available is of the highest quality but the cumulative errors 
resulting from small variations about the forecast mean winds can 
result in a pronounced change in the fall-out pattern. It appears de¬ 
sirable at this time to retain a crude system which is relatively in¬ 
sensitive to small variations in weather structure and which presents 
conservative upper limits to the fall-out hazard, 
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PART II 

A METHOD OF ESTIMATING RADIOACTIVE FALL-OUT 
T. N. Whit© 

Los Alamos Scientific Laboratory 


The method described herein is designed with an objective that is 
intermediate between operational requirements, and the requirements 
of a strictly scientific investigation. It is designed to include 
what are assumed to be the most important factors that determine a 
fall-out pattern, with the idea that we might find out enough about 
what is going on to produce a good simplified mettod for operational 
use. One simplified version that was used for local fall-out fore-* 
casting is described in Close-in Forecasting by New Techniques 
Developed after BRAVO, Tab D '’Fall-out Forecasting Techniques" of the 
Task Force Castle Report, It seemed good enough to justify further 
investigation of the basic ideas as applicable to any range of dis¬ 
tances over which a constant wind field could be assumed, 

‘ The basic assumptions of the mathod are as follows: 

The whole cloud, up to its height of stabilization, is formed 
instantaneously at the time of detonation. This is what we call the 
"initial cloud". 

In any height layer of the initial cloud, the concentration 
(radioactivity per unit volume) is distributed according to the 
Gaussian Law 

c(h,r # a Q ) » c Q (h) exp (-r 2 /a 0 2 ) E ARCHlVESf 

where c(h) is the central concentration at height h, r is the radial 
horizontal distance, and & Q is a "spread parameter" (analogous to 
standard deviation) that is also considered to be a function of height. 
From this assumption it follows that the total amount of radioactivity 
in a slice of unit vertical thickness is v e 0 (h) , 

Throughout all of any such layer, the radioactivity is dis¬ 
tributed normally with respect to the logarithm of the rate of fall 
of the particles. Thus at any distance r, the fraction of radio¬ 


activity that falls with speeds in the range f to f + df is given by 




where f(h) is the fall-rate for particles of greatest radioactivity, and 
cr (also considered to be a function of height) is the standard 
deviation of the logarithm of fall-rates, weighted according to radio¬ 
activity? f(h) and a- (h) are constant throughout the layer. 

The rate of fall of any particle remains constant until it 
reaches the ground. 

Any particle that starts from the central axis will follow a 
path strictly in accordance with the wind pattern, while all other 
particles that fall at the same rate from the same level will diffuse 
laterally from the central particle in such a way that the gaussian 
distribution is maintained, 

\ , _ 

During this process, the increase in the spread parameter is 
described by 



where S is the distance travelled by the central particle until it 
reaches the ground* (It is to be noted that S is not the straight- 
line distance from the origin to the landing point unless all winds 
at all levels are in the same direction. $ and m are parametric 
quantities that may be used to describe the amount of diffusion. They 
are not at present regarded as functions of height. (The quantity p 
is merely an abbreviation for the quantity in brackets). 


From these assumptions it follows that the dose rate on the ground 

is 

1 'S—f / Co(h) cxp 2 

jfWo 07® P ~r-r-7~ -iUr- dh df 

<r(h) */2 
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where K is dose rate per unit of surface concentration, H is the height 
of the top of the cloud, and r is now the distance from the point at 
which the dose rate is estimated to each of the landing points of cen¬ 
tral particles. These landing points will depend on the wind pattern 
below the level from which the central particle originated, so that r 
is a function of h. The landing points also depend on the rate of 
fall, so that r is also a function of f. Changing from rate of fall 

to time of fell, one obtain, POE ARCHIVES 

r 2 - [ t 5 (h) -X ] 2 + [t f (h) -T | 2 
where (X,Y) are the rectangular coordinates of the point where dose 




pate Is estimated, and ^ are the wind components in the same co¬ 
ordinate system, averaged up to the height h. 

We may also express 

p - 1 + -fr frjn- 

./3*o 

noting that w is the average speed regardless of direction, (u and v 
are not , in general, the components of {J ). x his expression is correct 
if one is satisfied that the diffusion depends on the total horizon tal 
distance travelled by a central particle. If one wishes to assume "that 
the vertical distance should be included, p becomes much more compli¬ 
cated 8 

The significance of /3 and m can now be visualised® If m ■ 2, then 
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so that the lateral, dimensions of any segment of the cloud will in¬ 
crease as if the segment had come from a point source located at a 
distance g a o upwind® The size increases linearly with distance 
travelled by the central particle. If m is greater than 2, the borders 
of the cloud will diverge mere rapidly, and if m is less than 2 they 
will diverge less rapidly® One can prevent any increase in size either 
by making @ infinite or by making m equal to zero. This method of 
describing the diffusive process is similar to that of Sutton, but not 
exactly the same. 

If m ** 2, then at sufficiently large values of t, tho area covered 
by a segment of cloud is proportional to the square of the time, as in 
Felt’s method. However, the proportionality factor varies, as £ 
varies with height, and further, the overall average proportionality 
factor changes with the overall strength of the wind field, and in 
these respects it differs from Felt’s method. 

Returning to the basic equation, the change of variable from f to 
t changes the argument of ths logarithm to 
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and d£ becomes dt 

T HE 


One notes also that concentration© in the initial cloud must be 
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reduced to those that would have existed at that time for which the 
dose rate is being estimated. 

The information that is needed for a calculation is then: 

The winds pattern at heights up to H, 

Hj the height of the top of the cloud, 
a , the initial spread parameter, or radiological radius, 
as a funct£on of height, 

C 0 , the central concentration as a function of height in the 
initial cloud, adjusted to the time of dose-rate estimation, 

f, the logarithmic mean rat® fall (weighted according to 
radioactivity), as a function of height in the initial cloud, 

v, the logarithmic standard deviation of this distribution 
as a function of height in the initial cloud* 

0 t diffusion parameter, described above, 
m, diffusion parameter, described above* 

Testing of the method requires the use of high speed computing 
machinery. With such machinery one can make many changes in the 
quantities described above, proceeding on a trial and error basis* 

In order to achieve some degree of objectivity, the following approach 
is adopted. 

The logarithm of the ratio of calculated to observed dose rate is 
estimated at a number of points for a given shot. This quantity is 
called y (gamma). Then the mean gamma, and the statistical variance 
of the individual gamma about the mean gamma, are calculated. This 
process is repeated for a number of values of some parametric quantity, 
beta for example. One then plots the variance against beta, and 
selects as the best value the one that gives the least variance. One 
then moves on to other parametric quantities, and treats them in the 
same way, hoping that there is not too much correlation between the 
effects of the different types of parameters. It will be noted that 
this application of the "least squares" method discounts the overall 
ratio of calculation to observation. In principle it is possible to 
get zero variance (an exact fit)when each calculated value is, for 
example, exactly ten times the observed value. One would then sus¬ 
pect that 90$ of the radioactivity had remained in the crater. If, 
however, one should obtain a good fit but with only 10% of the observed 
activity accounted for, one would have to consider other possibilities® 
One would first look to see whether any large fraction of the activity 
was excluded from the calculation. In numerical integration it is not 
practical to go all the way from zero to infinite time, and part of the 
activity might have fallen outside of the time range chosen. If this 
explanation fails, and if the fit is really good, one has to conclude 
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that the least squares criterion, as applied here, is not useful. 

We have not yet encountered this particular obstacle. 

The method of approach is subject to a hackneyed old criticism 
running as follows: by subdividing the height of the cloud at will, 
you can obtain as many discrete parametric values of a , c 0 , f, etc, 
as you wish, so that you should be able to fit any nurrSer of obser¬ 
vations exactly. This is true, in principle* But if on© can get a 
set of values that look reasonable, and can be scaled in a reason¬ 
able way with yield over a wide range, the method can serve a useful 
operational purpose even though the values might be scientifically 
unsound. 

For machine integration, using the IBM Model 701 "Defense 
Calculator", those parametric quantities (a, c, o- and f) and the mean 
wind components ( v, w ), which are functions of height, may be 
loaded as tables of data. The total height of the initial cloud, H, 
is subdivided into M equal layers, each identified by an integer i s 
0,1, 2, — (m-1). The time variable is laid on a logarithmic scale, 
each time being identified by an integer j ■ 1, 2, 3, — N, Storage 
limits the maximum value of M to 32, N may be any value that doesn’t 
take too much machine time, and the minimum and maximum limits of the 
time integration may be changed at will. The exponential factor in 
the formula is recorded as ecjto if the absolute value of the exponent 
exceeds a value A which may be as large as 10, 

The coding is so arranged that the time integration is performed 
first and the height integration second. At each location, the frac¬ 
tion of the dose rate that comes from each initial cloud layer is com¬ 
puted and may be printed along with the calculated and observed dose 
rates and the co-ordinates of the location. Or one may by-pass 
this printing and obtain only the statistics: mean gamma and variance 
for a preselected series of locations. 

The codes are not yet frozen, and additional features are being 
added from time to time. We have two codes (1) a fast "fixed-point" 
code as outlined above, and (2) a slower "floating-point" code that 
ie more precise and which is more flexible in some respects. 

Since starting on this problem about seven months ago, 
able part of the tim© has been spent on coding and de-bugging, which we 
undertook ourselves in order to learn how to use the Model 701, Using 
Bravo fall-out data, we demonstrated that least squares solutions could 
be obtained for the various parametric quantities involved. However, 
th© "best" values, as selected in this way, gave distant fall-out 
predictions that were only 20 to 30$ of the observed values, and 'the 
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"fit" was not good* We then turned attention to Nevada data for a 
while, and became interested in an approximation that seemed to offer 
a hope of eliminating on® of the two steps in the double integration* 
Before this possibility had been fully explored, Mr. Vay Shelton, 
Livermore Operations Division, joined forces with us, and we worked 
together for a week on UK-1 and UK-7* Mr, Shelton then took our codes 
to Livermore and continued working on the Nevada data, while we turned 
attention again to the Bravo shot* Mr. Shelton has reported recently 
that the method gives satisfactory results for UK-1 and UK-7, and he is 
continuing work on other shots. We have concentrated on the problem 
of predicting the Bravo fall-out on the assumption that practically 
all of the activity in the initial cloud was located above the tropo- 
pauso. To date, our method of calculation has not bsen able to give 
satisfactory results, even when the winds were arbitrarily twisted 
to make the fall-out occur In more nearly the right place. At this 
point w© feel, therefore, that we do not have any cloud modal in which 
we have confidence. We have merely a mechanism of calculation, the 
value of which has not yet been proven as far as Bravo is concerned* 


Lacking any satisfactory cloud model for Bravo, w© tackled ths 
’’homework” predictions on a guess-work basis, which does not Justify 
a description. Tho values that we used in the calculations are: 
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The values of c were taken as constant up to the tropopause, and 
thereafter decreased with air density. (The machine program requires 
only the entry of relative values, from which the actual values ar® 
adjusted so that the total radioactive content of tha cloud Is in 
accordance with toe yield). Logarithmic mean rates of fall were 
assigned to the 16 layers of the cloud as follows, counting frcs tha 
bottom, (Rates are in knots) 

Layers 1 thru h 20 

Layers 5 thru 8 7,b 

Layers 9 thru 12 2,7 

Layers 13 thru l6 1 


Sigma: 




S ?! ‘ 
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(Editor's note: The four contour figures which follow were received 
too late for inclusion in the comparison done by Cdr. Paine. 
The contours shown in that comparison for LASL were submitted 
prior to the symposiumj and, being scaled by simple yield de¬ 
pendence, did not include allowance for different cloud dimen¬ 
sions as those presented here do.) 
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PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS FUNCTIONS 
Of YIELD AND METEOROLOGICAL CONDITIONS, RAND METHOD 

PART I 

Dr. R. R. Rapp 
Rand Corporation 


The Formation of the Atomic Cloud 

Dr. Kellogg has already presented the physical picture of the atomic 
cloud upon which the Rand model is posed. There are, however, a few points 
to stress. First, as shown by the first sketch in Figure 1, there are 
large amounts of earth taken into the fireball in the very early stages 
of formation. Next, as is indicated in the other sketches of Figure 1, a 
large amount of earth is carried aloft with the fission products by the 
extremely stable toroidal circulation. This toroidal,circulation is 
caused by the intense thermal gradients of the fireball. Thus, large 
quantities of vaporized, liquified, and pulverized earth are in intimate 
contact with the newly formed fission products from a few milliseconds to 
a few minutes after detonation. By way of contrast, it is believed that 
the earth, which is seen to enter the fireball of an air burst, flows up 
through the "hole in the doughnut," circles around, and falls to the earth 
without actually encountering the radioactive material. For the surface 
hurst, where there are' large particles, after the cloud stops rising, 
there will be radioactivity present on particle sizes ranging from thou¬ 
sands of microns to about a very few microns spread through a region of 
the atmosphere. For a large yield burst, this region may look something 
like Figure 2. At this point, it is not necessary to quibble about the 
details of the configuration or the activity distribution with space and 
particle size, but simply to say that for a given bomb, there exists a 
.space and particle size distribution of the activity produced. It is 
interesting to note, however, that on the basis of currently available 
data, there is marked disagreement on the distribution of activity in 
space. In the Rand model, as will be pointed out later, the size of the 
active cloud is smaller than the photographic observations of the CASTLE- 
Bravo cloud would indicate. Thus, although the figure indicates an active 
particle falling from the edge of the cloud, this is schematic and has no 
significance in the Rand model. 

Doe Archived 

The thermal and kinetic energy which produced the cloud and caused “ 
it to rise through the atmosphere has, by the time of stabilization, been 
dissipated and the particles will begin to fall under the force of gravity. 

If particles in the size range from 100 to 1000 are considered, it be¬ 
comes necessary to consider the aerodynamics of the particle. No matter 
what value is used for a drag coefficient, the density of the air must 
enter the equations and the particle will fall more rapidly In the less 
dense atmosphere at high levels. 
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Figure 1 

Surface Burst - Schematic Representation of the Manner in which 
Earth is mixed with Fission Products and carried to great Altitudes 




the Idealized Path of One Particle 




Figure 3 shows some possible assumptions about the fall velocities 
of the radioactive particles- The Rand curves are based on the aero¬ 
dynamic fall of spherical particles as set forth in an article by Langmuir, 
whereas the NRDL curves are based on the formulae of Dalla Valla. The 
line labeled Stokes' Law is included mainly to show that it does not give 
the proper velocity variation with height and that it over-estimates the 
velocities of all except the very smallest particles. The variation of 
velocity with height is extremely important when computing the distance 
a particle will travel, since the more rapid fall through the upper layers 
means that the particle will spend a relatively shorter time in the upper 
wind flow than it will in the dense air near the surface. Referring back 
to Figure 2, it may be noted that the vertical section shows the particle 
falling faster in the high levels them it does in the low levels. The 
horizontal section of Figure 2 shows the erratic path a particle may take 
due to the action of the wind. It should be noted that the term wind 
implies all the components, scales and variations of the air movement. 

It may be informative to look at some of the wind factors which must be 
considered. 

The picture of the motions of the atmosphere that one gets is very 
much a function of the sensitivity and spacing of the measuring instru¬ 
ments. The established weather networks, in which the instruments are 
spaced at large intervals, indicate a rather smooth flow of air with 
perturbations from the mean motion which are the size of cyclones and 
anticyclones. These perturbations change with time so that time and 
space variations must always be considered. For example, a steady change 
of a velocity component of 2 knots per hour could make a position error 
of 24 miles in 12 hours. If more sensitive instruments are used at a 
clbser spacing, it is found that there are smaller perturbations imposed 
on the larger ones. There have been relatively few experiments with close 
networks, but the experiments which have been made show that perturbations 
may be found on almost any scale of measurements. The total motion of the 
atmosphere can, therefore, be represented by a spectrum of motions on 
different scales. Such a representation provides a useful mechanism for 
separating the motions into those which can be readily measured and those 
which must, for lack of observation, be treated statistically. 

DOE ARCHIVE^ 

An important feature of the spectrum of motions in the atmosphere is 
its variability in time and space. Large gaps in the spectrum apparently 
appear and disappear as the characteristics of a mass of air change with 
time because of mechanical and thermal effects. Sketch A, Figure 4, shows 
the effect of the large scale motions which move the cloud bodily. Sketch 
B show's the effect of motions on a scale comparable to the cloud dimensions. 
These motions tend to distort the cloud from its original shape and may 
cause an error up to 20 percent in the computed position of a particle 
on the ground. Sketch C indicates that motions on a scale, small compared 
with the horizontal cloud dimensions, tend to spread the cloud" by drawing 
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Figure 3 

Comparison of the Fall Velocities under Several Assumptions 

as a Function of Height 








out small wisps of.cloud matter from the edges and carrying clear air into 
the cloud; this type of phenomena is called eddy diffusion. Eddy diffu¬ 
sion spreads the cloud very slowly and may be neglected for short periods. 
Sketch D shows how the change of wind with elevation tends to shear the 
cloud into a long ribbon. The last sketch, E, is a schematic view of the 
integrated effects of these various motions on the cloud. 

The breakdown of the atmospheric motion into different scales is, of 
course, highly artificial. Such a breakdown is dictated by the available 
measurements and by our ability to analyze them. The gross movement and 
the gross deformation in the vertical can be conveniently, if not accu¬ 
rately, handled by recourse to actual weather observations. The horizontal 
deformation and the eddy diffusion are, however, on much too small a scale 
to be observed by our present network of stations. With such a small 
scale of motion, it is necessary to resort to some type of eddy diffusion 
theory or to rely solely upon empirical evidence gathered from the tests 
that have been made. Since neither the empirical evidence nor the eddy 
diffusion theories show much spread, diffusion has been neglected. 

A bit earlier, the distribution of activity with particle size was 
mentioned. There are two methods for arriving at such a distribution. 

One is to measure a large sample of collected particles, as Is done by 
KRDL, and the other is to choose, by trial and error, a distribution 
which fits an observed case. At Rand, the latter method was applied to 
JANGLE S to derive a distribution. AKDC uses basically the same tech¬ 
nique, but with different assumptions. Figure 5 shows a conparison of 
Rand, NRDL and ARDC distributions. The studies made by the Weather Bureau 
Group on reconstructing activity distribution with size and height for 
tower shots had found still another distribution. It is quite possible 
that there are changes in the distribution from shot to shot. Thus, in 
order to compute fall-out, it is necessary to know the space distribu¬ 
tion ofttebris, which is known only by Inference; the distribution of 
activity with particle size, for which no general law is known; the fall 
velocity of the active particles, which has not as yet been determined 
experimentally; and the wind in all its particulars, a requirement which 
is economically unfeasible. It is obvious that with the uncertainties 
inherent in the data now available, it would be unreasonable to expect 
any model to reproduce the observed fall-out with a great deal of accu¬ 
racy. It might be possible to compute within five to ten percent of the 
relatively close observation where the heavy particles spend the least 
time in the wind field. For more distant observations, it may be possible 
to be within a factor of two for the highly significant observations, 
but on the fringes of an area, fall-out may be computed where none was 
observed and vice versa due to the unknown factors of the wind variation. 
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Hie Model 


With all this uncertainty, it is still necessary to fix on some of 
the physical parameters in order to make computations. Therefore, a 
model appropriate to CASTLE-Bravo will be presented (see Figure 6). The 
assumptions used at Band are, first, the distribution of activity with 
particle size deduced from the JANGLE surface shot data. We are currently 
of the opinion that the high temperature and large volumes of soil pro¬ 
duce the characteristics of the activity distribution, but at present we 
have little to offer in support of our thesis. Our computations were 
based on early measurements of cloud diameters and heights, and have not 
been modified because the smaller sizes seem to fit the observation 
adequately. Our early model assumed homogeneous distribution of the 
material in the cloud and stem. We also assumed that the cloud and stem 
were cylinders. Recently, we have modified our concepts and adopted an 
exponential decrease in the activity within the cloud. 

The assumption on the fall velocities is that the particles fall as 
though they were spheres. 

For the winds, we have used two assumptions, both quite simple. 

First, we assumed that the winds at shot time apply to the entire descent 
of the particles. As the next more complex assumption, we have assumed 
that the time variation of the wind was more important than the space 
variation. These assumptions are shown for the Bravo shot in Figure 7* 

We have neglected the effects of medium scale motions and turbulence 
simply because there is no simple way of including these effects. They 
would tend to spread the cloud, but they would not spread it in a uniform 
way. Thus, if a spread parameter were included and an even spread of the 
material were assumed, it would tend to cover a large area, but would 
smooth out the value of the hot spots. Without a spread parameter, the 
concentration of material is about right, but the clumps are not distrib¬ 
uted over the proper area. Since we do not believe that we can place the 
debris accurately, we have decided to accept the position error and hope 
to keep the proper concentrations. 


The Synthesis 
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Having outlined the physical parameters of the Rand model, it is 
now necessary to go on to a method of computation. The machine method 
will be discussed later. This paper will be confined to the hand compu¬ 
tational scheme. 

Consider a section of the cylindrical cloud extending from h to h 
/ dh containing particles of a size range r to r / dr, Let b .± (r,h) be 
the activity contained on these particles in this height interval over 




Figure 6 

A schematic of the RAND Model showing the Spatial Dimensions of the Cloud, 
the Relative Distribution of Activity with Height, and the Distribution of 

Activity with Particle Size 


MODEL for CASTLE 'BRAVO 
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h for the CASTLE-Bravo Event. The Main Plot is Following the H Hour Winds, 
The added Sector Shows Time Variant Winds according to the Oahu Research 

Center Report. 
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a unit horizontal area of the cylindrical section. When all of these 
particles fall to the ground, the surface concentration will he simply 

er(r,h) ■ ai (r,h) 

Now it is necessary to determine ai(r,h) in the light of the assumptions 
of the model. First, it is assumed that, within the mushroom, the dis¬ 
tributions in the height and particle size are independent. That is, 

a-i( r >h) a K A (r) G (h) dr dh 

Where K is the total activity per unit horizontal area of the cloud and 
A and G are distribution functions. Next, it is assumed that A(r) is 
the JANGLE distribution. This A(r) is now a density function and is the 
derivative of the curve, presented as the Rand curve in Figure 5* For 
G(h), it has been assumed that the distribution of particles follows 
the density distribution of air so that G(h) »» G Q e-&(H - Hfl) where H is 
the height of a point in the mushroom and Rg is the height of the base 
of the mushroom. 

Thus for the small cylindrical section 

«r(r,h) ° a^h) - KA(r) G c e"® (H " H B> * 

This is the surface density of activity for a point on the ground due to 
the particles of size r to r / dr from the height h to h / dh. 

It is possible to find the coordinates of this point on the ground 
as follows: If V x (h) and V y (h) represent velocity components, then 

dx = V x (h) dt and dy a V y (h) dt 

dh dh 

But W(r,h) - ^ thus dt ” v ^ere W(r,h) is the vertical 

velocity. These are derived for all (r,h) from the aerodynamic fall of 
spheres, a sample of the variation was shown in Figure 3. Now 

0 , V 

D _ /• V y(h) 

y " / vfoT DOE ARCHIVES, 

h 

D x and D y are distances traveled by a particle of radius r, from a height 
h. The actual position of a particle will depend on its initial position 
in the cloud. If the origin is placed at GS and (xj_, y-jj is the initial 
displacement from GZ, 


Dx '/ 


Vjh) 


dh 


X - Xi / D x 
K 

a x ± / Z 

j = 1 

For any given point, say (xq, y 0 ), there will be many sets of r,h values 
which can reach the spot. If a grid is made of the x,y positions from 
the center of the cloud for a large number of r,h values, any r,h value 
from the center of a cloud which lies within a cloud radius will contribute 
to the activity at the point (xq, y 0 ), since the spread of the cloud has 
been neglected. 

Since tr (r,h) is the surface density of activity due to one particle 

Le total density at XQy^ will be 

h2 

/ A(r) dr G(h) dh DOE ARCHIVES’ 

h l 

A(r)dr 


G(h) dh 


h 2 

/ d H ± 

h l 

a constant, depending on the size of the cloud 
and the total activity in the mushroom, times the integral which is an 
area on an Ri, Hi plot. 


In order to compute the activity at a point, assuming all the activity 
has fallen out, the (x,y) positions of a set of (r,h) pairs starting over 
GZ are mapped onto the ground. A circle the size of the mushroom is 
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drawn about the point. The (x,y) positions of the edge of the circle 
are converted to (r,h) values. These (r,h) values are then converted 
to (Ri>Hi) values by means of appropriate graphs. These are then plotted 
on any convenient piece of rectangular graph paper and the enclosed area 
is measured. This area represents a fraction of a unit area column 
through the cloud which arrives at the point x^y 0 , multiplying this 
fraction by K gives the surface density of activity at this point. 

CASTLE "Bravo 

Referring again to Figure 7> the time variant winds were used to 
make the (r,h) plots shown across Rongelap. The 12-hour dose rate at 
the northern tip of Rongelap from the time variant winds computed by 
the method just outlined was 52 R/hr. The time of arrival of the debris 
was from 5 — 1^- hours. By using the same system on the Bikini H hour 
winds, a 12-hour dose rate of 45 R/hr was computed at a point with 
similar (h,r) values. The observed 12-hour dosage at this island was 
reported to be 96 R/hr. Thus, with the best winds available, our com¬ 
putations are well within a factor of two of the observations. Using 
the Bikini winds, we are still nearly within a factor of two, but are 
about 55 miles off in position. The calculations, based on H hour 
winds which should be valid for this close-in fall-out, close to G-Z 
appear to be very accurate. 
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PART II 

UTILIZING THE IBM 701 ELECTRONIC DIGITAL COMPUTER 

S. M. Greenfield 
Rand Corporation 

The previous paper by Dr. R. R. Rapp has indicated that it is 
possible to construct a model which will adequately represent the 
physical phenomenon of fall-out. In working with the model he pre¬ 
sented, we at Rand have found that it provides us with as much of the 
fine structure of the fall-out pattern as one might reasonably expect 
when considering the accuracy of the various parameters that go into 
it. By these parameters, I mean such things as the accuracy of the 
wind observations, the dimensions of the cloud at stabilization, etc. 

As one might imagine, however, having once established what can be 
considered a representative model, it is not a sinple task to use it 
for computing a fall-out pattern containing the fine structure that can 
be expected to result from the quality of the input data under normal 
conditions. One does not just sit down and, in a short half-hour or 
hour, construct a complete pattern. Rather, it may take as many as two 
people as much as a day to finish one of these. While a method requir¬ 
ing this many man-hours is certainly adequate if one were only interested 
in one or possibly two or three patterns, it falls far short of desir¬ 
ability when one enters into a general fallout study. This is true 
whether the study involves an investigation of the phenomenon of fall¬ 
out or an investigation of the implications of the phenomenon in an 
operational sense. 

We at Rand are faced with both problems in our fall-out study. We 
would very much like to examine how the patterns change as one varies 
the wind structure with time and space, the distribution of activity 
with particle size, and the concentration of activity with altitude. 

Also, due to the peculiar nature of Rand's work, we have a very definite 
need for fall-out information as an input into an overall weapons employ¬ 
ment study. As such, we have been painfully aware of the shortcomings 
of the hand computation of the fall-out pattern for many months. However, 
at Rand we are very fortunate in having available both an excellent 
digital electronic conputer like the IBM 701 and a very capable staff 
of problem programmers for this machine. Working with two of these 
people, Miss Shirley Marks and Mr. Dave Langfield of our Numerical 
Analysis Division, it was recently found possible to place the model 
that has been described previously, or one with % slight modification, 
on the IBM 701 for computational purposes. As . the case of the hand 
computation method, the machine model recognizes the fact that the r,h 






plot (that is, the positioning of particles, of different radii starting 
from different altitudes, on the x,y grid representative of the ground) 
is the basic bit of information necessary for the pattern calculation, 
and it is toward this end that the initial (and major) segment of the 
program calculation is slanted. However, it should be noted that, due 
to the computational speed of the 701, ve have at our command the capa¬ 
bility of making a much finer grid of (r,h) than was feasible by means 
of the tedious hand computations. Because of this advantage, we have 
chosen to subdivide the atmosphere into 1000 ft. intervals, thiB being 
a limit imposed on us as being about the smallest intervals at which 
winds are reported. In the case of radii of the particles, one-hundred 
intervals of particle size were chosen as being adequate to treat the 
problem. These intervals were chosen so as to contain equal fractions 
of the total activity in the cloud. It was initially assumed that within 
either the stem or the mushroom, the activity distribution with particle 
size was invariant with altitude. We then have, essentially, two matrices 
of particle size and altitude, one representing the stem and one the 
mushroom. Choosing the value at the center of both the particle size 
interval and the altitude interval as representative of both the intervals 
themselves, we can then tag each r,h combination with its appropriate 
fall velocity as described in the previous paper, assuming that this 
fall velocity taken for each r,h applies over the entire interval. Due 
to the fineness of the interval divisions, it is found that this assump¬ 
tion is true to within a maximum of one or two percent. These matrices 
provide us with the first part of our basic input for the machine. It 
should be noted that by choosing the particle size intervals so as to 
contain equal fractions of the activity, we have made these matrices 
independent of bomb yield. In other words, if we can fix on a distribu¬ 
tion of activity with particle size, then each matrix can be constructed 
so as to include all necessary altitudes and then be included in the 
program as a semi-permanent part of the memory. For a specific case 
then, providing the machine with information as to the vertical dimen¬ 
sions of the cloud, and the variation of wind velocity and direction with 


altitude, will allow it to compute the position of the r,h grid on the 
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For this confutation, it is assumed that all the particles are 
distributed along a line extending vertically upward from ground zero 
(the particles present at each altitude are fixed by the distribution 
assumed for the stem and mushroom). Under this assumption, the machine 
then proceeds to solve the following analytical equations: 
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W(r,h) 
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Where V x (h) and V y (h) are the x and y components of the wind as a 
function of altitude, A H is the 1000 ft. altitude interval and W(r,h) 
is the fall velocity applicable tc a given particle size and altitude. 
Choosing an r and an altitude at which to start, the machine then holds 
the particle size constant and merely sums over the altitude to determine 
its x or y position on the ground. This process is presented schematic¬ 
ally in Figure 1 where each dot represents the ground position of the r,h 
combination. At the sane time, consideration is given to the expression 
w/r^K T whiel1 is simply the time spent by a particle in each altitude 
interval. By summing this expression over altitude, we therefore have 
not only the ground position of each r,h combination, but also the time 
that it arrives at the ground. 


As was stated previously, each r,h combination represents a given 
fraction of the total activity. Recognizing further that this fraction 
is contained in a slab whose thickness is 1000 feet and whose cross- 
sectional area is the cross-sectional area of the mushroom or stem, then 
the machine, having been given the yield of the bomb and the dimensions 
of the cloud, can calculate the "infinite plane" gamma radiation effect 
of this slab when deposited on the ground. To do this, we utilize the 
well known relationships for gamma radiation from fission products 


Number of 1 hour megacuries/KT - 300 


( 3 ) 


fgrlf? at 1 hour = r i = (3.09) 

A 


( 4 ) 
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Where F(r,h) is the fraction of the total activity associated with * 

a given particle size at a given altitude, and A is a cross-sectional 
area of the mushroom or stem in square nautical miles, and 3*09 is the 
converstion factor for megacuries per square nautical mile to roentgens 
per hour at 1 hour. 


r l t 


■ 1.2 


( 5 ) 


Where t is the time in hours, and r^. is the dose rate at any time t. 
And finally 


At 


■/ 2 


dt 








where is the dosage accumulated over a time interval between 
t^ and tp- 

It was decided that the machine would calculate the dose rate at 
t ■ 6, 12, 2b, and 48 hours after the bomb explodes, and the dosage 
accumulated between 0 - 6 , 6 - 12, 12-24, and 24 - 48 hours. Implicit 
in this is the understanding that for any time or time interval occurring 
earlier than the computed time down for a particle, the dose rate and 
dosage that the machine will report will be equal to zero. 

Having available many thousands of r,h combinations, positions on 
the ground and their radiation contribution as a function of the time, 
we can them proceed to calculate radiation levels at specific points 
within this pattern. Due to the fact that it has been found that about 
100 points scattered representatively throughout the r,h grid are suf¬ 
ficient to determine the fall-out pattern, it is pointless to attempt to 
compute the radiation level at each of these many thousands of points 
and present them as a final output. What was done then was to choose 
125 of these (r,h) combinations in such a way as to make them representa¬ 
tive of the entire pattern. As illustrated schematically in Figure 2, 
this was accomplished in the following way. It is assumed that a repre¬ 
sentative set of grid points would have a density that was directly 
proportional to the gradient of density of radioactive material on the 
ground. Further, it is felt that the variation in time down for particles 
at the same altitude is a fair approximation to their relative positions 
on the ground. Under these assumptions, par-tide sizes were chosen at 
each 10,000 feet of altitude shot for the following times down: 

Every .2 hour for the first hour. 

Every -5 hour for the next four hours. 

Every hour for the next five hours. 

Every 2 hours for the next ten hours. 

The r,h combinations resulting from this type of selection were taken 
as grid points with some degree of confidence that they truly are repre¬ 
sentative of the entire pattern. It should he noted that having once 
chosen our grid points, or r,h combination, ve need not concern ourselves 
with the problem of how the pattern changes as a function of different 
wind inputs. This is due to the fact that grid points chosen in this 
manner will appropriately arrange themselves within the pattern with 
the wind so as to maintain their representativeness. 

It is in the final stage of this machine calculation that the 
dimensions of the cloud really come into play. It should he remembered, 
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SELECTED PATTERN OF GRID POINTS AT 
WHICH RADIATION WILL BE COMPUTED 



Figure 


SUMMATION OF THE RADIATION CONTRIBUTION FROM PARTICLES 







as vas stated previously, that each r,h point in the grid actually 
represents the center of a circle that is the projection of a slab of 
material onto the ground (slab 1000 ft. thick, cross-sectional area 
equal to the cross-sectional area of cloud). Furthermore, it is assumed 
that every ground point within this circle receives an equal radiation 
contribution from this r,h interval. Conversely, it can be seen that if 
a circle of radius equal to the radius of the stem or mushroom is drawn 
around any point on the ground, it will include all r,h combinations 
that will contribute to the radiation found at this point. This is 
illustrated in Figure 3- Figuratively, then, the machine draws a stem 
circle and a mushroom circle around each of the 125 grid points and 
simply sums the contributions from the appropriate r,h combinations. 

The final output from this machine calculation is, therefore, the co¬ 
ordinates of 125 points on an x,y grid (whose origin is ground aero), 
the dose rates found at these points at H / 6, H / 12, H / 24, and 
H / 48 hours, and the dosages accumulated from 0 - 6, 6 - 12, 12 - 24, 
and 24 - 48 hours after burst time. As a secondary output, we have 
available the many thousand of r,h combinations, their locations on the 
ground, and their associated radiation contributions. This means, in 
effect, that any other set of grid points can be utilized in conjunction 
with this data and their radiation levels calculated. 


After the winds are fed into the computer, total machine time for a 
problem of this type is approximately one hour. Since plotting and 
analyzing the results takes about one-half hour, we can estimate that 
the total time from problem to pattern is something under two hours. 

doe 

Let us now examine some results of these calculations. The f3 
problem that was tried on the machine was the CASTLE-Bravo Shot. This 
was essentially our test of the machine program in that this particular 
case had been done by the original fend computing method reported in Rand 
Report R-265-AEC. An attempt was made, therefore, to duplicate the input 
data used in the initial calculation. In essence, then, the pattern 
computed by the machine was as wrong as the pattern in R- 265 -AEC in that 
it, too, used only the H-hour winds from the Curtiss (the ship near 
Bikini). This machine-computed pattern for integrated dosage to H / 24 
hours is presented in Figure 4. As was hoped, the machine computation 
duplicated the hand computation surprisingly well. In fact, due to the 
greater accuracy of the fall velocities input into the electronic com¬ 
puter and the greater detail possible with the machine, the computed 
pattern duplicates very closely the observed pattern in the vicinity of 
Bikini where the H-hour Curtiss winds may be assumed to apply with some 
degree of confidence. This is illustrated in Figure 5 which presents 
the observed dose rates in the Bikini Lagoon and the pattern computed by 
the above method. It should be pointed out that the computed pattern 
does not fit the far out fall-out observed from CASTLE-Bravo. This is 
to be expected in that one should not hope to duplicate a pattern that 
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took 17 - 24 hours to deposit on the ground stretching over a distance 
200 - 300 miles downwind from ground zero utilizing one set of winds for 
a particular time and place. As was pointed out in Dr. Eapp's paper, 
however, when this variation of wind in time and space is taken into - 
account, the Rand model does, indeed, produce an adequate reproduction 
of the observed Bravo fall-out. 

Figures 6 and 7 are the 48-hour dosage patterns for respectively a 
50 MT bomb under condition "A 11 and a 1 MT bomb under condition "B" . Con¬ 
ditions "A" and "B" are the wind conditions established by AFSWP for test 
purposes. Since these will be discussed later in this report, they are 
presented here just as a further illustration of machine computed pat¬ 
terns. It should be pointed out, that following the suggestion of Dr. 

F. Henriques of Technical Operations, we have gone over to the assumption 
that the concentration of activity in the mushroom decreases exponentially 
with altitude. This assumption has been used in all our calculations for 
conditions "A" and "B". One final point should be made. As was indicated 
above, to be truly representative, the computation must take into account 
the variation of wind in time.and space. Although this capability las not 
been utilized in the past, the 701 program is flexible enough to take 
this into account with a very small modification. Very simply, sets of 
winds representative of the variation in time and space are placed in 
the machine memory (these variations being in terms of finite increments 
of time and x and y distance) in solving equations (l) and (2). Then, 
for a given r and h, the machine examines the step-wise buildup of x, y 
and t over the summation, and switches to the appropriate wind set at 
the proper time and distance. It is hoped that this added refinement 
will be utilized in the future wherever feasible. 
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PREDICTION OF DOSE-RATE AMD DOSAGE CONTOURS AS FUNCTIONS OF YIELD 
AND METEOROLOGICAL CONDITIOtS, U.S. WEATHER BUREAU METHOD 

Mr. K. M. Nagler 
U.S. Weather Bureau 


Introduction 

Under contract to the Santa Fe Operations Office of the Atomic 
Energy Commission, the Special Projects Section of the Weather Bureau 
has attempted to devise a method of forecasting radioactive fall-out 
far tower bursts at the Nevada Test Site. In the past, fall-out fore¬ 
casts have been made by various approximate empirical methods. It 
was felt, however, that enough data had accumulated during the TUMBLER- 
SNAPPER and UPSHOT-KNOTHOLE tests to provide a basis for a more logical 
forecasting scheme. Specifically, we had useable amounts of observa¬ 
tional data from nine 300-foot tower shots, ranging from about 11 to 
about 52 kilotons in yield. 



The first phase of our work involved the preparation of maps 
of the observed deposited radioactivity. We have used primarily the 
ground monitoring data and have considered the aircraft estimated of 
radioactivity at the ground merely as approximate information which can 
provide a rough independent check on our method of determining fall-out. 


For roost of the data outside the test site itself, the monitors 
drove along the roads taking readings at intervals of 1-5 miles. 

Where only a single monitoring run was made along a road, we accepted 
its reports; but for most roads, there were several runs by different 
monitors, and the concentrations which we used on the map3 of observed 
radioactivity were a compromise - often a difficult one to make. The 
discrepancies between two or more monitoring nans along the same road, 
may be due to several causes: redeposition of debris that originally 
landed elsewhere, removal of debris by wind, instrumental errors, a 
different rate of decay than that assumed, errors in positioning, or 
to various other human errors. Figure 1 shows the results of several 
monitor runs along one road, following the ninth UPSHOT-KNOTHOLE burst. 
This illu strates many of the uncertainties which are encountered in 
using the observational data. Two runs give similar peak values but 
differ markedly to the left of the peaks, probably because the fall-out 
was not complete when the one run was made. A third peak suggests 
much lower concentration and is displaced from the otter two. The two 
partial runs obviously do not fit the others very well. One seemed to 
indicate that the reported positioning was wrong, perhaps due to a 
faulty speedometer on the monitoring automobile. 
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VARIATIONS IN OBSERVED ACTIVITY FROM MONITOR RUNS 



Figure 









This series of monitor runs is somewhat more confusing than mo3t, 
but not as confusing as some. It must be remembered that these obser¬ 
vations of deposited activity were not made with the care and accuracy 
desirable in a scientific experiment. They were made to determine 
whether there was any radiological hazard. 

At any rate, for each tower burst in the TUMBLER-SNAPPER and 
UPSHOT-KNOTHOLE series, the observations we selected were extrapolated 
to a common time (12 horn's after the burst) and were plotted on maps 
in units of milliroentgens per hour. 

Fall-out Plots 

The second phase of our work called for the preparation of fall¬ 
out plots. A fall-out plot or radex plot is an analysis vhich shows 
vhere particles of various sizes falling from various heights will 
reach the ground, based on the winds and on the assumed rates of fall 
of particles of the various sizes. A fall-out plot may also include 
an analysis of the times at which the different particles reach the 
ground. 

Figure 2 gives an example of a fall-out plot for UPSHOT-KNOTHOLE-9. 
Each curving line shews the estimates of where particles of a certain 
size should fall, and each radial line indicates the line along which 
particles originating at a certain height in the atomic cloud are 
deposited. A line indicating constant time of reaching the ground 
would pass through large size, hi£+i altitude areas to small size, lower 
altitude areas. 

Our fall-out plots differ from the type used previously at Nevada 
in two basic ways. In the first place, rather than a uniform rate of 
fall, we used an aerodynamic fall velocity, derived by means of the 
formulae given in the Rand report on Project Aureole. For the various 
altitudes, we used air temperatures and densities obtain®! from a mean 
sounding during the UPSHOT-KNOTHOLE series. As an example of the effect 
of the use of aerodynamic fall velocity, a particle in the 150-200 
micron rarge falls only about 7/10 as fast at 10,000 feet as it does 
at 40,000 feet. DOE ARCHIVES 

The second basic difference from most earlier fall-out plots is 
that in computing the path of a falling particle, we used analyses of 
the space and time variability of the wind, net merely the wind reported 
at the time and place of the burst. In brief, we prepared map® showing 
the mean winds for each 5,000-foot layer for each three-hour time when 
observations were available. On the maps each individual wind obser¬ 
vation (for a 5,000-foot layer) was plotted at the point over which the 
balloon actually was when it passed through that layer, rather than at 






the station from which the balloon was launched. Also, the original 
drift of the rising atomic cloud was taken into account in determining 
the starting point for the various elevations. 

Initially trajectories of particles of six sizes ranging from 60 
to 833 microns were drawn, starting from each 5,000-foot level - 
10,000, 15,000, 20,000 feet, etc. - and from the reported top of the 
cloud. Using the space and time changes in wind along the descending 
trajectory, particles of each size were tracked to the position on the 
ground or to the 5,000-foot surface (MSL), which with a few exceptions 
was close to the height of the terrain. 

For this initial fall-out plot we interpolated to get convenient 
intervals of height and particle size, as shown in Figure 2. Our 
height lines are at 2500-foot intervals, starting with the height of 
the top of the cloud. The particle size interval was chosen such that 
the time of fall from 35,000 feet to 5,000 feet is in multiples of an 
hour. Thus, a particle 291 microns in diameter takes one hour to fall 
this distance, a 162 micron particle takes 2 hours, etc. 

We took considerable care in preparing fall-out plots for the 
various tower bursts; nevertheless there are still several sources of 
error. 

The trajectories of the falling particle may be in error, due to 
errors in the wind observations and the analysis. The error in the 
location of a point is probably less than 10^ of the distance away 
from ground zero. However, the relative error between adjacent points 
on the fall-out plot is much less. 

A further possible source of error is that the rate of fall may 
be incorrect because the particles are not spherical, have a density 
different from the 2.5 assumed, or follow a different settling law 
than the one assumed. However, as far as we have been able to deter¬ 
mine, the radioactive particles produced by a tower burst are nearly 
spherical, particularly those smaller than 200 microns, which contri¬ 
bute mainly to the off-site fall-out. Moreover, data for the few cases 
for which times of arrival were available agreed with our fall-out 
calculations. Also, for the second UPSHOT-KNOTHOLE burst, there are 
particle size observations that seem to verify the aerodynamic law 
rather than Stoke*s law, or uniform descent. 

Construction of a Model &OE ARCHIVES 

Equipped with naps of the observed radioactivity along the roads 
and with fall-out plots for the various shots, we came to our basic 
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problem, namely to determine the distribution of the debris in the 
original clouds which would give rise to the ground radiological data. 

An explanation of our procedures is in order, however, before we 
discuss our model of the distribution of debris in an atomic cloud. 

This fall-out plot, strictly speaking, represents the position where 
various sized particles from various heights would reach ground, if 
they originated along the central axis of the cloud or in other words, 
if the cloud were very narrow. We know, however, that the cloud does 
have a finite width. We have assumed a cloud to be cylindrical and 
five mil es in diameter at all levels and that the diffusive spread of 
the cloud over the distances here being considered is negligible com¬ 
pared to the spread by wind shear. As will be pointed out later, the 
exact diameter chosen is not very important to our computations, so 
that our use of the same width for both stem and mushroom can be 
justified. 

Let us examine the importance of using a cloud of finite cross- 
section. If the cloud were not of finite width but were a line source, 
then each height and particle size in the cloud would correspond uniquely 
to one point on the ground. With a finite cloud diameter, however, each 
point on the ground could be affected iy a range of particle sizes and 
heights. 

We take the finite dimensions of the cloud into account in a 
fashion similar to that of the Rand and other groups. We first cen¬ 
ter a circle of the diameter assumed for the cloud at the point on the 
fall-out plot for vhich the fall-out is to be computed. The particle 
sizes and initial heights vhich contribute to the activity at the 
point are those that lie within the circle. 

In our model the amount of radioactivity is assumed to be uniform 
throughout each 2500-foot layer for each particle size range - a, b, 
c, etc. The mechanical computations of fall-out at a given point are 
of three types: first, where the circle lies entirely within one of 
the polygons; second, where the circle contains areas from several 
polygons; and third, where some or all of the polygons become squeezed 
into lines or nearly so - that is, the case of little or no directional 

shear# &QE ARCHIVES 

In the first case, if we know the concentration in appropriate 
units for each layer and size range, and if the circle lies within the 
sector, we merely need the ratio of the area of the initial cloud to 
the area of the sector, as shown by the upper point in Figure 2. This 
fraction multiplied by the amcunt of radioactivity assigned to the 
sector will give an estimate of the concentration at the point. 
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In the second case, we must find the ratio of each part of the 
circle to the area in which it lies. Each of these ratios or fractions 
is multiplied by the amount of radioactivity assigned to the appropriate 
height and particle size range. The sum of these products gives the 
estimate of the radioactivity at the center of the circle. The third 
case - little or no directional shear - will be discussed later. 

Here it should be pointed out that the assumed area of the cloud 
has no effect on the concentration as we would compute it at a point 
when the circle lies entirely within the bounds of one polygon. And, 
when the circle covers parts of several sectors, the assumed area is 
important only when the gradient of radioactivity from one sector to 
another is large. 

The procedure I have just described is the one *hich would be used 
when the amount of radioactivity in each height layer and particle size 
range is known. However, it is this information i*ich wa3 sought. 

To get this information, we first obtained the activities which 
would best fit individually the seven cases with best observational 
data. Each of these sets of numbers had to be a compromise, as no 
reasonable model could fit every observation well, and the various 
sets differed widely. Then, the main problem was to prepare a com¬ 
posite model which would give a best fit to all of the cases. 

First, we considered all cases relative to the top of the cloud. 
Second, to scale for yield, we assumed that the amount of activity is 
proportional to the area of the ground intersected by the fireball. 

(This scaling method is a simple one and one which seemed to work 
reasonably well. Further tests will be made with other scaling methods 
which are physically more acceptable. Changes in the activity assigned 
to various ranges of height and particle size may result from changes 
in the scaling method. The revised values, together with any refine¬ 
ments or simplifications in our procedures, will be given in a forth¬ 
coming report.) 

DOE ARCHIVES! 

When we completed our initial composite model, we used it to com¬ 
pute a grid of values - one for each sector on the fall-out plot, for 
each case. Then, considering how patterns and activities fitted the 
various cases, we continued to change the model to give a better fit. 
Figure 3 shows the resulting model, which gives the average activity 
in the various height and particle size ranges for a typical 20 kiloton 
bomb. In this connection, it should be pointed out that the total 
deposited radioactivity is not dependent solely on the yield, but also 
on the nature of the bomb itself. We have considered this effect only 
qualitatively, but feel that it may be an important one. 
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The model assumes a cloud top at 40,000 feet and a mushroom 
15,000 thick. It will be noted that a little more than half of the 
activity in the particle sizes greater than 58.3 microns occurs in the 
stem (7760 compared to 6680 units). Most of this, however, reaches 
ground within about 10 miles of the burst site, much of it within 
three or four miles, particularly that activity on large particles. 

It must be understood that these numbers are only approximate. 

There are enough observational uncertainties and perhaps enough real 
differences from one atomic cloud to another, so that errors as large 
as a factor of two or three are to be expected in parts of the model 
when it is applied to individual cases. . 

The units are those which will yield milliroentgens per hour at 
12 hours - if a cloud of 5 mile diameter is assumed. However, these 
units can be converted to units which will give the lifetime dosages 
(or dosages for any other period), 3ince for each initial height and 
particle size the tine of arrival at the ground can be calculated. 

In Figure 4 the units are for lifetime dose, in roentgens. 

At any rate, we have computed the activity on the ground accord¬ 
ing to this composite model and our assumptions, for eight of the 
bursts. Figure 5 shows the observed as well as the computed fall-out 
for one case, the ninth UPSHOT-KNOTHOLE burst. The data plotted along 
the road3 are a representative selection of the observed values of 
radioactivity in units of milliroentgens per hour at 12 hours after 
the burst. The data plotted near the dots are the computed values. 

The analysis showing activities of 5, 25, and 100 milliroentgens per 
hour at 12 hours is based on the computed values from our model. This, 
incidentally, is the case which our model fitted least well of all the 
cases, largely because of the abnormally great amount of activity 
observed over fairly wide areas. As can be seen from the figure, we 
are in error by a factor of about 3 over much of the area, but the 
general features of the computed and observed data are similar. 

Figure 6 shows a similar diagram for UPSHOT-KNOTHOLE 2. In this 
case, the model, again scaled for the appropriate yield, overestimated 
the activity as given by the observations by about a factor of 2. 

Also, the axis of computed peak concentrations differs from the observed 
peak by about seven cr eight mile3. This could, of course, be due to 
meteorological errors in preparing our fall-out plot or to an anomalous 
distribution of radioactivity in the cloud, with very high concen¬ 
trations existing near the very top of the mushroom. 

t»GE ARCHIVES 






AVERAGE DE5TRXBDTICW OF RADIOACTIVITY FROM A 20-KT, 300 FOOT TOWER BURST 
ASSIMSG A CLOUD DIAMETER OF 5 MILES. UNITS ARE THOSE WHICH GIVE 
LIFETn® DOSE AT THE GROUND, IN RCENTGENS. 
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Cases with Little Lateral Wind Shear 


I hav'e described our physical approach to faU-out prediction 
only insofar as it applies to wind fields in which there is consider¬ 
able lateral shear - that is, where the fall-out plot is spread out. 

Figure 7 shows a partial fall-out plot for a case (UPSHOT- 
KNOTHOLE 7) with little lateral spread. A complete fall-out plot 
would be practically meaningless. The "2 , s n indicate the position 
of the 88 micron line and you can imagine the elongated shapes of the 
polygons. The first and sixth UPSHOT-KNOTHOLE bursts had similar fall¬ 
out plpts with little lateral spread. 

In order to utilize the observational data from these bursts, and 
to test our composite model, we had to devise an alternate method of 
treating these cases. What we have done is to assume that there was 
no directional shear at all. Then, from the distances to the various 
selected sizes, we prepared a working graph shown in Figure S. In this 
graph, we have essentially compressed the fall-out plot into a line. 

One can see from this graph at what distance along this line a particle 
of a given size from any initial height will reach the ground. To 
determine the concentration at any point along the line, using our 
assumption of a cloud five miles in diameter, we first draw a vertical 
column five miles wide centered at the appropriate distance from ground 
zero. Then, to get the contribution, for example, of the top layer of 
the cloud and our "b” particle size range, we determine the fraction 
of the total area, (1, b) vhich lies within the column. In this case, 
it turns out to be about 17$. Thus, 17$ of the radioactivity in the 
height range 40,500 to 43,000 feet and in the particle size range 162 
to 291 microns in our model would contribute to the deposition at the 
selected point. To.estimate the total deposition at the point, we have 
to compute similar percentages for each sector cut by the column, mul¬ 
tiply each by the appropriate concentration from'the model, and sum up 
the products. This resulting sum must, of course, be multiplied by a 
scaling factor based on the yield of the bomb. Now, what we get by 
this method is the deposition at the appropriate point on the ground 
which would exist if the cloud were five miles wide and if there really 
were no change in direction of the wind with height. We know that 
there is some directional shear. Therefore, we would normally expect 
the value computed by this method to be too high. There is one 
exception, namely that a narrower cloud than the assumed five mile one 
exists together with no directional shear. Figure 9 shows the observed 
fall-out from UPSHOT-KNOTHOLE 7 and shows the ccraputations based on the 
method just described for several points where the fall-out line 
crossed roads along which observations were taken. The computed values 
are shown in italics and the observed values are plotted along the roads 
as before. In this case, it. can be seen that the computed values in 
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all cases exceed the observed peak values. As a test of our model, 
however, the observed peak values are not appropriate to use since we 
have assumed that the fall-out was oompressed into' a five mil e band as 
shown in the figure. Accordingly, we have integrated the activity along 
the various road3 and have converted to the concentration which would 
have existed if the cloud had been compressed to a five mile width. 

These ficticious peaks are shown in parentheses below the observations. 
It is apparent that in this case our model still gives excessive values 
but they are generally within a factor of 2 of the observed values. The 
position of the computed axis of high concentration fits the observed 
peak values very well. 


Proposed Procedure for Future Tests 


I have discussed the post mortem analysis of the fall-out from the 
300-foot tower bursts in previous test series. Now, let us consider 
hew we would use this information in forecasting for a future test. 


First of all, unless we had a convincing forecase to the contrary, 
we would assume that the cloud would rise to 40,000 ft., and that the 
mushroom would be 15,000 ft. thick. We would assume the amount of 
debris available far close-in fall-out to be proportional to the area 
of the fireball intersecting the ground (or we might use sons other 
scaling law if our present investigations indicate that it would work 
better). This can be estimated from the probable yield and from the 
tower height. Although in our post mortem investigation, we used 
trajectories and not merely the winds at the site, the time and effort 
to prepare such trajectories would not generally be warranted, consider¬ 
ing the uncertainty which would occur in forecast winds. Thus, we 
would base the fall-out plot on the forecast winds at the site only and 
the height lines would then be straight lines. 

Figure 10 presents a hypothetical fall-out plot showing consicfel^®? 
able directional shear in lower levels, but little directional shear 
in upper levels. From the lower part of the cloud, where there is 
considerable lateral spreading by the wind, the fall-out would be 
determined by the method which we described for the UPSHOT-KNOTHOLE 9 
shot. First, we would calculate the size of the various areas, possibly 
by planimetering, or by some more approximate but quicker method. Then 
where our assumed cloud lies entirely within a polygon, the ratio of 
the area of the 5-mile circle to the area of the polygon would yield 
the fraction of idle activity which would fall at a point within the 
area. Multiplying this fraction try the proper contribution from the 
model, we get an estimate of the fall-out to be deposited at any point 
who3e circle lies within that polygon. From a grid of these estimated 
depositions, generally one for the center of each sector, as isoline 
map should then be prepared. 
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From the higher altitudes of the cloud, where there is little 
lateral shear, the fall-out can more safely be predicted by the method 
assinning no directional shear, namely by assuming that the height lines 
are actually coincident. It must be remembered, of course, that the 
predicted values are apt to be on the high side unless the initial 
cloud has a diameter considerably less than five miles. At any rate, 
we would reduce the upper part of the fall-out plot to a graph similar 
to that previously shown, but only for the top 17,500 ft. (7 layers), 
to this graph, we would put columns five miles wide, centered at various 
distances from the burst site, and we would then determine the approxi¬ 
mate fractions of the concentrations in the various layers and particle 
size ranges. Then, after multiplying the fractions by the appropriate 
numbers in the model, we would have a series of forecast concentrations 
in the area to be contaminated by the upper part of the cloud. Now, 
from the predicted concentrations for both parts of the cloud, an 
analysis of forecast activity can be made. 

As with any forecast, the possibility of a shift in the wind, from 
that forecast and hence a different area of peak concentration must be 
allowed for as a safety factor and there may also be errors due to 
incorrect estimates of the cloud height end yield and, of course, 
there are bound to be errors in the details of the fall-out, 

Conclusions 

If anyone expects to find a model which can predict the fall-out 
from Nevada shots to within a factor of two consistently, we feel he 
is bound to be disappointed. The errors in the radiological observa¬ 
tions, meteorology and doubts about our various assumptions have already 
been dwelt upon at length. 

Our attempt to fit a single model to all the fall-out data, how¬ 
ever, uncovered certain interesting features. Two cases were very 
much different from all the others. Almost any reasonable composite 
model gives too much fall-out for TUMBLER-SNAPPER 8 and far too little 
fall-out from UPSHOT-KNOTHOLE 9. If it were not ifor an attempt to 
adjust the model to fit UPSHOT-KNOTHOLE 9, all the other cases would 
verify much better. While there was heavy cloudiness and possibly soma 
rain following UPSHOT-KNOTHOLE 9, we did not feel justified in discard¬ 
ing this case. DOE ARCHIVES 

We have in rough form, maps of our computed and the observed 
fall-out, in the event that some of you would like to inspect them. 

The verification of our model by ground monitoring observations leads 
one to the conclusion that about half of the predictions verify by a 
factor of two and probably 80 or 90^ verify within a factor of three. 







PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS 
FUNCTIONS OF YIELD AND METEOROLOGICAL CONDITIONS, 

USNRDL METHOD 

PART I 

Mr. C. F. Ksanda 

U. S. Naval Radiological Defense Laboratory 


Introduction 

The Military Evaluations Group at NRDL first got into the business 
of predicting dose rate and dosage contours 3omevhat by default. The 
group was established about four years ago for the main purpose of 
analysing the effects of radiological contamination on military opera¬ 
tions . It had been assumed that the basic fall-out data would be 
available elsewhere—but we soon found out it wasn’t. We were faced 
with the dilemma of getting quantitative answers to problems without 
having even qualitative input data. 

Before we could really do very much, we needed at least sons kind 
of rough estimates of expected dose rates and doses for a variety of 
detonations. At the time, the underwater shot at CROSSROADS was the 
only significant contaminating shot that had been fired. Also we were 
concerned with several Navy problems, and the Navy is particularly 
interested in detonations in and around water. It was natural, then, 
that the first attempt at prediction was for shallow underwater bursts, 
and that the predicting method was based on CROSSROADS data. The first 
approach was to devise a method for scaling the CROSSROADS data to 
other yields. In principle, this was like the scaling of physical 
effects—crater formation, for example, where one predicts th® influence 
©f yield without necessarily knowing the exact mechanism involved. 

At the time of Operation JANGLE the scaling method was extended 
to land surface and underground detonations. The JANGLE results, far 
1 KT, showed considerable area3 of significant fall-out downwind beyond 
the blast damage region. Scaling methods indicated that the areas 
enclosed by a given dose rate contour would increase more rapidly with 
yield than the areas of blast damage. So it appeared that the fall-out 
problem would become increasingly more important with higher operational 
yield weapons . 

This scaling 

intended only to get from 1 KT to 20-100 KT. No attempt, was made to 
predict twists and distortions in the fall-out patterns which might be 
caused by wind shear. Instead the main object was to estimate areas 
contaminated to given levels, and the distances downwind to which these 
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process was necessarily quite idealized and was 




levels would extend. The area concept Is useful when one is consider¬ 
ing a target region which is large compared to the fall-out contours. 
The downwind distance is useful when one wants to estimate the maximum 
radiation level which might be encountered at any distance from some 
assumed detonation. 


While this kind of information has a number of uses, it was 
apparent that some kind of quantitative model of the fall-out process 
was needed, in order to predict dose rates and doses for a wide 
variety of situations. Some work was started on (the development of) 
such a model—which I will mentioned at the end of ny talk. After 
that Mr. McCampbell will describe the recent development of a fall-out 
model. 

Meanwhile, for planning purposes, we have been using the scaling 
methods—and since they still seem to be useful, I would like to give 
a brief review of their development. 



The general problem required extrapolating a very limited amount 
of fall-out data to other yields and to other detonation conditions—- 
particularly wind velocities, soil types, and detonation depths. 
Scaling represents a limited solution; one starts with a known fall¬ 
out pattern and extrapolates to other yields. 


The method of scaling fall-out contours from one yield to another 
was developed in the following way. Shortly after the detonation, we 
visualize a region in the air containing all the radioactive particles. 
They have just stopped rising, and this region we will call the radio¬ 
active cloud. In order to scale fall-out patterns from one yield to 
another, we make five basic assumptions about this cloud and the par¬ 
ticles in it. POE ARCHIVES 

First, for equivalent scaled depths of detonation, the radio¬ 
active clouds from different weapon yields are assumed to be geometri¬ 
cally similar. Eoth the dimensions and heights of corresponding volumes 
are assumed to vary as the yield to some power * . 


the yield. 


Second, the total activity is assumed to vary directly with 


Third, for a given soil, the size distribution of active 
particles is independent of yield. 

Fourth , the relative spatial distribution of particles of any 
given size within the cloud is independent of yield. 
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Fifth , the rate of fall, for particles of a given size, is 
independent of the yield. 

These assumptions are undoubtedly idealizations of the actual facts. 
Later, scat® known departures of the facts from the assumptions will be 
mentioned. It appears, however, if the scaling is not done over too 
large a range of yields, these assumptions should lead to reasonable 

results. 

One further assumption is made vhen field dose rate patterns are 
to be scaled, rather than activity concentrations on the ground—that 
the dose rate at any point is directly proportional to the activity 
per unit area (on the ground). 

The scaling method follows directly from these assumptions. It can 
be stated in this way: 

If all linear dimensions of a contour to be scaled are 
increased as the yield ratio to the power a , then the 
associated dose rate will vary as the yield ratio to the 
power 1 - 2 a . 


The best value to assign to the cloud scaling exponent® is 3omewhat 
questionable. However, the numerical results are not particularly 
sensitive to the exact value chosen; and in most calculations we have 
used a value of l/3. In this case, the scaling method simply requires 
that all linear dimensions of the contours, as well as the associates 
dose rates, be scaled as the cube root of the yield ratio. 

The restrictions on this scaling method are that it applies only 
for: equal wind speeds —or if shear is appreciable, equal wind veloc¬ 
ities at equivalent scaled heightsj equivalent scaled depths of 
detonation; and the same soil . 


DOE ARCHIV ES 

We have tried to make some tentative estimates of the influence 
of detonation conditions on the fall-out patterns for a given yield. 


Wind : quantitatively, tends to 'longate ideally-circular patterns 
for no wind. HE data indicates areas of contours of given fall-out 
concentration quite insensitive to wind speeds; this generalization, 
together with an empirical estimate of the variation in downwind dis¬ 
tance with idnd spaed, was used to control the wind correction. 


expected to influence size distribution of radio¬ 
active particles formed, and hence the resultant contamination patterns. 


Examination of radioactive particles at JANGLE indicates they are - m 

formed from fused soil. It appears therefore that HE tests, where 
the dust fall-out consists of original soil particles, cannot be used 
to estimate influence of soil type. 

Depth ; provided cloud dimensions can be estimated as a function 
of depth of detonation, it is possible to make extrapolations by , 

arguments similar to those used to develop the method for scaling with 
yield. The main difference is that in the basic assumptions, the total 
amount of activity is constant, and the cloud heights and dimensions 
are assumed to vary in some way with depth. If the ratio of cloud 
dimensions for two depths is denoted by ft , then it can be inferred 
that the dose rate will vary as 1/at distances scaled as ft . 

Departures from Idealized Scaling 

The scaling method was first applied, for lack of any other data, 
to dust fall-out measurements from high explosive tests at Dug way 
Proving Grounds. Dust concentration contours from 320 pound charges 
were scaled to those from charges as great as 320,000 pounds—or over 
a charge -weight ratio of one thousand. Generally the results were 
fairly good, the scaled areas being off in the worst case by a factor 
of'less than three. 

It is now possible to make a comparison between the JANGLE surface 
diet and CASTLE Bravo. In this case the yield ratio is about 10,000. 

The nature of the soil over ’-hich the detonations occurred is, of - * 

course, quite different, as were the meteorological conditions. 

For the 2000 r/hr at 1 hour contour, the Bravo downwind dE&BihAsRCKIVES 
is about twice the scaled JANGLE distance. For the 100 r/hr contour, 
the observed and scaled downwind distances are equal. In all cases, 
the Bravo contours are considerably broader than the scaled JANGLE con¬ 
tours—by an amount which depends on vhich estimate of the Bravo con¬ 
tours we look at. On the basis of AFSWP-507 and the preliminary report 
of Project 2«5a, the CASTLE contours are about 2 to 3 times as wide as 
the scaled JANGLE contours. The area enclosed by the contours are * 

roughly 4 to 3 times as great as the scaled JANGLE areas. On the 
basis of the revised 2,5a contours, the discrepancy is considerably 
greater. The 2000 r/hr contour is about 9 times as wide as the scaled 
JANGLE value; and the 100 r/hr contour about 5 times as wide. 

archives * 

Two possible reasons far discrepancies between scaled JANGLE and 
Bravo results have already been mentioned: the difference in soil and 
in meteorological conditions, especially the wind structure. These 
differences really represent a failure of the tests of conform to the 
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necessary conditions. There are, however, two additional reasons (at 
least) where the basic assumptions used to derive the scaling method 
are not exactly correct. First, the rates of fall of the particles 
are not entirely independent of yield. The particles fall faster at 
higher elevations; hence as the yield is increased and the cloud 
height increases, the average rate of fall far a particle of given 
size increases somewhat. Second, the clouds are not geometrically 
similar. Empirical evidence shows a progressive flattening of the 
cloud with increasing yield. Thus the 1 KT JANGLE cloud was nearly 
spherical, while the width/thickness ratio of the Bravo cloud was 
close to 7#1. Scaling measurements as the cube root of the yield, the 
height of the center of the Bravo cloud corresponds almost exactly to 
the scaled JANGLE height. However, the depth of the Bravo cloud is 
only about one-half the scaled JANGLE value, while the width is about 
3 times as great as the scaled JANGLE width. 

Refined Model 

Some approximate corrections could be made to account for the 
variation in cloud shape and rates of fall of particles, in order to 
improve the scaling method. But the only way to predict actual fall¬ 
out patterns for a variety of conditions is by way of a dynamic model 
of the fall-out process. 

Exploratory work on such a model was begun nearly two years ago, 
to see what kind of agreement with the JANGLE data might be obtained. 
The first attempt assumed a spherical cloud, uniform wind velocity, 
and contained a very rough estimate of the distribution of activity 
with particle size. Somewhat surprisingly, the dose rates predicted 
from this model looked fairly much like the observed values. 

Now it is apparent that any model of the kind that has been dis¬ 
cussed here contains quite a few parameters, to which some numbers 
have been assigned. These par acre ters are closely related, to the things 
which were mentioned earlier in the five basic assumptions used to 
develop the scaling method. 

First , the total activity 
straightforward calculation. 

Second , the cloud height and dimensions must be known. 
Generally these can be estimated from photographs and extrapolated with 
the help of some theories---although there is no guarantee that the 
radioactive cloud is actually coincident with the visible cloud. 


S^OE ARCHIVER 

must be known—a relatively 




Third , activity in the cloud must be known as a function of 
the particle size. 

Fourth , distribution of the particles in space must be 

known. 

Fifth , rates of fall of the particles must be known. 

Regarding these last three points, there is some controversy. 

Rates of fall depend on the properties of the atmosphere as 
well as those of the particles. Actually, particle properties vary 
considerably—those from Nevada being sometimes spherical, sometimes 
very irregular, and showing density variations of about a factor of 2, 
mainly because of varying amounts of included air bubbles. In our 
numerical work, we have assumed that all the particles are irregular 
and of uniform density. 


Distribution in space is clearly not well known. Here we 
have assumed that practically all the activity is in the mushroom 


cloud and is uniformly distributed. 


Finally, the relation of activity to particle size is not 
well known either. Rather than try to guess at this, or adjust the 
distribution to give the best agreement between the predictions of the 
model and experimental data, we tried to nake an independent evaluation 
from some particle studies at JANGLE. 


A fairly large number of fall-out samples were separated into 
various size fractions at NRDL, and the relative activity of each size 
fraction was measured. The data were not sufficient to draw activity 
concentration contours for each size fraction, which i3 what we would 
have liked, but it was possible, for the underground shot, to draw 
curves of activity concentration as a function of distance in the down¬ 
wind direction. By integrating these curves, it was then possible to 
estimate the percentage of activity in each size range which was 
initially present in a unit width slice of the cloud. 


This distribution tunned out to be a normal distribution with 
respect to the logarithm of the particle diameters—and this is the 
distribution we are using in our numerical calculations, jjqj. ARCHIVES 

This gives soma of the background of the development of our fall¬ 
out model, which Mr, McCampbell will describe. 
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PART II 


Mr. James M. McCampbell 
U« S. Naval Radiological Defense Laboratory 


During this Symposium, it has been amply demonstrated that the 
problem of devising a mathematical model of fall-out consists of two 
major parts; first, tracing the trajectories of the fall-out particles 
to determine where they land, and second, summirg all the activity that 
lands at any given point on the ground. These two basic ideas will be 
developed in their mathematical form. 

Concerning the first part of the problem, we consider a particle 
falling through the atmosphere with a velocity u and drifting in a -wind 
of velocity v as shown in Figure 1. During a small interval of tine 
the particle drifts along the vector dD whose components dx,dy, and dz 
refer to a convenient coordinate system. The nagnitude of the time 
interval equals the distance the particle travels divided by its rate. 

- dz 

In the vertical direction, dt = u. The distance traveled in the 
direction of the wind is d/ie v dt and by substitution, equals v 
But dx is a component of dp , hence dx-=» dp sin 0 where 0 is the 
angle between the y-axis and the direction toward which the wind is 
blowing. Theta is 180° from the conventional wind angle. 

As the particle traverses successive layers of the atmosphere, 
its path is defined by the sum of the vectors dD. Therefore, the total 
drift in the x direction is the sum of the x components; that is, Djx 
equals the integral from the surface to the altitude of origin of 
v sin 8 divided by u. A similar expression involving cos 8 will of 
course give the y component of the total drift. The parentheses indi¬ 
cate that the wind speed and direction vary with altitude and that the 
particle falling rate is also a function of the altitude. In addition, 
the falling speed of the particle is a function of its size, 1. 

A convenient coordinate system is the one with the origin at 
ground zero and the y-axis directed north as in Figure 2. Consider 
particles of a given size originating in the infinitesimal volume dV 
at the point (x,y,z)» These particles drift along some irregular path 
indicated by the heavy dotted line, and land at the point (A,B). The 
final displacement of the particles from the origin is the sum of their 
initial displacement in the cloud and the displacement they experience 
while drifting with the wind. The x-component of the final displace¬ 
ment is A, and it equals x plus D^. Transposing the terns gives 

x = A - D x and by a similar derivation y = B - Ey 









How what is the meaning of these two equations? For a given point 
(A,B) and a given wind structure, any value of z determines Dx and Dy, 
and in turn determines a cofrespondii^ x and y. Thus, these equations 
are the parametric expression of a curve in space. It is along this 
curve that particles of a given size can possibly originate, traverse 
the wind structure, and land at the point (A,B). Mathematically, this 
locus of origins has an infinite extent, but physically, only a segment 
has significance; that is, that segment which lies within the cloud. 

We must therefore define the boundary of the cloud. 

As shown in Figure 3» an ellipsoidal cloud is assumed. Its semi¬ 
axes are a and b and its center is directly above the origin at the 
altitude H. Such a model will accomodate the very flat clouds from 
high yield weapons, and will degenerate into a sphere for clouds from 
low yield weapons. The locus of origins is shown penetrating the cloud 
at the altitudes zp and z 2 . The values for zp and z 2 are different for 
each particle size. The points of intersection aire found by the simul¬ 
taneous solution of the equation of the curve and the equation of the 
cloud boundary. For a spheriod, the latter is 

r Z—Hi* r X-i* r Y-i 2 

L“r~] + Lt 3 + hd = 1 

Our equations for the curve were in terms of x and y, so substitution 
gives 
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The values of zp and z 2 are implicitly given by the last equation. 
Any point in space which satisfies this equation must lie on the curve 
and also on the cloud surface; it is therefore a point of intersection. 
These points may occur in real pairs, imaginary pairs, or singly. 

Real pairs indicate penetration of the cloud, imaginary pairs indicate 
that the curve does not touch the cloud, and single points indicate 
that the curve is tangent to the cloud. The physical meaning of tbs 
imaginary case is that particles corresponding to that curve can not 
reach the point in question. The curve may penetrate the doud, fold 
back and penetrate it again; and perhaps again. However, for most wind 
conditions, the simple case as shown in the figure will probably pre¬ 
vail. Whether or not it does is immaterial—as will be shown later. 

Because the doud height, shape, and size vary with the yield, 
the cloud parameters H, a, and b will take on values depending upon 
the yield. Therefore, substituting emperical or theoretical relations 
between these parameters and the yield would make the last equation 
applicable far all yields. 
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(I)dl dxdydz 



RU,y) = kj P*(i)e(i)[z 2 (i)- Zl (i)] 
















We must now turn our attention to the second part of the problem, 
that is, to summing the activity that lands at a given point. In 
Figure 4, P is the number of particles per unit volume in the cloud. 

If<£ ( 1 ) is the size distribution of the particles then P<£ (l)dl is the 
number of particles per unit volume of sizes 1 to 1 plus dl. Multi¬ 
plying by dxdydx gives the number of particles of that size in the 
volume dV. 

Let £ (1) denote the radioactivity of a particle of size 1, assum¬ 
ing that the radioactivity of particles is in some way related to their 
size. Then multiplying by £ (1) gives the radioactivity of the par¬ 
ticles of size 1 in dV. But it was previously shown that these par¬ 
ticles were deposited at the point (A,B). Assuming that there is no 
dispersion of the particles in their flight, they will be spread, over 
an infinitesimal area dxdy. Therefore, the activity per unit area is 
found by dividing by cbcdy as in the last equation in the figure. It 
oust be pointed out that da^ is a second order infinitesimal; the two 
remaining differentials are dl and dz. It is apparent that subsequent 
integration must be performed over the altitudes and the particle sizes. 

For particles of each size, the integration with respect to z must 
be confined to the applicable altitudes. In Figure 5 we see that a]_, 
the activity due to particles of size 1 , is found by summing over the 
range z^ to Z 2 * ^ our curve of origins had penetrated the cloud twice, 
we would add another term with limits Z 3 and z^. Having integrated 
with respect to altitude, we must now sum over the particle sizes as 
shown in the middle equation, where the previous integration has been 
carried out. 

Finally, the radiation field produced above a contaminated surface 
is proportional to the activity on the surface. Therefore, the radia¬ 
tion field R is obtained by multiplying the activity per unit area by 
the proportionality factor K. The value of K includes the effect of 
shielding by rough terrain, the mean energy of the gamma rays, and the 
height above the surface that the radiation intensity is desired. It 
should be pointed out that the product^ ( 1 ) £ ( 1 ) is the activity 
distribution. It is sufficient to know the product of the two 
functions; they need not be known separately. 

It turns out that solving the problem in this form 
impossible. However, it can be readily solved by taking finite incre¬ 
ments of particle size and altitude. When formulated in terms of 
finite intervals, it is ideally suited to machine computation. A com¬ 
putational scheme based on this theory is currently being used in 
computing fall-out patterns on a Univac. The machine work is being 
done by the Applied Mathematics Laboratory of the Bureau of Ships. 
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You will notice that no mention has been made of a very important 
variable, namely, the time elapsed since the explosion. This variable 
enters the problem in four specific ways: 

In the decay of radioactivity. 

In the build-up of the activity at any poiiit and the conse¬ 
quent growth of the fall-out pattern. 

In the changes of wind structure that occur during the fall¬ 
out process, and finally. 

In the change of the mean gamma energy with tins. 

Recent extensions of the present theory permit a machine solution 
including the variable time in all four aspects. The output table from 
the machine could provide the following information: 

The coordinates of the point in question. 

The radiation intensity at the point at any time. 

The dosage accumulated at the point up to any time. 

The time in question. 

The time that fall-out began to arrive at the point, and. 

The time that fall-out ceased at the point. 

If the preliminary computations now being made confirm the model, 
there will then be available a useful tool for studying fall-out 
problems. In particular, investigations could be carried out in such 
subjects as the: 

Scaling of radiation fields. 

The effect of different wind shears, and, 

The effect of different particle and activity distributions. 
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PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS 
FUNCTIONS OF YIELD AND METEOROLOGICAL CONDITIONS; 
PREDICTED CONTOURS FOR THE CASTLE DETONATIONS 

E. A. Schuert 

U. S. Naval Radiological Defense Laboratory 


At Operation CASTLE Project 2.?a had as a primary responsibility 
the definition of the fall-out patterns resulting from the 6 deto¬ 
nations. A great deal of analytical work has been expended in the 
reconstruction of the contours that defined the primary fall-cut 
pattern resulting from Shot 1. Since the analysis that was employed 
combined experimental data with particle trajectory evaluations, much 
of the work is similar to the model developments presented here this 
week. I will discuss the assumptions made and the techniques employed 
in this development as well as the relevant experimental data obtained 
and how this data fits the resulting fall-out contours. 

The data available after Shot 1 included the measured fall-out in 
and about the Bikini lagoon, the negative data from the free floating 
sea stations west of the atoll, the outer island survey measurements 
made by Dr. Scoville’s team, the observed time of arrival of fall-out 
at Bikini and on the outer atolls, as well as the measured radioactive 
particle siae, particle density, and particle configuration data for 
both the Bikini area and the outer islands. There was also data ob¬ 
tained on the time of arrival of fall-out particulate at Bikini. 


An attempt to reconstruct the fall-out contours solely from the 
measured gamma field data offered the analyst a great deal of latitude, 
for there existed no information that would control the lateral width 
of the fall-out pattern. Because of this problem an attempt was made 
to define the area of fall-out by use of particle trajectory analysis. 
Since we were primarily interested in singularly defining the Shot 1 
pattern, data from this detonation was used in defining the trajectory 
parameters. The parameters necessary to apply the trajectory analysis 


were: 


A source of activily. 
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The initial spacial distribution of the source material. 


A mechanism for distributing the source material. 


Some of the above parameters were supported by experimental data; 
others required basic assumptions. The assumed parameters were: 
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Figure 1 

General Area of Fall-out, CASTLE-Bravo 








The spacial distribution of the source material was described 
by the visible dimensions of the cloud and stem at time 0 + 10 minutes* 
This assumption was further simplified by describing the Shot 1 stem 
and cloud as cylinders. 

The distribution of the source material was homogeneous 
throughout the cloud and stem with all particle sizes at all locations. 
The particle trajectories cloud be defined by their rate of 
fall under the influence of the existing wind speeds and directions. 

The experimental data obtained from analysis of the fall-out 
material as well as other physical measurements was used as follows: 

The measured physical characteristics of the particles were 
used to determine their rate of fall and to describe the particle size 
distribution that was actively contributing to the primary fall-out. 

The cloud dimensions as measured by cloud photography were 
used to describe the physical dimensions of the source. 

The observed meteorological data both on windspeed and 
direction and temperature distribution as obtained by the Task Force 
weather central were used to calculate the particle trajectories. 

The trajectory analysis was then defined on both assumed and 
measured parameters. 

In determining the rate of fall of the particles, a size range of 
from 2000^-to 25>* diameter was used. The shape of the particles had 
been experimentally determined to be almost entirely irregular in 
shape. The apparent density of 70 particles was measured resulting in 
a mean density of 2,36 g/crP with a standard deviation of 9 par cent. 

The air temperature used in determining the viscosity of the atmosphere 
with respect to height was taken from that measured at Bikini Just 
prior to shot time. After determining the falling speeds of the 
particles over the stated size range using the above experimentally 
determined parameters and applying aerodynamic particle settling rates, 
a table of average falling speeds was constructed for each 5000 foot 
increment of altitude to a height of 100,000 feet. 
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With the above falling rates, 231 particle trajectories were ****< 
computed. These computations were done manually, for at that time 
application to machine analysis was not available. Particle sizes of 
2000, 1500, 1000, 750, 500, 375, 250, 200, 150, 100, 75, 50, and 25ft 
in diameter were placed at 5000 foot intervals to an elevation of 
100,000 feet and their intersections with the surface of the earth 
were determined by following each particle as affected by the wind 
structure. In constructing these trajectories the changes in the wind 
pattern both with respect to time and location of the particles were 
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Fall-out Contours, CASTLB-Bravo 











considered. This analysis resulted in the plot shown in Figure 1. 

Here there is a picture of the general area of fall-out based on the 
stated assumptions. The outer dashed line was constructed by expand¬ 
ing each point by the cloud or stem radius from which the particle 
originated. The dimensions used were 66 nautical miles cloud diameter 
and 6.6 nautical miles stem diameter, with the cloud base at 60,000 
feet. The dashed arcs show the maximum distance the stated particle 
size could travel. Because of the wind structure existing at shot time 
the area of fall-out had a rather well defined axis along which the 
larger particles were closely grouped. This axis is shown by the heavy 
solid line in the figure* Because this axis did not change noticeably 
by varying the particle size distribution, it was reasonable to repre¬ 
sent it as an axis of symmetry of the fall-out pattern* 

The reconstruction of the fall-out contours for Shot 1 were then 
made, based solely on the directive axis of fall-out as determined 
above and the measured gamma field readings. The analysis resulted in 
the fall-out contours as shown in Figure 2, These contours show gamma 
activity in r/hr at 1 hr. 

The most significant result of this analysis was the unusual width 
of the fall-cut pattern. For comparison, the surface Jangle fall-out 
pattern was scaled to IS megatons using the cube root relationship. 

The most significant observation was again the great difference in the 
configuration of the two patterns a A reason for this discrepancy was 
discussed by Mr. McCampibell earlier this morning. 

This analysis based on a fall-out axis and the measured gsnma 
fields, created a pattern that fits well with the lateral limit of the 
pattern as defined by the expansion of the trajectory points as pre¬ 
viously shorn in Figure 1* 

Two material balances were made on the reconstructed pattern* 

The first material balance based on calculated parameters as described 
by Dr. Werner earlier this week result,ed in from £0 to 70 par cent of 
the device accounted for within the 100 r/hr at 1 hr. contour. The 
second analysis based on extrapolation of a radio-chemical analysis of 
the fall-out material collected at one location, indicated 30 per cent 
of the device accounted for within the 100 r/hr at 1 hr. contour. 

In order to shew pictcrially the growth of the fall-out" 
with time and the actual existing levels of activity, I have prepared 
five figures (Figures 3”7), showing the existing gamma field at 1, 3, 

6, 12 and 18 hours. 
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Figura 3 

H+l Hour Fall-out Contours, CASTLE-Bravo 








Figuro 4 

Bf3 Hours Fall-out Contours, CASTLE-Bravo 














I MAXIMUM EXTENT OF 
PATTERN AT H + 12 HRS 



Contours, CASTLE-Bravo 
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Figure 7 

H4-1S Hours Fall-out Contours, CASTLE-Bravo 





I would like to discuss the experimental data collected that is 
relevant to the analysis of the fall-out based on trajectory 
evaluations* Figure 8 shows the computed time of arrival of fall-out 
based on the 231 trajectories evaluated* It indicates an arrival time 
of 8,2? hours at Rogerik where the measured time of arrival was 0 
hours* It shows an arrival time of l?*it hours at a distance of 302 
nautical miles when the observed time was 18 hours. The analysis 
appears to satisfy the time of arrival parameter very well. 

Figure 9 shows the particle size distribution obtained experi¬ 
mentally from fall-out at Bikini atoll. This is a histogram of some 
7000 radioactive particles. The trajectory analysis agrees well with 
deposition over this size range at Bikini atoll. The data on these 
particles analyzed with respect to their time of arrival shows positive 
evidence that of the particles over 5CX3^in diameter ?0 to 90 per cent 
had to come from above 60,000 feet. 

Figure 10 shows the measured particle size distribution of the 
fall-out that arrived at the outer islands of Alinginaie, Rongelap 
and Utirik. The computed size distribution at Rongelap atoll based on 
the trajectory analysis indicates that particles from ?0 to 1?0 
diameter would arrive at this location and the 2$0 /x. particles would 
fall just 10 nautical miles north of the atoll. The measured distri¬ 
bution shows that only it particles arrived at Rongelap with diameters 
larger than 200 ft * 

The analysis also indicates no particles above 7? p. diameter 
could arrive at a distance of 302 nautical miles from ground zero 
where the measured distribution at Utirik had a geometric mean diameter 
of kS with only 3 particles being over 100 p.in diameter. 

The agreement here is surprisingly good and certainly suggestive 
of continuation of work along these lines. The major discrepancy in 
the analysis is the inability to account for fall-out at Utirik. This, 
may be resolved by the use of Dr. Palmers analysis of the wind data 
existing at and after shot time* 

In conclusion, it may be well worth discussing the fields in which 
serious work should be added both from an experimental as well as 
theoretical viewpoint. A portion of this work necessary to a more 
complete understanding of the phenomenology may be listed as follows: 

The relation of activity to particle size appears necessary 
to most model evaluations. At this time data is seriously lacking in 
this field, 
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Figure & 

Computed Time of Arrival of Fall-out 
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Little is known at this time on particle formation 
phenomena. Certainly more thought should be given to this field, 
especially with respect to understanding the complete mechanism of 
fall-out. 

The configuration of the radioactive cloud source as well 
as the particle size distribution within the source are at this time 
unknown. This knowledge seems paramount to any resolution of the 
basic concepts of the models presented at this symposium. 

There still remains the small particle anomaly where 
measured fall-out with respect to time indicates the very small 
particles arriving at times much too early to account for their 
trajectory by any falling rate equation. 
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PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS 
FUNCTIONS OF YIELD AND METEOROLOGICAL CONDITIONS, 

ORO AND TECHNICAL OPERATIONS, INC. METHOD 

F. C. Henriques 

Technical Operations, Incorporated 
Introduction 

Three months ago Technical Operations was requested by the Opera¬ 
tions Research Office of the Johns Hopkins University (ORO) to provide 
1 Hr-R/hr dose rate and two-day dose contamination patterns which would 
result from the surface detonation of a 15 MT thermonuclear weapon. 

These patterns were computed for the thirty wind conditions which re¬ 
sulted from considering wind conditions at six different locations along 
the Eastern Seaboard during five different but “typical 11 meteorological 
situations. Previous to initiation of these computations, conferences 
-were held with those groups at Sandia, NRDL and Rand that are developing 
fallout models and prediction procedures. Based upon these discussions, 
a "fallout" model was evolved and a "hand" computational procedure was 
developed which enabled the calculation of sixty fallout patterns in a 
one-calendar month with a reasonable number of man months. 

The same fallout model and procedual method were used to compute in 
a two-week period the thirty-six radiological contamination patterns 
given in this report. These thirty-six patterns result from the calcula¬ 
tion of dose rate and two-day dose contours for 1, V~> and 50 MT surface 
detonations for the six AFSWP wind conditions. 

The underlying constraints required for these confutations are 
markedly similar to those which have been described in the Symposium 
Report by Rand, the Signal Corps, and NRDL. The details of the computa¬ 
tional methods employed by the various groups do differj those recently- 
developed by Rand and NRDL are well suited to machine confutations, while 
those developed by the Signal Corps and Technical Operations have been 
specifically tailored to "hand" calculation procedures. It is not the 
purpose of this report to detail cur computational procedure; this is 
reported elsewhere. 

If the input assunftions used by the above-mentioned groups were 
identical, identical fallout prediction patterns should result. Unfortu¬ 
nately, since there is a lack of experimental data, assunftions must be 
made to? (1) the radioactive cloud spacial dimensions; (2) the size dis¬ 
tribution and fall times of entrained particles; (3) the fraction of 
activity associated with a given particle size; and (U) all the conversion 
factors which ultimately convert this activity once it reaches the ground 
to roentgen dose rate and biologically effective cumulative dose. Each 
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of the four above-mentioned groups use different sets of assumptions. 

* 

In Section I we outline the constraints -which are believed to be 
common to each of these four methods and list the assumptions inherent 
in the Technical Operations model. We then discuss the areas and shapes 
of the 100, 500 and 1500 fission product 1 hr-R/hr dose rate and two-day 
cumulative dose contours which result from the surface detonation of 1, 

15 and 50 MT weapons for the six AFSWP wind conditions. This is followed * 

by an analysis of the similarities and scaling relationships of both the 
areas and shapes of these contours. Finally we propose a rapid scaling 
procedure that is believed to be reasonably valid. 

In Section II we assume that the 1 hr-R/hr and 2-day dose roentgen 
readings at any ground location are known input data. We then derive 
the equations required to compute the biological effective dose from 
fission products and the contribution of U-239 and Np-239 thereto. The 
utility of the graphs and tables generated by these equations is illus¬ 
trated by solving three typical practical problems. 

The detailed justification, of the information stated in Sections I 
and II is available elsewhere (ORO report in press). 

Radioactive Fall-out Constraints. Assumptions and Results 

To avoid duplication, it is assumed that the Rand and Signal Corps 
papers which comprise a portion of the Symposium Report have been read. 

The constraints inherent to their procedures as well as our approach may 
first be listed. Some of these could be dropped but only at cost of more 
complicated calculation procedures. 


A. There i3 some instant ("instant of stabilization") at which 
the distribution of activity vs particle size and the spa¬ 
tial distribution of such particles is known. 

B. The range of particle sizes is such that pure diffusion of 
particles can be neglected and yet the particles as they 
fall will follow wind currents at each pertinent altitude 
with no important lag. 

C. The range of particle size is also such that for any two 
given particles the ratio of their fall times from various 
altitudes is constant (independent of altitude) to suffi¬ 
cient approximation. 

D. At the instant mentioned in A, the distribution 
vs particle size is independent of spatial position. 

E. At the instant mentioned in A, the spatial distribution of 
total activity is laterally uniform (within some area) at 
ary one altitude. 

F. There are no vertical winds, and the horizontal wind pat¬ 
tern vs altitude remains constant over a time sufficiently 
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long for all those particles destined to arrive at loca¬ 
tions of interest to complete their fall. 

These list only those constraints required to yield a rapid '•hand” 
calculation technique. Before discussing these, it is perhaps worth 
pointing out several items that need not be assumed to achieve compu¬ 
tation rapidity. The vertical distribution of activity in the cloud at 
the instant of stabilization can be arbitrary, and the cloud can have 
an arbitrary physical shape. Likewise, the distribution of activity vs 
particle size can have any shape ; the only restrictions here are on the 
range of sizes, which must satisfy B and C. 

Constraint A merely states that the problem is restricted to falling - 
particle calculations and does not concern the processes which establish 
the radioactive cloud. 

An inertial lag and diffusion analyses show that particles at least 
in the range 10 microns to 2000 microns satisfy Constraint B, The 
lower limit is certainly unnecessarily high, but appreciable raising of 
the upper limit would require more careful investigation. However, the 
effects of large particles in the model are rather insensitive to this 
constraint since they all fall close to ground zero in any case. 

The rather surprising property postulated in Constraint C is actu¬ 
ally satisfied quite well over the same range, 10 to 2000 microns. 

This is subject to direct calculation. The results of these calcula¬ 
tions are given in Tables 1 and 2. 

Constraint D is rather ad-hoc; it enables the problem to be broken 
into independent components. Until actual measurements at the "instant 
of stabilisation" are available, it must rest entirely on the argument 
that continual mixing occurs during the relatively brief period (10 
min) when the "boiling" cloud is rising and that in this way the chances 
of a given particle being located at a given position are independent 
of its size and of the activity which it carries. It is doubtful that 
this is literally true, but the argument seems intuitively reasonable 
enough to justify its use until data can be shown to seriously refute 
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Constraint E simplifies calculation to the extent that only geo¬ 
metric coincidence data need be taken off the plots used in calculation. 

It could be lifted at the cost of considerably more data-handling dur¬ 
ing calculation, but it does not affect the over-all scheme of calcula¬ 
tion. Very little data seem available on this point, but the few 
measurements which have been made with drone airplanes at Greenhouse 
tend to give confirmation. 
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CLOUD PROPERTIES "AT STABILIZATION 





















Constraint F is required to avoid plotting the map position of more 
than one particle size. Vertical winds are probably negligible in a 
great many cases. Constancy of the horizontal winds over long time 
periods is probably not too realistic. Provided that the meteorologi¬ 
cal data were sufficiently well known, Constraint F could be lifted at 
the cost of plotting in detail the final positions of all particle 
sizes from all altitudes. Subsequent calculation would not be affected 
however. 

We now proceed to discuss the Technical Operations data assumptions 
which result in our specific fall-out patterns. These assumptions are 
not specific to the method of calculation but must, of course, be given 
before a calculation can be performed. 

The cloud properties are shown in Figure 1. The plotted points are 
experimental surface shots where the radius is the maximum radius of 
each cloud. fHe^experimental points have been connected by rough curves 
which have been extrapolated in both directions. These extrapolations 
are based on the following arguments. 

For larger weapon yields, in lieu of data, it seems reasonable to 
assume that the radius will continue to grow at about the sans rate as 
before; thus the radius is extrapolated as wV 2. x positive steep 
atmospheric temperature gradient starting at about 120,000 ft (the 
“ozone absorption layer") can be expected to limit the cloud top to 
about this altitude. This is because the cloud must expand as it rises 
in order to maintain pressure balance with the surrounding air. Be¬ 
cause of the expansion, it cools more or less adiabatically. At lower 
altitudes, however, the surrounding air temperature also drops with 
pressure approximately adiabatically and the cloud, once much hotter 
than its surrounding air, remains hotter as it rises and thus continues 
to have bouyancy. This process reverses, however, at about 120,000 ft 
and the cloud should very rapidly lose boi^ancy and cores to thermal 
equilibrium with the atmosphere. The atmosphere exhibits another posi¬ 
tive temperature gradient in the neighborhood of U0,000 - 60,000 feet 
(the tropopause). However, clouds resulting from nuclear detonations 
in excess of 1 MT still have sufficient internal energy when they reach 
the tropopause to punch through it. The cloud top has therefore been 
assumed to have a ceiling of 120,000 ft. The data for the cloud bottom 
also appear to indicate a leveling off. This would be expected on the 
basis of the other two curves, for the mass of the original fireball is 
known to be nearly proportional to W. If this is true also of the cloud, 
the cloud bottom must level off aa shown at 60,000 ft. The data seem 
to support this conclusion. 

The extrapolation toward lower energies is based on data for the 
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surface shot of 1.2 KT which nearly fall on the straight-line extrapo¬ 
lations shown. In addition, these extrapolations yield a total cloud 
mass very nearly proportional to yield - assuming the mass-distribution 
within the cloud is the same as that in the atmosphere it replaces. 

An additional property of clouds which is required is their shape. 
Here data are quite sparse for the "late times 1 ' (about 10 min) when 
the cloud finally stabilizes. 

Early calculations were done assuming a cylindrical shape but this 
leads to rather low peaks of contamination* and does not agree with pho¬ 
tographs and descriptions scattered through the classified literature. 
The shape finally used was taken from a written description of the 
Castle-Bravo shot and is substantiated by Figure 2 which is a photo¬ 
graph** of Ivy-Kike shot (H.2 MT) taken essentially at the "instant of 
stabilization". The result (illustrated in Figure i*) is an inverted 
cone truncated so that the cloud bottom has a diameter lj.0$ of top diame¬ 
ter given by Figure 1. There is nothing at all sacred about this partia 
ular shape; any shape could be used in the calculations but a truncated 
cone seems to be the best possible assumption at the present state of 
knowledge, especially since the limiting action of the "ozone ceiling" 
must enlarge the cloud top with respect to the bottom. 

Another piece of required data is the distribution of radioactivity 
vs. particle size. Such measurements have been made only*** for the 
Jangle surface and underground shots of about 1 KT. These data (NKDL, 
private communication) are plotted in Figure 3 which shows for each 
particle size, plotted vertically on a logarithmic scale, the fraction 
of activity associated with particles of that size or less, plotted 

* This fact has also been recognized by the Rand group. For their 
calculations they are apparently now using artificial cloud radii 
which are about 60 $ of the reported cloud radii. 

** Observation of the motion picture of the 15 MT Castle-Bravo shot 
indicates the same shape at this late time. A still photograph 
is not available at this writing. 

*** Other measurements have given the size-distribution of particles 
with some activity, while still other measurements have yielded 
activTfy - versus particle size but unfortunately for different 
weapon tests. Neither alone is sufficient. The combination of 
these is what is required. 
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horizontally on a probit scale. The straight line plot means that the 
activity has a Gaussian normal distribution in the logarithm of particle 
diameter. The data in Figure 3 may he expected to hold for surface 
shots over most land areas since the intrinsic particle sizes of various 
soil-particles and concrete are not greatly different. The data may also 
be expected to hold for a wide range of weapon energies for the follow¬ 
ing reason. Within the fireball the conditions of formation of radio¬ 
active particles must be quite localized and thus depend only upon 
the energy density per unit mass of debris. The mass of soil blown out 
of the crater which becomes part of the fireball in a surface shot is 
proportional to weapon energy. Therefore the energy density per unit 
mass of soil within the fireball is independent of yield. This leads 
one to expect that the conditions of formation of the radioactive 
particles are independent of the weapon energy. 

It must be realized, however, that the distribution curve of Fig¬ 
ure 3 does not necessarily hold for all other soils. In particular, 
the coral found in the Pacific area is very different; it is very likely 
decomposed under the intense heat with CaCO^ giving off CO 2 to become 
CaO which may then be converted to hydroxide by moisture in the air. 

These processes should still lead to a distribution relatively indepen¬ 
dent of weapon energy but the curve may be vastly different from the 
Nevada results. Also shown in Figure 3 is the curve which Rand found, 
necessary to duplicate the results of Castle shot Bravo. Our calcula¬ 
tions lead to a similar result. 

The distribution of total concentration in the cloud must also be 
known. (Figure 3 has only specified the breakdown of this total into 
particles of various sizes.) For reasons discussed under Constraint 
E, we have assumed the total activity density to be constant within 
the cloud at any one altitude, but the vertical distribution has not yet 
been specified. 

The first assumption that generally springs to mind 
tion is asked is that the activity is probalfly uniformly distributed 
in the vertical direction as well. Consider, however, the fact that 
a cloud from a large weapon extends from altitude 60,000 to 120,000 
ft. Over this range the air density and pressure change by a factor of 
8. Now the pressure in the cloud must equal the external air pressure 
if motion is to cease and the temperature cannot long stay very different 
from that of the surrounding air. If this is so and if the cloud main¬ 
tains a stable shape, then the (mass) density must (by the universal 
gas law) also be near that of the surrounding air. This also balances 
the bouyancy forces due to gravity. 
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Thus, the mass of the whole cloud must he distributed in the same 
way as that in the surrounding atmosphere. (As already mentioned, 
this checks the expected proportionality of total cloud mass to yield, 
given the data of Figure 1.) If the further assumption of violent 
mixing during the short (10 min) cloud rise is invoked, the cloud debris 
must also be distributed in the same way as the density of the atmos¬ 
phere. 

In lieu of detailed data, then, it would appear that the most 
reasonable assumption is that total activity density is proportional 
to air density. There is an indication from early drone flights at 
Greenhouse that ion chamber readings were independent of altitude 
within the cloud. This is precisely what would be expected on the above 
assumption since the gamma ray path lengths would vary with air density 
in such a way as to exactly cancel the variation of activity corcentra- 
tdonj in other words, an ion chamber measures activity per unit air- 
mass which is what we have assumed to be constant. Further data on this 
point would, however, be most desirable. 

. In order to compute the fall times of the radioactive particles 
from various altitudes, they were treated as irregular spheres subjected 
to the standard aerodynamic fall equation. For very small velocities 
this equation reduces to Stokes Law and for very high velocities to 
the law for the turbulent region. At very low air densities, when 
the mean free path of the air molecules becomes comp arable to the parti¬ 
cle diameter, the fall velocity must be increased according to a slip¬ 
page equation. The fall velocities presented as the time in hours to 
fall 10,000 ft as a funtion of particle size and altitude are given 
in Table I. The computational details are presented elsewhere. 

Table II has been prepared from Table I. It demonstrates that the 
relative time for irregular spheres of various sizes to fall to ground 
is markedly invarient to the originating altitude. Thus Constraint 
C is satisfied and it is this fact properly exploited which results 
in the tremendous simplification of the Technical Operations computa¬ 
tional method. 

Finally, we present the conversion factors used to obtain the 
fission product dose-rates expressed in roentgen units. These are shown 
in Table HI. In comparing the fallout contour predictions of the 
various groups, due consideration must be given to the numerical 
values used by these groups which correspond to those presented in 
Table HI. For example, Rand at rhis writing uses 1200 in column (A) 
and unity in column (C). 

DOE ARCHIVER 

All of the constraints, assumptions and data required to predict 
fallout patterns have been enumerated above. At this point we will in- 
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Relative Particle Fall Time to Ground from Altitude, 
Normalized to 100 Units for 70 Particles 
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dicate how this information is employed, through reference to Figure U. 

In Figure U the cloud sizes containing the radioactive debris at the 
instant of stabilization for the three sizes of weapons are drawn to 
scale with the top and bottom altitudes and top-radius given by Figure 
la This shape is then divided into horizontal slabs having equal activ¬ 
ity per unit area. (Ten slabs suffice for 15 and 50 MT weapons, and 
five slabs for the 1 MT.) Due to the assumption that activity density 
varies as air density, the pressure drops across all slabs are the same. 
This results in the slabs having the radii and mean altitudes depicted 


TABLE in 

Roentgen Contour Constants 
(Fission Products Only) 




Flat Ground 
1 hr-R/hr per 

sq mi/fission Radioactive 
products on Decay Law 
Ground per KT (ii) 

of 

Thermonuclear 

Explosion 


(AxBxC) 



(ii) This t" 1 *^ law is a good approximation for moderata t but is'^nown 
to be inaccurate for times beyond 200 days; however, until better 
data are available this law seems to be adequate for purposes at hand. 

in the lower portion of Figure L. In the "zero wind case" all SSPfh^CHlVES* 
activity from each of the slabs must ultimately fall upon the ground 
with zero lateral displacement. The conversion factors listed in Table 
III and shown in Figure U result in the 1 hr-R/hr dose rates for various 




DOSE RATE PATTERNS FROM THERMO-NUCLEAR WEAPONS 

ZERO WIND CASE 



Figure 




















distances from ground zero. Thus the highest possible effective 1 hr-R/hr 

wina conditions is 
__ 1 MT weapon. The 

additional contamination near ground zero due to the small amount of 
radioactive debris in the stem is not treated in this model. Although 
the two-day cumulative dose for the zero wind case is not shown, it is 
readily calculable since: (a) the fall.times to ground of the various 
size particles from each horizontal cloud slice is predetermined by 
Table I; (b) the portion of activity associated with the various size 
particles is given in Figure 3j and (c) radioactive decay is assumed to 
obey the t“l*2 law. 

The effect of wind on zero wind patterns is to change the ground co¬ 
ordinates but not the fall times of the particles contained in each 
horizontal' "cl oucf slice. This distorts the shape and alters the area 
of the roentgen contours. The degree of distortion and alteration 
depends upon the meteorological conditions and is calculable from as- 
sumptions and constraints discussed above by a method similar in prin¬ 
ciple to that described by Rand in this Symposium Report. The results 
of these calculations are the eighteen 1 hr-R/hr and two-day roentgen 
cumulative dose contours compiled in another section of the Symposium 
Report. To facilitate comparison, our contours have been renormalized 
to a terrain factor of unity. We are of the opinion that a terrain fac¬ 
tor of 0.55 for fiat land is far more realistic. Thus with respect to 
what would actually be measured with a roentgen dose-rate meter 3 ft 
above ground, the roentgen levels of our contours as depicted (as well 
as those shown for the other groups) should be reduced twofold. Table IV 
tabulat-es the contaminated ground areas enclosed by fission product 
roentgen level contours for 1, 15 and 50 thermonuclear surface bursts 
for all six AFSWP winds. 

From the data contained in Table IV, we can draw the following con¬ 
clusions, some of which in view of the complexity of fallout phenomena 
could not have been anticipated: 

These contour patterns were conputed for a variety of rather 
widely varying wind conditions. Therefore, for a given burst energy, 
the contaminated area enclosed by any 1 hr.-R/hr contour is not strongly 
dependent upon the meteorology; this is also true for each two-day 
cumulative dose contour. 

For each wind condition, except for some of the 1500 R con¬ 
tours resulting from the 15 and 50 MT surface bursts, the areas enclosed 
by the same numerical value of the 1 hr-R/hr and two-day dose do not 
differ greatly. 
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dose rates- which can result from this model under any 
5?00 for the l5 MT and 50 MT weapons and R200 for the 


Areas of Contaminated Ground inclosed by the Roentgen Levels Noted for the Six AFSWP Wind Conditions 

(Areas Tabulated in 103 sq. mi.) 
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A Castle-Bravo type weapon (15 MT) detonated on land areas sim¬ 
ilar to that of the Eastern Seaboard would result in 1 hr-R/hr dose rate 
and two-day cumulative roentgen dose contours of 100, 500 and 1500 
whose contamination areas would encompass about 10,000, 3,000 and 600 
square miles, respectively. 

For any given wind condition the area encompassed by any dose 
rate or two-day dose, respectively, appears to be nearly proportional 
to the weapon size. 

Note that Conditions A, B and C are all moderate wind cases 
with varying degrees of shears Conditions D and E are of low wind veloci¬ 
ty; and Condition F is a high wind case. Thus, if meteorological con¬ 
ditions are classed solely according to wind velocity, namely as low, 
moderate and high, the phrase, effectively independent of , can be sub¬ 
stituted for the phrase, not strongly dependent upon , 'contained in 
the first Conclusion. 

The last two Conclusions enable Table IV to be reduced to Table V* 


TABLE V 

The Effect of Low, Moderate and High Wind Velocities on the Amount of Area 

Enclosed by the Roentgen Levels Noted 
(Tabulated values are average ground area in 10^ sq, mi/KI) 


Fission Product 
Roentgen Level 


Wind Velocities 



100 

0.6 

0.35 

0.8 

0.7 

0.7 

1.0 

5oo 

0.2 

0.2 

0.2 

0.3 

0.1 

0.3 

l5oo 

0.08 

0.1 

0.025 

0.1 

0#0 

0.06 


R is 1 hr-R/hr dose rate; D is two day cumulative roentgen dose* 
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IDEALIZED CONTOUR SHAPES 
APPLICABLE TO BOTH 2 DAY DOSE AND I hr. ~ 
100 TO 1500 ROENTGEN LEVELS 


CURVE A 





Table V should be a fairly good guide to the roentgen level flat 
ground area coverage to be expected from the surface detonation of all 
thermonuclear weapons under most meteorological conditions, provided 
that the land area is not comprised of coral. 

The fact that the contaminated contour areas seem to scale almost 
directly as the thermonuclear burst energy means that any known contamina¬ 
tion pattern for a given energy can be transformed at once into a pattern 
for another energy if the meteorological condition remains unchanged. 

Draw vectors from ground zero to various points along the periphery of 
the contour in question, and scale the length of these vectors by the 
square root of the ratio of the weapon energies. This working rule is 
subject to direct verification by measurement of our contour plots for 
each of the six AFSWP winds. This has been done with excellent agreement 
between the 15 and 50 MT weapons and fair agreement between the 1 MT and 
the 15 or 50 M? weapons. Further inspection shows that all of the dose 
and dose rate contours (somewhat idealized) are lemniscatic in shape, as 
depicted in Figure 5* Figure 5 shows two curves. Curve A is typical of 
contours when the wind 6hear is low or moderate (Conditions A to E). 

Curve B typifies the contour shape for high wind shear (Condition F). 

For contour curves both A and B, the terminus of maximum cross-wind 
dimension, d, is the end of a polar vector which makes an angle of about 
60° with the average downwind trajectory line. The relative ratios and 
scaled lengths of polar vectors a, b and c depend,of course, upon the 
wind velocity, but appear to be relatively independent of wind shear. 

The numerical values of these three polar vectors for low, moderate and 
high winds (see Conclusion (e) ) are tabulated in Table VI. 

In view of the data on the contour areas (Table V) , the R and D 
polar vector numerical values for each case should be nearly identical, 
except for the 1500 R contours under high winds. This is borne out by 
Table VI. Table VI indicates that a CASTLE-Bravo type weapon (15 MT) 
detonated on land areas along the Eastern Seaboard under a moderate wind 
condition would result in both a 1 hr-R/hr dose rate and a two-day cumu¬ 
lative dose 500 R contour which would extend about 100 miles downwind, 

10 miles upwind, and have a maximum cross-wind extent of about 40 miles; 
the total area encompassed would be about 3000 square miles. (See Table V). 
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TABLE VI 


The Effect of Low, Moderate and High Wind Velocities 
on Figure 5 Polar Vector Distances for the Roentgen Levels Noted 
(Tabulated Values are Polar Vectors in rai/(MT )2 


Wind Velocities 


Fission Product Polar Low Moderate High 


Roentgen Level 

Vector 

R 

D 

R 

D 

R 

D 

100 

a 

27 

15 

80 

58 

69 

70 


b 

17 

13 

8 

9 

5 

8 


c 

7 

6 

k 

6 

2 

li 

500 

a 

11 

10 

25 

25 

21 

36 


b 

8 

9 

5 

7 

2 

h 


c 

5 

12 

l 

h 

-7 

-2 

1500 

a 

8 

8 

9 

13 

0 

17 


b 

6 

6 

2 

5 

0 

1 


c 

k 

U 

-1 

2 

0 

-3 


R is 1 hr-R/hr dose rate. D is 2-day cumulative roentgen dose. 

(-) sign for polar vector c denotes that c is above ground zero 
in Fig. 5. 

a is maximum downwind distance to contour. ®Ql 

b is 60 ° polar vector* 1.73b is maximum crosswind distance, 
c is maximum upwind distance to contour« 
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Figure 5 and Tables V and VT can be used as a rapid method to esti¬ 
mate 2-day dose and 1-hour dose rate contour patterns over contaminated 
land areas as follows: 

Use Figure 1 to determine altitudes of bottom and top of 
initial cloud resulting from a surface thermonuclear burst of energy in 
MT. 

To characterize wind shear, qualitatively examine wind vector 
data for only those altitudes which are bounded by cloud bottom and 
top. If wind shear is high, contour patterns will be similar to Curve 
B, Figure 5; otherwise it will be like Curve A. Since for altitudes 
above 1|0,000 ft, high wind shear is rare, most contour patterns will 
be like Curve A. 

Use an air density table to determine a mean altitude, H, 
which corresponds to the center of mass of an arbitrary cylinder 
bounded by the cloud top and bottom altitudes. If the burst energy 
is 1 MT, H is 15,000 ft; if the energy exceeds 10 MT, H is 70,000 ft. 

Consider a 10 Qa particle that was originally located H feet 
directly above zero ground. Use fall times given in Table I together 
with the wind data to ascertain the location of this 100^ particle when 
it reaches the ground. 

Draw a line through ground zero and the location of this part¬ 
icle. This is the average wind trajectory axis. Measure distance in 
miles from ground zero to this particle location. If this distance is 
less than 100 miles, classify as low wind; if distance is between l5o 
and 750 miles, call moderate; and if distance exceeds 750 miles, it is 
a high wind case. 

Now draw polar vector distances for the 100, 500 and 1500 
roentgen levels as given in Table VI for appropriate wind velocity case 
and burst energy. Sketch in contour lines as depicted in Figure 5 for 
appropriate wind shear case. Keep in mind that enclosed areas must 
satisfy Table V. 

D °E ARCHIVES 

Roentgen levels of these contour lines are valid for the con- ' 
version factors listed in Table III. If you are using other values for 
these factors, alter roentgen levels accordingly. 



Fission Product Biological Damage Dose and Contribution of U-Np 
-, _ « Activity Thereto 

- Introduction 

Biological damage from exposure to radioactivity is complicated by 
the tendency of the body to repair itself. Thus the maximum effect from 
a given total dose depends upon the length of time over which the dose was 
received. If the duration of the dose is less than two or three days, 
the effect of self-repair on the danger level is slight, but times 
longer than this are of current interest. 

Data on this subject are still rather sparse. However two studies 
have been made which give two estimates of the effect* The results 
differ principally at periods of the order of six months or more, where 
data do not yet exist in sufficient quantity to support either theory. 

The mathematical forms which summarise the two studies differ considera¬ 
bly, and both will be considered in what follows. 

Method I 

A British study states that exposure to a constant dose rate d 
per day for t days is equivalent, at the end,of that period, to a one- 
day dose ("Damage Dose") given by 

D = d t°* 6i4 (t > 1 day) (la) 

Manif estly the restriction "t >1" day is required since for very short 
times, the effective dose must be more nearly 

D *= d t (t < 1 day) (lb) 

Equations (la) and (lb) fit data taken with a constant dose rate 
d per day. To get relations for dose rates which vary in time due to 
radioactive decay, consider the effect of extending t by a small amount 
At which is added at the beginning of the interval t. The additional 
damage dose due to this increment is easily seen to be: 

AD = . if t < 1 day DOE ARCHIVES 

*■0.64 t" 0-36 dAt if t > 1 day 

This relation says that the effect of a pulse (dAt) of radiation 
is (a) unreduced for times, less than 1 day, and (b) reduced by a factor 
(0.6k t-0.3o) after times, t, longer than one day. If then, exposure 



occurs from time t. to t^ at variable fission product dose rate of 
c/t^* 2 , the final damage nose will be (provided 1 2 > t^ + 1 day): 


/ t + / 


If t 2 < ti 4- 1 day: 


V 1 


t 2 -i 


0.64 cdt 

<Vt/>-5V- 8 


( 2 ) 


-/ 


cdt 

TT^ 


(?) 


Therefore if t^ < t^ + 1 day 


5 c(t. 


- 0.2 


t; 0 - 2 ) 


(4) 


In equations (2) and (4) time is in days and the dose rate from 
decaying fission products is taken as c/t 1 * . The constant c is there¬ 
fore in roentgens per day at one day. To convert to more usual units 
let R be the dose rate due to flsBion products in roentgens per hour 
at oni hour. It is then easily seen that 


c = R 0 /(24) 0 ’ 2 - (0-530)R o 

In equations (2) and (4), time is still measured in days. 


(5) 


The second integral of equation (2) has been computed numerically. 
The results along with those from equation (4) have been plotted in 
Figures 6& and 6b which give directly the ratio , for fissicei 

products alone and various entry and exit times, as implied by the 
British results, equations (la) and (lb). An example will be given after 
discussing other approaches and refinements. 


Method II 
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An American investigation of the same questions has yielded a 
different expression. The concepts used here are those of a "permanent 



BIOLOGICAL EFFECTIVE DOSE FROM FISSION PRODUCT DECAY (BRITISH DATA) 
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Figure 6a 





























































BIOLOGICAL EFFECTIVE DOSE FROM FISSION PRODUCT DECAY (BRITISH DATA) 



p Of e Oomnge Dose From Fission Products only in Roentgens 
R = Fission Product Dose Rote. R/hr ot I hr. 


t| s Entry or Arrival Time in Days 


1 2 weeks 2 Irfla 4 s 2feo? 8 910 4mo. 1 

to 100 


IMMiESMi 


1 yr. 4 b fc 2'yr? 9 3yr. 5yr. 2 
1000 

Abscissa: Esi t Time in t g (days) 


10 yf. 5 6 7 8 9 1 

10000 


Ordinate: R 


DOE ARCHIVES 


Figure 6b 
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injury* 1 and a "recoverable injury 11 . The former is empirically 20% 
of the "instantaneous injury" while the latter decays exponentially 
with an apparent recovery half-life of about 2.U days. In symbols, 
if I is the permanent injury dose, I r the reparable injury dose, and 
R the dose rate then 


■ 0.2 R and r “ 0.8 R - k I 


where k 


Q*fi9_3 B 0.2? per day. 

2.4 


The damage dose at time, t, is then 

D(t) = I_(t) + I (t) (7 

F r 

For the case of fission product decay R ■ c/t^ * 2 and exposure from 
t, to t OJ equations (6) and (7) may be integrated. The result i3 


n r 

<°- 2 > °] 7I7Z ♦ <°- 8 > c J ^-n 


k(t - tj 


dt (8) 


Here time may also be measured in days and equation (5) again 
applies: 

c = (0.530)R o (5) 

where R 0 is the dose rate in H/hr at one hour due to fission products. 

Equation (8) has also been computed numerically. The results are 
given in Figures 7a and ?b which are completely analogous to Figures 
6a and 6b c 


)ose Without 


mission Product Cumulative Dose, 


To emphasize the effect of biological recovery, there is displayed 
in Figures 8a and. 8b tne integral dose, which would correspond to the 
case of no self repair. The equation for the integral or cumulative 
dose is simply equation (k) without the restriction that t 2 < t +1 day. 
It will be noted that such a computation is in error in two respects. 
First the time of maximum danger would always appear to be the exit time. 
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Secondly, the maximum dose may be much too high. These points are il¬ 
lustrated in Figure 9 which shows the integral and damage doses versus 
time for an entry or fallout arrival time of 10 hours. For long times, 
the integral dose is at least a factor of 3 above the actual biological 
damage dose. The spread between the "American" and "British" curves in 
Figure 9 gives some idea of the accuracy (or rather, paucity) of the 
original data as a function of exposure time. 

Neptunium Activity 

Before discussing further the application of results obtained above, 
it is necessary to consider an additional factor. The curves developed 
above concern only the gamma activity from fission products. There is 
in addition another activity which has a different decay law. This 
arises from capture of neutrons by 0-238 to yield U -239 which decays to 
Np-239. Both give off biologically significant radiation. U-239 
gammas have a mean energy of Q.07U ME? while that of Np-239 (including 
the effect of internal conversion) is 0.223 ME? per disintegration. The 
half lives are 23.3 min. and 2.36 days, respectively. 

Calculations lead to the following results. If R ? is the dose rate 
in R/hr at one hour for fission product activity, the integral dose 
from U and Np assuming unity capture to fission ratio is ' 


I NP = 0.313 [ *1 - e" 1 * 2 * 5 ^ 2 ] + e ~°' 29 11 - e “ 0 * 2 ? * 2 (?) 

R o 

The units for I„ are R/day and t-^ and are the entry (or arrival) 
and exit times in days. The biologically effective dose (damage dose) is 



( 10 ) 

0 . 232 e“ 0 o 29 t 2 (t 2 -t 1 ) 


where the units are the same and 1^ is given by DO£ ARCHIVES 

The damage dose (10) is to be added to that previously computed for 
fission product gamma rays. To give an indication of the relative impor¬ 
tance of the two terms (for unity capture to fission ratio). Table ?II 


p /* 0 1 


is given by (9). 


gives the relative dose rates from the two processes as a function of 
time for a unity capture toTfission ratio. Table VII shows that from 


1 to 10 days after shot time, U-Np-239 make a significant contribution 
to the dose rate. 
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INTEGRAL DOSE FROM FISSION PRODUCT DECAY 
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COMPARISON OF PREDICTED DOSES FROM FISSION PRODUCTS 
ONLY, FOR 10 HOUR ENTRY OR ARRIVAL TIME. 



EXIT TIME (DAYS) 



TABLE VII 


Dose Rate Ratio of Np-239 and U-239 to Fission Products 
at Various Times Following Detonation 
(Unity Capture to Fission Ratio) 

Hours Days 

Time 0.1 0.3 0.5 2~~ 5 10 2k 2 U 7 14 30 

Dose Rate 


Ratio Np 0.03 O.lU o.ll 0.06 0.08 0.17 O.U1 0.66 0.8^ 0.71 0.20 0.0U 
and U/F. P. 



In order to simplify the addition of this extra activity to that of 
the fission products, there is presented in Figure 10, at, the ratio 
Of the biological damage dose due to Np and U to that due to fission 
products alone (for unity capture to fission ratio): 

a “ —SE (Ratio of equation (10) to equation (8) (11) 

D f 

Figure 10 presents calculations for a for various entry and exit times, 
•where D„ is given by (10) and D f is the ("American'*) damage dose given 
in Figures 7a and 7b. 

To use the information in Figure 10 the damage dose due to 
fission products alone is first determined from Figures 7a or 7b for 
given values of t, and t ? . These same times are then used to obtain a 
value of a from Figure ID. The complete damage dose D is then found by 

D"D f (l + FG) (12) 

where F is the capture to fission ratio for the weapon in question. 

(F ® number of neutrons captured in U-238 for each fission which takes 
place at the instant of detonation.) 

To complete the presentation. Figure 11 shows the analogous ratio, 
P,for integral doses as computed from (?) and the standard integral- 
dose formula (Figures 8a and 8b). 
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Use of the Data in Fallout Calculations 

In calculations such as discussed in Section I, the final results 
are presented in the form of 1 hr-dose-rates and 2-day integral (cumulative) 
dose. The employment of the data presented in this section to modify 
dose rate patterns is straightforward but the effect on integral dose 
is more complicated. We will first consider the effect on dose rate. 

Dose Rate. - The 1-hr dose rate -contour cannot be used to 
predict actual dose rates on the ground until most of the radioactive 
debris contributing to the roentgen rate has reached the ground location 
of interest; depending upon the distance from ground aero, this may be 
as long as twenty hours following shot time. Thus the dose rate contours 
should not be used to predict dose rates as would be measured on the 
ground by a roentgen dose-rate meter for times less than twenty hours 
after explosion. The following example illustrates the use of these 
contours to predict the dose rates at a given location: 

Example. - 

Surface detonation of 15 MT weapon with capture to 
fission ratio of unity. 

Contour predicts 1500 1 hr-r/hr due to fission 
products, 70 miles NW of ground zero. 

What is dose rate as measured by roentgen dose meter 
2k hours after detonation? 

fission product dose rate * 3-500 = 3-500 = 33 R/hr 

From Table VII, (U-Np/F. P.)^ ^ = 0.1a 

U-Np dose rate = 0.U1 x 33 * lU R/hr 

Total roentgen dose rate at 2h hrs * Ifl R/hr 



Integral Dose. - It must first be emphasized that all graphs 
in Figures 6 through 11 are plotted on the assumption that the amount of 
contamination is constant. No further fallout particles are considered 
to arrive during the time of exposure. The modifications due to arrival 
during exposure will be discussed below; for simplicity we sh all first 
suppose that, as predicted in Figures 6 through 11, all contamination is 
already present at time t-^.. 
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With this proviso, consider first a simple example: Exposure 
occurs at a location where the fission product dose rate is 500 R/hr at 
1 hr. Exposure starts at t-j_ = 10 hr = 0.42 days and exit occurs at 
= 2 weeks a 14 days. 

From Figure 6b we find that the "British" curve for t-^ * 0.42 
has the value Dj/R = 0.47 when tp = 14 days. Thus the damage dose 
due to fission products only is 0.47 R Q * 235 roentgens at the exit 
time of 2 weeks. The American data predicts approximately the same 
biological effect; from Figure 7b the damage dose is 0.42 R 0 = 210 
roentgens at 2 weeks. 

Note, however, that both curves have an earlier maximum. From 
Figure 7a, for example, the damage dose (from fission products only) at 
2.6 days was 0.66 R = 330 roentgens. In effect, the mechanism of 
biological repair was unable to "keep up" with the damage at early times 
when the radioactivity level was high. After 2.6 days, however, the 
activity decayed to the point where repair mechanisms could somewhat more 
than make up for the damage. 

Consequently, the danger is measured not by the final damage 
do se but by its maximum value 'during the ‘ time of ex po sure . The (possibly 
delayed 1) manifestations of injury will in general correspond to this max¬ 
imum and not to the final value. 

So far, we have considered only the damage dose due to fission 
products alone. This must be further increased by the effect of U-Np 
activity. From Figure 10, for example, the t, = 10 hr curve shows that 
a - 0.56 at tp = 2 weeks. The total damage dose at the exit time is 
thus 210(1 * 0.56F) * 322 roentgens if the capture to fission ratio F of 
the weapon was unity. 

Recall, however, that it is the maximum dose which represents 
the actual danger. We must therefore maximize D^-(l + a) (if F = 1) 
where D^. is obtained from Figure 7a and a from Figure 10. Calculating 
this quantity for the 10 hr curve in the two figures, one soon finds 
that the maximum is very flat and occurs at about 3 days where 
(D|./R )(1 + gc) - 0.99. Thus the maximum danger in this exposure occurs 
at 3 days and the total damage dose is then 500roentgens. (Assuming 
T^~i) ‘ ‘ 

This result is to be contrasted with integral dose which may 
be determined from Figures 8b and 11. The former gives a fission product 
dose of 1.6 R Q for the same times while the latter leads to a correction 
factor of (1 * 0.55F). If F * 1, the total integral dose is therefore 
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12l|0 roentgens. Not only is this 2$0% too high in terms of actual 
biological damage but the maximum always occurs at the exit time, thereby 
leading to the false conclusion that earlier evacuation would be helpful. 
We saw above, that unless evacuation occurs before 3 days, there is no 
reduction of the danger. 

At this point, it may be noted that even moderate shielding can 
provide considerable protection. The maxima which meas\ire the danger 
of biological damage all occur earlier than about one week and their 
high values are due almost entirely to the very high dose rates occur¬ 
ring at times earlier than this. Protection during these early times 
is therefore very important and helpful. Thus retreat to a shelter for 
a period of one week under the circumstances already used in the above 
example changes all calculations to the curves for t-j = 1 week. Suppose 

the shielding reduces the earlier 3-day maximum by a mere factor of 10 
leading to a dose of $0 roentgens while in a shelter. The one week 
curves in Figures 7 and 10 then show that the quantity (D*-/E 0 ) (1 + Fa) 
has a maximum (assuming F = 1) of 0.16 at the exit time, 5 weeks. The 
maximum damage dose is then 80 roentgens due to exposure after leaving 
the shelter plus (0.ii2)(l + 0.56)(R^IO) = 33 roentgens lingering damage 
dose from that received while in shelter. The net result is thus a dose 
of 113 roentgens instead of $00. If the arrival time is smaller than 
10 hrs, the effect of even moderate shielding is even greater. 

Finally, we consider the result of lifting the restriction 
stated earlier on the arrival of contamination. All of the above dis¬ 
cussion assumed that fallout was complete before exposure started. 

We will consider the contribution to the 2-day biological ef¬ 
fective dose from radioactive debris that arrives at a ground location 
five hours after shot time. It will be comprised of two parts, that due 
to fission products alone and that due to U-Np. We will again assume 
a unity capture to fission ratio. If the fission-product-only 1 hr-R/hr 
dose rate of this debris is R q then from Figure 7a, . . 


the fission product contribution to the 

2-day biological effective dose ■ 0.95 R c 
while from Figure 11 the Np-U contribution 


is given by 0.3U x 0.95 R r 


0.32 R 


Hence the total 2-day biological effective dose 
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1.27 R„ 
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We will now compare this total value with the fission product only 
2-day integral dose from 5 hrs arrival time debris whose 1 hr-ft/hr dose 
rate is also R Q . We refer to Figure 8a and find this to be 1,3 R 0 , and 
note that this is identical to the value for the 2-day biological 
effective dose. If we repeat these calculations for any other arrival 
time between 0.5 and 15 hrs and again make the above comparison, we 
reach the same conclusion. Now it can be proved that the location of 
the perimeter of any two-day dose contour equal to cr greater than 100 
R is determined solely by radioactive particles which arrive at a ground 
location in 15 hours or less. 

In view of the above, the following conclusion can be stated, 
namely: 

Conclusion (A). - 

Irrespective of arrival times of the radioactive 
debris, for all thermonuclear surface detonations where the capture 
to fission ratio is about unity, the fission product two-day contour 
is the two-day damage (biological effective) dose contour, provided 
that the contour in question is at least 100 roentgens. 

We will now consider damage dose for periods greater than two days. 

Biological damage is apparently additive. The effects of self¬ 
repair require one to successively down-grade" earlier doses before 
adding them to obtain the total effect, but the addition otherwise 
appears to reproduce the data on damage. (This additivity is implicit 
in both the British and American damage dose relations, equations (la, 
lb, 6 and 7).) Further, two days after shot time, the 1 hr-R/hr dose 
rate is a real measure of the actual dose rate at any ground location 
within the two-day dose 100 R contour, since by this time all radio¬ 
active particles have arrived. 


The information in the preceding paragraph plus Conclusion (A) 
means that the biological effective (damage) dose at a non-shielded 
reasonably flat location at any time, t, is given by the following 
equation: 


For t > 2 days, D„, ^ 100 roentgens and capture to fission * lj 

then the damage dose at uhe time, t, is 

B t " B 2 d {°- 2 * °- 8 e -° <29(t - 2) } * R 0 < jp- ) ( 1 ♦ « ) (13) 

o 

where t is tin® in days. D , and R are the respective fission product 


®GEaj? 



C ’~ Tn/ Ss 




ft 
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two-day dose and 1 hr-R/hr dose rate levels in roentgens at the ground 
location in question* The right hand term is the "down-grading" of the 
2-day damage dose with time and comes from equation (6). The left hand 
term is the damage dose which results from an exposure starting at 2 
days and ending at t* It comes from equation (12) where and a 

at the time, t, can be found from the curves labeled two-day entry time 
on Figures (7) and (10), respectively. 

Equation (13) can be rearranged to read: 

D t " D 2d M(t) + R o Q(t) < ll 0 

where the functions M(t) and Q(t) are cbvious by comparison with 
equation (13) • For convenience they are tabulated in Table VIII. 


TABLE VIII 


Values of M(t) and Q(t) of Equation (ll:) 



M(t) 1*00 0.80 0.65 0.5b 0.39 0.28 0.20 


Q(t) 0.00 0.28 0.1*0 0.50 0.5l 0.1:5 0.35 

These values are based upon continuous exposure to radioactive 
fallout for time periods of at least two days duration where time 
period is measured from shot time. See equation (13) for other 
restrictions. 


From Table VIII we see that the time after shot at which maximum 
damage dose due to continuous exposure to radioactive fallout can 
never exceed 7 days. The time at which the maximum damage dose takes 
place and its numerical value depend only upon the 2-day dose and the 
1 hr-R/hr fission product roentgen contour values of the ground lo¬ 
cation in question. To illustrate the use of equation (lb), we will 
now give two numerical exmplesj Example (i) leads to a practical 
working rule, while Example (ii) illustrates the use of equation (lb) 
to solve a problem which involves protective shielding. 


LM2 


DO E 






Example (i).- 


Surface detonation of 50 MT weapon with capture to fission 
ratio of unity* 

Contours predict 300 R 2-day dose and 350 1 hr-R/hr due to 
fission products 270 miles NW of ground zero# 

What is maximum biological effective (damage) dose due to 
continuous exposure and when after shot time will it occur? 

From equation (lU) and Table VIII we obtain 


Days 



2 300 x 1.00 + 350 x 0.00 - 300 R 


3 300 x 0.80 + 350 x 0.28 * 31*0 R 

1* 300 X 0.65 + 350 X o.Uo - 335 R 

5 300 x 0.5U * 350 x o.5o » 335 R 

7 300 x 0.39 + 35o x 0.51 - 295 R 

Thus maximum damage dose is 3lt0 R and will occur 3 days 
after shot time. 

In Example (i), where ratio of R /D„ . is 1.2, the maximum damage 
dose is only 10$ greater than D«.. Examination of all of the contour 
patterns for roentgen values of interest for all of the AFSWP winds 
shows that R AUj is usually of the order of unity, and essentially 
never exceed! 1.2. This fact results in the following working rules 

Conclusion B,- 


For a surface burst thermonuclear weapon with a 
capture to fission ratio of about unity, the maximum biological 
damage dose that can ever be received at any geographic location 
(with no shielding) is essentially the 2-day fission product integral 


dose and hence can be read directly from a 2-day fission product dose 
prediction contour map. 


We new proceed to Example (ii).- 
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Example (ii).- 


Surface detonation of 15 MT weapon with capture to fission 
ratio of unity* 

Contour predicts 1500 1 hr-It/hr and 2000 R 2-day dose due 
to fission products, 70 miles NW of ground aero. 

A person is in shelter of shielding factor 10, for the first 
two days after shot. After first two days he leaves shelter and is ex¬ 
posed in open flat area for 8 hours each dayj remaining l6 hours spent 
in shelter. 

What is maximum biological dose received by the person and 
when will it occur? 


Effective D 0 - 2000 

2d " Yq " 


200 R 


Effective R_ - 1500 x 8 + 1500 x l6 

O ’ r '’AfT~ ¥: m 1T ’ ■WTjMyj-. "nr - 


600 1 hr-R/hr 


From equation (lb) and Table VIII we obtain 


Days 

after Shot 


2 200 x 1.00 + 600 x 0,00 E 200 R 

b 200 x 0.65 + 600 x O.bO 3 370 R 

5 200 x 0.5b + 600 x 0.50 ■ UlO R 

7 200 x 0.39 «■ 600 x 0.5l - 380 R 

Thus the person will receive a maximum damage dose of blO 
roentgens and this will occur 5 days after shot time. 
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QUESTIONS AND ANSWERS FOLLOWI1C DR. HENRIQUES 1 PRESENTATION 


Question:* Dr, Kellogg, Rand Corp. 


1 


This last slide Dr. Henriques showed brought up a question. 
Would you include the attenuation factor far the terrain or 
are they based on the infinite flat plain? 


Answer: Dr. Henriques, TOI 

Mine include both the 
terrain factors. 

Question: Dr. Genevese, AFOIN 



DELETED 


Answer: Dr. Henriques, TOI 


When you have a nuclear explosion, two main things happen. 
One is you get fission products that fall to the ground; 
and you get some U-239 which goes into Neptunium-239 and 
these are normally not considered with the fission products. 
These curves provide a means to do it if you know what 
fraction of the bomb goes to fission products and vhat goes 
to U-239. It makes no contribution to the fission products 
dose rate early but it does in time, because of the dif¬ 
ference between the 1.2 law and the exponential decay. 


Question: (Unknown) 


Could you amplify your remarks about the degrading effects 
of terrain? How much does it vary with different types of 
terrain? What is the cause of it and so forth? 


Answer: Dr. Henriques, TOI 


There are magic numbers you can use and they are all avail¬ 
able in reports and all roughly the same. I just don't know 
how to answer it. It depends on whether you are in the city 
and it comes down on buildings or if you are on the ground. 
The reason for using .55 instead of .7 is not basically 
important. Any open terrain or open field will not be as 
flat as the Salt Lake flats. 
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Question: 


Answer: 


Comment: 


Question: 


Answer: 


Comment: 


Comment: 


Dr. Krey, CRL 

You think the different types of soil around ground zero 
will have an effect? 

Dr, Henriques, TOI 

Yes. Two different things, coral and Nevada dirt, there is 
a real basic difference. If you use the NRDI, curve, it is 
impossible to duplicate Bravo with the NRDL particle. 

Dr, Kellogg, Rand Corp. 

There was one point you mentioned - the geometry of the 
cloud and that altitude is an important point. The height 
of the cloud as being constant through a rather large range 
of yields - this seems not to be borne out when you plot the 
CASTLE results with IVY results. The entrainment at each 
altitude is relatively insensitive to yield. When you put 
in more energy, the cloud would become higher. The dif¬ 
ference might not be very large, but I think it would be 
incorrect to assume that the altitude is not a function of 
yield. 

Dr. Scoville, AFSWP 

Isn’t the bottom of the cloud likely to be at the bottom of 
the tropopause? 


Dr, Kellogg, Rand Corp. 
I don’t know. 


Dr. Pelt, LASL 


The Chemical or physical nature of the soil are more impor¬ 
tant to the fall-out. Dr. Henriques mentioned coral and I 
think I was talking mainly about soil. Some have proposed 
that the physical nature is more important. 


Dr. Rapp, Rand Corp, 
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I would like to clarify our position. Our way of looking at 
it is the earth that will eventually become active must get 
into the fireball almost instantaneously. 



Commenti Dr. Krey, CRL 

At CASTLE, we found quite a lot of particles over 100 microns 
to have activity principally on their surface. It was hard 
to see how they could get into the fireball and have activity 
strongly on their surface. These were actually found at 
CASTLE. 

Comment: Dr. Rapp, Hand Corp. 

I will quote some work that was done at Rand. They investi¬ 
gated this condensation principle. According to this model, 
the silicates come out and then some of the fission products 
are condensed at a later time, so this model definitely 
accounts for not having the fission products evenly distri¬ 
buted throughout the fall-out particle. Most of it is from 
the condensation particles. 

Comment: LTGOL Lulejian, AEDC 

I think I remember reading in a GREENHOUSE report where they 
found a small piece of chlorophyll which had fission products 
on it. Apparently everything is happening. Condensation, 
vaporization, and other things. Hence, a leaf coming down 
with a lot of activity on it. 
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PREDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS FUNCTIONS 
OF YIELD AND METEOROLOGICAL CONDITIONS; RECENT DEVELOPMENTS 
IN WEATHER FORECASTING TECHNIQUES FOR THE PACIFIC PROVING GROUND 

PART I 

LTCOL C. t). Bonnot, USAF 
Joint Task Force SEVEN 


Introduction 

We have been talking the last few days about fall-out, and now we 
will talk about weather. We know we can only have as good fall-out 
analysis as we have weather forecasts. We have discussed that pretty 
thoroughly. From the various fall-out models given, I thought it would 
be well to go back and discuss the Bravo weather situation. I want to 
give credit for this work to the Oahu Research Center in Hawaii. The 
work was not done by me personally, but by Ohmstead and Dean under 
Dr. C. E. Palmer's direction. The report is unclassified and it was 
deliberately made that way; there is no reference to Bravo or to shot 
times. They selected this date because there is a lot of good weather 
data, so the report states. The Oahu Research Center is part of the 
Institute of Geophysics, UCLA, and operates under contract to the Air 
Force Cambridge Research Center, They have both military and civilian 
personnel assigned and the Task Force works with them. They assist dur¬ 
ing the test period and then we give them all the raw weather data. 

They now have data from all the Pacific tests. They are a storehouse 
of all the raw data and weather information for the tropical Pacific. 

It might be well for everyone interested in Bravo to get this report. 

It has now been published in its final form as Special Report No. 1, 

Oahu Research Center. The report is unclassified and available to any¬ 
one who wants it. Just ask AFCRC at Boston-, or contact me and I'll see 
that you get a copy. 

This presentation will be in three sections. I am going to talk 
about just the horizontal field, or streamline analysis, for Bravo day. 
Commander Rex will speak of the vertical motions for Bravo day, which 
are also analyzed in the report mentioned above. He also will discuss 
five other situations which he has analyzed for vertical motions, and 
then I will close with a short discussion of the proposed rawinsonde 
observing net for the Southwest Pacific which we feel is necessary if 
adequate analyses and forecasts of the weather and wind fields are to 
be furnished to the Commander, Joint Task Force SEVEN. 


Horizontal Fields 

Hie synoptic situation changed very little during the three day 
period between February 28 and March 2, 1954. At any given level, the 
wind directions and speeds did change slowly but these changes were due 
to systematic internal developments within the cyclones and anticyclones 
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in situ. No. synoptic systems moved through the area; no new synoptic 
systems developed in the area during this three day period. 

In the vertical, two transition levels (shear levels) happened 
to intersect various parts of the analysis levels. The transition 
levels occurred near 10,000 and 60,000 feet. The apparent synoptic 
changes at these analysis levels are due to small height changes in 
the shear levels together with some slow internal synoptic intensi¬ 
fications. 

During the three day period, at 10,000 feet, there was a gradual 
intensification of two cyclonic vortices that lay in a trough oriented 
east-west between Wake and Rongerik. (Figures 1, 2, 3* Note: These 
maps have been analyzed by the isogon-isotach method on a 1 : 7 , 500,000 
scale. All time changes of the winds in the horizontal as well as in 
the vertical have been incorporated into the isogon-isotach analysis.) 

Only one cyclonic vortex appears on the analyses. It is located north 
of Rongerik; the other lies in the area covered by the legend box. 

These vortices slowly intensify during the period, but do not descend 
to the surface. They are reasonably well fixed by reconnaissance 
reports (10,000 feet only) on February 28, March 1, and later again on 
March 3. The gradual intensification of the trough in which' they are 
imbedded is apparent by the westward movement of the neutral point from 
Eniwetok on February 28 to the Guam area by March 2. In other words, 
the higher level westerlies have descended into the analysis level. 

An anticyclone at 10,000 feet is located near Kwajalein on 
February 28 (Figure 1), and during the period it moves slowly south- 
westward to Kusaie. At the same time another anticyclone appears to 
the north of Ponape (based on off-time reconnaissance reports). It 
also moves slowly southwestward towards Truk during the three days. 

Both systems descend further into the 10,000 foot surface during the 
period. The combined effect of the slow intensification of the cyclones 
and the equatorial anticyclones at 10,000 feet over the Marshall Islands 
is to maintain northerly wind components over the Rongerik area below 
20,000 feet. 

The subtropical anticyclones are located north of 20°N at 10,000 
feet (just off the map). (Figures 2 and 3.) They weaken vertically 
north of Wake, completely disappearing as wind singularities above 
10,000 feet, and cannot be found in the speed field above 20,000 feet. 

The systems have very little effect on the weather of the Marshal-1 
Islands. DOE ARCHIVES! 

At 30,000 and 40,000 feet the How pattern around the equatorial 
anticyclones is the same with height and with time, giving southwesterly 

























components, over the Bikini area. (Figures 4, 5,6,7) A gradual 
southwestward shift of this belt occurs with time, in agreement with 
the southwestward displacement of the anticyclones at 10,000 feet. 

At 50,000 feet southwesterly flow appears on February 28. (Figure 
8) On March 1 an anticyclone descends into the analysis level at map 
tiros, but ascends again after a 6 hour period, and is not evident at 
50,000 feet on 2 March. (Figures 9 and 10) 

The horizontal picture of the wind flow at 60,000 and 70,000 feet 
is complicated by the presence of the tropopause and the low level 
stratospheric easterlies. (Figures 11 and 12) Any horizontal analysis 
in these layers represents only the height changes in all of these 
shear levels, and therefore has not been performed except for March 1. 
(Figures 13, 14, 15) Wind speeds at these levels are light. 

Steady and strong easterly winds persist at 80,000 feet through¬ 
out the period. 

In conclusion, analyses for all levels from the surface through 
80,000 feet for 1 March 1954 are shown (Figures 16, 17, 2, 18, 5, 6, 

9 , 11 , 12 , 14 ). 
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Figure II 



















































PART II 

VERTICAL MOTION IK THE MARSHALL ISLANDS AREA 

CDR Daniel F. Rex, USN 
Joint Task Force SEVEN 


In the Task Force we are interested in the field of vertical 
motion for two reasons] first of all, through an analysis of vertical 
motions we my be able to gain a better understanding of the circula¬ 
tion in the tropical atmosphere and secondly, if the vertical components 
are sizeable they will certainly affect the fall-out pattern produced 
following an atomic explosion, G. A. Dean and V/, D, Ohmstead, in Oahu 
Research Center Special Report No. 1 under Air Force Contract No. AF 
19(604)-546, have computed the vertical motion at various levels for 
'1 March 1954 at 0300Z from streamline analyses of the horizontal wind 
field. In order to extend their interesting results it was decided to 
compute the same and other cases using a different technique. So my 
purpose is first, to describe the vertical motions "typically" observed 
in the Marshall Islands area and second, to demonstrate the computa¬ 
tional method which was used—a method which is ^perhaps simple and 
objective enough to be applied successfully under operational conditions. 

We have available six cases in which the vertical motion in the 
area of interest has been computed] 0300Z 1 March 1954 by Dean and 
Ohmstead; 2100Z 23 February 1954, 0300Z and 1500Z 1 March 1954, and. 

0300Z 2 March 1954 by Rexj and one case during GREENHOUSE computed by 
Dr, C. E. Palmer and published in the GREENHOUSE meteorological report. 

It is obvious that no general conclusions may be drawn from such a 
meager statistical sample. Nevertheless, one can perhaps obtain a 
gross, qualitative picture of vertical fields common in the Marshall 
area. Figures 1 though 6 illustrate the fields of vertical motion 
obtained at the 2000, 1000G, 20000, 30000, 40000 and 50000 foot levels 
by Dean and Ohmstead and are generally typical of the computational 
results for all cases both with respect to the magnitude of the values 
and to the geometry of the patterns obtained. A comparison between 
the vertical wind speeds drawn in Figures 1 through 6 with those 
obtained in each of the other cases computed, gives a maximum varia¬ 
bility of ± 25/6 of the values illustrated for ary particular level. It 
may be seen from the figures that the vertical motions are organized 
into upward and downward moving currents or cells having horizontal 
dimensions of many hundreds of miles. It is also apparent that the 
magnitude of the vertical component increases with increasing altitude, 
and there is a tendency toward larger cell dimensions at the higher 
levels. In a very brief way this summarizes the computational results. 

For comparison purposes it should perhaps be noted that all vertical 
speeds in Figures 1 through 6 are given in centimeters per secondj one 
centimeter per second is equal to one hundred and twenty feet per hour 
or thirty centimeters per second is equal to 3,600 feet per hour. 

DOE ARCHIVES 
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Figure 2 















































Because the fields of vertical motion are obtained indirectly 
through computation of the horizontal divergence it is important to 
investigate the reality of the computed patterns. For this purpose 
two checks were used: the distribution of observed weather was com¬ 
pared with the computed motion; and vertical motions, computed for the 
same tine by two different methods, were compared. In Figure 7 ia 
shown the distribution of weather in the Marshall Islands area at 
2300Z on 28 February 1954; the observed cloud amounts and types and the 
location of precipitation elements are indicated by the following 
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Blank areas are cloud and precipitation free DOE ARCHIVER 


Figure 7 was taken without modification from the original weather 
briefing chart used just prior to the Bravo event; it is 3imply a repro¬ 
duced section of the original and was constructed many months before 
the vertical field computations were begun. Figure 8 shows the field 
of vertical motion computed for the 2,000 foot level (Figure 9 for the 
30,000.foot level) at 21Q0Z on 28 February 1954 or just two hours prior 
to the tin® of Figure 7» Comparing Figures 7 and 8 we find the greatest 
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convective activity (showers, ’thunderstorms and cumulonimbus) in the 
areas of upward motion. The similarity between the two patterns is 
quite striking. In a similar fashion comparing Figures 7 and 9 we 
find high cloud amounts greatest in the areas of upward motion at the 
30,000 foot level. Figure 10 shows the field of vertical motion at 
the 30,000 foot level at 0300Z on 1 March 1954 as computed by two 
independent methods, that of Dean and Ohmstead which is described in 
their report mentioned earlier and is based upon the calculation of 
horizontal divergence from streamline analyses and a method, suggested 
by Dr. C. E. Palmer, which I have U3ed and which will be described 
later. Inspection of Figure 10 will show that the vertical wind speeds 
obtained by the two methods are of the same order of magnitude and 
further, that the geometry of the two patterns is quite similar. From 
the studies which have been completed and were briefed in the preceding 
paragraphs, it appears to be possible to compute fields of vertical 
motion which are at least qualitatively representative of the actual 
atmospheric motions. 


Now I should like to discuss briefly the method of computation 
which we have used. There are in general two approaches to the compu¬ 
tation of vertical velocity; the thermodynamic and the kinematic. 
Thermodynamic methods are derived from the fact that changes in the 
vertical temperature and humidity structure are largely the result of 
vertical motions. Such changes as may be produced by non*-adiabatic 
processes such as radiative effects, the release of precipitation, etc., 
are either neglected or treated empirically. It is for this reason 
that the thermodynamic approach is not especially useful in tropical 
areas where non-adiabatic processes are maximized. Kinematic methods 
are derived from the continuity equation; they involve the computation 
of the divergence of the horizontal wind field at various levels and a 
vertical integration to determine the vertical velocities. The simplest 
kinematic methods involve the calculation of horizontal divergence with¬ 
in a polygon (usually a triangle) having observation stations at its 
vertices; the observed winds being averaged along each side of the 
polygon. This method produces good results in an area of dense obser¬ 
vations where the sides of the figure are not long but cannot be 
applied in the Marshall Islands area with any assurance of success 
because of the extremely sparse observing network. The method used by 
Dean and Ohmstead, as has been mentioned before, required a streamline 
(isogon-isotach) analysis at various, rather closely-spaced, levels. 
Divergence fields are computed from the streamline charts and the 


vertical velocities computed over a set of grid points from smoothed 


values of the divergence at each point. 
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developed by C. E. Palmer during GREENHOUSE. In Figure 11 is shown the 
derivation- of the computational function as I have used it. Using 
standard nomenclature the of continuity (l) is integrated in the verti¬ 
cal direction and by substitution the variable is changed to pressure 
resulting in equation (4). Defining the vector components u and v as 
indicated in equation (5) we may write the vertical velocity at the 
level of interest, w-j_ is equal to the horizontal divergence of the 
vector field 7 divided by the product of density times gravity. The 
computational, procedure used was developed from this relationship, 
equation (6). The two assumptions made in the derivation of equation 
(6), i.e., that the local rate of change of density and the vertical 
velocity at the surface are both negligible, have been shown to be 
essent ially valid over low oceanic stations in the tropics. 

To determine the u and v values for any particular case, the u and 
v components of the observed winds at each station available within the 
area of interest are plotted as a function of pressure as illustrated 
in Figure 12, In the figure, which is for Eniwetok at 21002 on 28 
February 1954, the heavier curves represent the u and v values. These 
curves are then graphically integrated with respect to pressure by the 
method of equal areas to determine the corresponding u and v values; a 
pressure interval of 5 centibars was used in the computation. The u 
and v values are then plotted as a function of pressure; the lighter 
curves in Figure 12 represent these graphs. Finally the value of u and 
v at any particular height can be read from the graphs by means of the 
pressure-height relationship. In Figure 12 the average values of the 
wind components u and v in each 5 centibar pressure interval are 
denoted as BJ and /yj respectively. In order to interpolate between 
the few stations available in the Marshall Islands area it has been 
found best to plot and analyze the U and v fields. This step is 
illustrated in Figure 13 showing the v field as solid isolines and the 
u field as broken isolines. An analysis must be made for each level at 
which it is desired to compute the field of vertical motion, but it 
should be emphasized that the u and v fields may be analyzed as scalar 
fields thus greatly simplifying the labor involved. From the charts 
of the u and v fields the field of Vertical motion at any level can be 
computed as indicated in Figure 14. A square grid was drawn on a trans¬ 
parent overlay in order to facilitate reading off the required u and v 
values; values were read off at each solid black dot shown in Figure 14» 
The average value of u aloqg the east side of square 3 is denoted as 

other averages are indicated by a similar nomenclature. Also shown 
in Figure 14 are the locations of the nine observing stations which were 
available in the Marshall Islands area for the cases investigated. 

From the preceding discussion, we conclude that if any appreciable 
percentage of the active material producing radioactive fall-out Is 
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carried above the 30,000 to 2^,000 loot level in the atomic cloud and 
if such material, or any large part of it, is carried on particles 
having a falling velocity less than about 5>000 feet per hour, then any 
quantitative calculation of tile fall-out pattern geometry and activity 
will certainly have to take into account the existing field of vertical 
motion within the area in question. Furthermore, if computations of th$ 
vertical field are required in connection with future test operations, 
it would appear to be feasible to compute them on a routine basis by 
means of the technique described in the preceding paragraphs. 
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CONCLUDING REMARKS 

LTCOL C* D. Bonnot, U5AF 
Joint Task Force SEVEN 


Figure 1 is a map of the Southwest Pacific area, centered on 
Eniwetok. This map presents the presently operational rawinsonde 
stations, as well as several new sites which we are proposing. You 
w ill note that Yap, Koror, Truk, Panape, and Majuro are to give us 
additional information; this was obtained through contract with the 
Weather Bureau ty the Task Force. There are several sites south of 
the equator which we are asking the British Commonwealth to establish. 
Also, there are two ocean station weather ships we are requesting the 
Coast Guard to operate for us on a continuing year-round basis. 

Lastly, you will note that we plan to establish temporary stations 
during test periods at Kapingamarangi, Kusaie, Bikati, and Taongj. 

This net will give us what we feel is an absolutely minimum observing 
net to support the Task Force. 

If such a weather observing net is operating, we feel that ade¬ 
quate information will be available to those required to predict the 
fall-out and radiation hazards at the Pacific Proving Ground, Work 
is being continued at the Oahu Research Center by Dr. Palmer toward 
improving our weather forecasting procedures for this area. 
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HtEDICTION OF DOSE-RATE AND DOSAGE CONTOURS AS 
FUNCTIONS OF YIELD AND METEOROLOGICAL CONDITIONS; 
EFFECTS OF SPECIFIC WING STRUCTURES USING THE ARMY 
SIGNAL CORPS METHOD 

Kenneth Barnett 
Office Chief Signal Officer 


Dr. Swingle has explained to you the Signal Corps procedure for 
making forecasts of the radiological fall-out patterns that result 
from l£ megaton explosions in the tropics# I would like to show yon 
three things. First, how this forecast procedure was calibrated by 
comparison with CASTLE Bravo observations. Second, the forecast made 
by using this procedure with different wind conditions, and third, some 
suggestions for making better fall-out forecasts in the future. In 
Figure 1, we have the trajectory of the 68 micron particles which fall 
out from the center of the cloud. These winds are the winds taken at 
shot time at CASTLE Bravo with some changes in the higher altitudes to 
account for the changes that occurred eight hours after the shot time. 
Since it required about one hour for the 68 microns particles to fall 
through £,000 feet, this trajectory is very similar to a hodograph. 

In Figure 2, we have the iso-c hr ones of the arrival time of the first 
particles end the ending times of the particles. Both are in hours, 
and using the predicted data IVos CASTLE Bravo. In the upper part 
we have two points where the fall-out time, the arrival time of the 
first particles, was ctoserved. We have one point at £.7 hours; it ia 
right on the £ hour predicted line. We have another point at 7*8 
hours that is almost exactly correct. The comparison of these few 
pieces of observed data fit quite well and give us confidence in the 
dimensions of our model. It is interesting to point out -that some¬ 
where between 30 and 35 hours the particles begin to drift backwards, 
as would be indicated by the X80 degree wind shift around £5,000 feet. 
Likewise in the lower curve, the ending time of particles, we also see 
how the ending time begins to change. In Figure 3, we have at the top 
the predicted dose rate, and this requires a calibration. As you will 
recall, we determine our dose rate by bringing out each of the disks 
or slices from our model according to the wind velocity and according 
to the rate of fall of the median particle size for that particular 
disk. We let these fall over a square, cross-sectional grid. We 
make a tally of the total number of cylinders that fall on each grid 
and then we draw lines for total nun±>ers of each of these tallies. 

In order that these total numbers for each square will be in roentgens 
per hour, we have compared the top figure of Figure 3 with the bottom 
figure, which is the observed data from AFSWP Report £07. You will 
notice that in the top figure, our predicted figure, as drawn here, 
our 100 line encloses an area a little larger than the area enclosed 
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by the 100 line in the bottom or observed figure. Likewise In our 
predicted pattern, the area enclosed by the l£00 r/hr line is slightly 
analle r than the observed area enclosed by the 1J>00 r/hr line. We 
used a calibration factor of 100 to obtain this tip figure and some¬ 
where, probably around the 600 roentgen per hour line, the areas are 
the same in the predicted model and in the observed areas. Looking at 
oer predicted dose rate, you will see that there is one problem rather 
obvious, that is the lumping of it. For example, out towards the end 
of the 100 line there's almost a perfect circle which indicates that 
OB& disc alone fell at that point. And right next to it, to the left, 
are the other circles. Later on 1*11 discuss some ways we have we 
hope will improve this. All of our calculations have been for the 15 
megaton model and we did not have the time to make models for one and 
megatons for this symposium. So ve used a simple scaling factor 
to give sane idea of what this procedure should give for a one megaton 
and 50 megaton yield. You'll see the two scales here in Figure 3 and 
in later Figures. In Figure U we have a comparison of our predicted 
total dose pattern with the idealized observed total dose pattern from 
CASTLE Bravo. The total area included by our predicted 100 total dose 
line and the idealized observed ICO total dose line is almost the 
ss®B» However, the area included by our IJjOO total dose line is a 
little smaller than the observed area. We think we have some pro¬ 
cedures that will improve this, that is, that will make our 1900 total 
dose line a little larger relative to our 100 predicted total dose 
line. Again here in Figure 1* we have these scaling factors to go from 
2_5 j megatons to either one megaton or $0 megatons. I would now like to 
go on to the examples that were furnished to us by AFSWP. You know 
their example 1 was the condition of an approximately 90 degree shear 
at a height of about It0,000 feet at Dodge City, Kansas. In the lower 
fi0,000 feet, the mean vector Is generally towards the south, from the 
ground seroj above that the vector goes out toward the east, and I 
think in our resulting patterns we should see some indications of these 
two vectors. In Figure 5, we have the isochrones for the arrival times 
and the ending times. The arrival time curves are quite symmetrical 
and go out steadily. Thi3 suggests that we might develop a short-cut 
method in which we extend one vector out along one azimuth and then 
space the isochrones at equal spaces there and draw them off, more or 
lass, as circles. I think we could use that to get an approximation 
teat would be very close to this and would save us a great deal of 
Mass. Looking at the predicted ending values we see that at one hour, 
there is a little irregularity which I think is the lower vector which 
points toward the south* In Figure 6, we have our predicted dose rates 
arid our predicted total dose* You will see that in both of those there 
is a little extension of the 15>0G area toward the south, which is a 
reflection of the fall-out in the lower 40,000 feet which is extended 
toward the south. And again we have our scaling factors here to give 
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the approximate dimensions for a one megaton and a 5?0 megaton yield. 
Going on to Case B, which is a gradual wind shear of approximately 90 
degrees (Figure 7), I think we should find a fairly symmetrical pattern 
because the mean wind vector from the surface up to various altitudes 
is almost constantly towards the east-southeast. In our Figure 3, we 
see that the arrival tints of the first particles, the isochrones for 
those and the ending time are quire symmetrical and extend out towards 
the east-southeast. Likewise, in Figure % the dose rate contours and 
total dose contours were quite symmetrical and point straight out 
towards the east-southeast. Case C should give vs very symmetrical 
results too, since the winds are from almost the same direction at all 
levels. And Figures 11 and 12 give the isochrones for ending and 
arrival times for dose rates and total dose. Case D, in my opinion, ie 
the most interesting one, and that is the condition of relatively low 
wind speeds at all levels. In fact, the effective wind from the base 
of the fireball on down to the ground is only about six miles per hour. 
Figure 14 shews the arrival time of the first particles and you notice 
that up to about 25? or 30 hours it is quite regular and then due to an 
almost 180 degree wind shift, the lines bunch up together and begin to 
go down towards the south. We likewise have a pattern for the siding 
time (Figure 15), which might be the super-posit ion of two symmetrical 
patterns, one on the other. Figure 16 shows by the dose rate contours 
that, at least the 1500 contour is almost a circle, with a long piece 
sticking up towards the northeast where the winds had begun to shift 
around and there was an unusually large amount of fall-out at that 
distance out in this particular case. Figure 17 shows us the total, 
dose and the contours here are nearly circles and up towards the north- 
cast those we have are very email, much ana Her than the dose rate 
values in that area. This is due to the deterioration factor because 
it is further out from the source. Figures 19 and 20 show the results 
of condition E, which is a case of relatively high wind speed above 
10,000 feet. (Figure 18) As would be expected, it is a long, thin 
pattern in the dose rate and total dose. Condition F is a case of a 
gradual shear of approximately ISO degrees (Figure 21), and again this 
might be broken down into two major vectorsj the lower half and the top 
half. Then we see that reflected in Figures 22, the idochrones of 
arrival time that are fairly symmetrical out to about 10 or 12 hours 
and after that the upper shear which goes down toward the south and 
makes itself apparent. In Figure 23 we have the dose rate and the 
total dose. Notice that on the right hand side there is quite a big 
lump out by itself In both cases. This is due partly to the wind 
shift at high levels going down towards the south. It is also the 
result, partly, of the lumpiness problem that we have in this model. 

In Figure 24 we have a summary of the dose rate and total dose areas 
for all of the conditions as well as for the observed CASTLE Bravo 
results and our predicted results. One interesting thing that should 
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be noticed is that for very high wind speeds, the area enclosed by the 
100 roentgen total dose line is much larger than in the low wind 
speeds* However, in very high wind speed cases the area enclosed by 
the higher roentgen values, such as 1500, are much smaller than those 
enclosed in the low wind speeds. I mentioned to you earlier that we 
had some suggestions for improving this forecast procedure and we 
have two. The first one is to get rid of the lumpiness by increasing 
the total number of discs in particle 3ize categories. To make any 
significant correction in this direction is going to require a great 
deal more work. For example, in working up the example of this type, 
if we use a 15 mile grid for tallying the dose rate, it requires 
three or four hours for an example. If we use a smaller grid, such as 
the 3^ mile grid which we employed for these examples, it requires 
from 15> to 20 man-hours. Now if we were to increase the disc and 
particle size categories, say ten times, it would increase our work 
ten times which is beginning to get unreasonable. The obvious answer 
to this is to apply it to machine methods which we hope to do* The 
second improvement we would like to make to our procedure applies to 
the case of the CASTLE Bravo whore we made predicted values using the 
composite winds in that test. Our results would have been better if 
the area enclosed by the 1^00 roentgen total dose lines had been better 
and larger relative to the area enclosed by the 100 roentgen total dose 
line,. I believe that this is partially the result of the assumption 
we made in developing our total dose formula. This assumption was that 
between the arrival time of the first particle and the ending time, at 
any one point, the accumulation of particles is constantly timed. Not?, 
this does not hold, and particularly near ground zero in the first few 
hours. I believe that our total dose values near ground zero would be 
significantly larger if we integrated instead of once, if we integrated 
several times, and particularly with each particle disc as it came 
down. Naturally this is going to be a great deal more work and again 
it implies the necessity for machine methods. One additional advant¬ 
age that might come from a procedure like that is that we can use a 
different deterioration rate with timej that is, in the first few hours 
we might use a or -1, and later on, after 8 to 10 hours, we might 

use values t“ ° or something of that nature. In ary case, we are 
attempting to do this and see if we can compute the CASTLE Bravo 
pattern by electronic calculating machines employing the improvements 
that I have just discussed. There are several obvious advantages to 
electronic calculationsj first it is much faster and this allows one 
to compute a great many cases of fall-out patterns whifch are necessary 
if we are to determine realistic probabilities. For example, suppose 
we wanted to know what the probability was of a lethal dose occurring 
50 miles west of a certain town. The only really reliable way to do 
that Is to compute daily patterns for several years and then add up 
these patterns and see what the probabilities are of having the lethal 




dose f>0 miles to the west. At the present we can make some approxi¬ 
mations to .this by using statistical wind data but this can lead to 
misunderstandings, if not used very carefully. Another advantage of 
computing these faster is that it admits the possibility of very de¬ 
tailed forecasts for tactical purposes; that is, ones that can be 
made in perhaps an hour or two. A last advantage to electronic 
computations i3 that it will make it practical for us to consider more 
factors in our procedure. For example, we should take more accurate 
account of the time space variations of wind during fall-out. We 
should take more accurate account of the stability of the atmosphere, 
dispersion, the terrain effects, and perhaps even precipitation. All 
of these things will greatly complicate the calculations and would 
only be practical with electronic computations. 

If we can look into the crystal ball for a few minutes here I would 
like to make some suggestions on something that might develop as a re¬ 
sult of the electronic calculations. Suppose we develop a series of 
wiring panels, each of which contains a different model; another series 
which contains a number of particle size distributions. Once we have 
those and have tested them to the point where we feel they are reliable 
then it would be possible to get a prediction of the fall-out pattern 
for my yield bomb for a wide variety of soil types and for any part of 
the world. For example, the yield of the explosion and the atmospheric 
stability would indicate the model size, so that the appropriate wiring 
panel for that model size could be chosen. Likewise, if we knew the 
soil type, w© could choose the particle size distribution for that soil 
type, we could choose the particle size distribution for that soil type 
and put that into the machine. Then our dose rate calibration factor, 
which you recall in our CASTLE Bravo prediction was 100, could be 
varied with the condition of the soil and the height of the explosion. 
That is, if the soil were frozen, there certainly would be much less 
soil sucked up in the explosion than if the soil were not frozen. 
Likewise, if the explosion were 500 or 1,000 feet above the surface, it 
probably would suck up much less dirt. The dose rate calibration 
facto*s take account of these things. DGE ARttZnrr 


One final point canes up and that is the question: Are we trying 
to do something that is more accurate than the data we are working 
with? I think perhaps we are. The data we have of the actual re¬ 
sulting fall-out patterns, some of the wind conditions, are not as 
accurate as we are trying to do. However, we can surely expect that 
this data will be more accurate in the future and ultimately we are 
going to need the most accurate prediction procedures possible. So I 
think we are justified in going ahead in trying to develop a fall-out 
prediction procedure as accurate as possible. 
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In closing I would like to say that I believe that we have a 
procedure that gives relatively good results and we have immediate 
prospects of some improvement. However, there are problems that 
need to be solved along with problems of this procedure. We need 
more information about soil types, particle size distribution, 
biological factors. We need more information about the winds to high 
levels. We need more accurate forecasts of winds. Last but not least, * 

we need more observed data of fall-out patterns to thoroughly test 
forecast procedures. 
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COUNTERMEASURES,* AN ANALYSIS OF NEW WEAPONS EFFECTS 
INFORMATION REQUIRED FOR FALL-OUT COUNTERMEASURES 

Paul C« Tompkins 

U. S. Naval Radiological Defense Laboratory 


A symposium on fall-out would not be complete if the need for 
information were not analyzed in the light of countermeasures needed 
to offset the effects. It is the purpose of this paper to consider 
that type of information which will lead to reliable countermeasures 
systems. 

The general conclusions derivable qualitatively from an analysis 
of existing information may be summarized as follows: 

Radioactively contaminated areas can be expected to be en¬ 
countered with sufficient frequency and at sufficient gamma intensity- 
levels to demand direct countermeasures. 

Existing facilities and equipment are not adequate to provide 
the necessary protection. 

Areas which can be affected are too large to permit use of 
evacuation or avoidance of a radiological area as the predominant 
control measure. 

Protection of both personnel and installations from the initial 
and subsequent effects supersedes recovery as a primary objective. 
Alternatives to protection and recovery from the effects of the event 
decrease as areas affected increase. 

The primary countena-sasur es system must be based on provision 
of sheltered living areas with shielding for protection against the 
gamma radiation, and recovery through pre-protection of an installation 
or its subsequent decontamination. 

The objective of countermeasures studies is to arrive at a deter¬ 
mination of the provisions that must be mads to permit people to live 
through a contaminating event, and still perforin essential functions. 
This implies that not only must one deal with the effects of radio¬ 
active contamination under emergency conditions, but one must also 
know how to live with radioactive contamination under long-term ex¬ 
posure conditions. 



Let us now turn to a consideration of the essential elements of 
countermeasures systems. These elements aret 

Reduction of gamma radiation intensities by shielding* 
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^ Reduction of garana radiation fields caused by loose radio- 
active weapon debris by removal of such debirs. This includes both 
post event reclamation methods and protective measures such as the 
"washdown" system employed by the Navy. 

Reduction of ingestion, inhalation, and beta radiation burn 
hazards by the use of proper doctrines, clothing, respirators and 
effective decontamination. 

Reduction of dosage by control of exposure time* This re¬ 
quires an intimate knowledge of time and motion involved on the part of 
personnel, and an equally intimate knowledge of the radiation intensity 
and radiological conditions under which these personnel must operate. 

Exclusion of contaminated material from critical areas such as 
those designed for occupany, food and water sources, etc* 

The first requirement of countermeasures is to provide enough pro¬ 
tection to permit survival during the acute phase immediately following 
a detonation. This must be done primarily through an adequate shelter 
program. Such shelters must provide a shielding factor of the order 
of 1000, if military forces are to retain the capability of operating 
despite the presence of fall-out. 

The second requirement is to establish the capability of recover¬ 
ing and occupying the affected regions as normally as possible* To do 
this it is necessary that "living areas" be designated. These areas 
must have maximum shielding available and provisions for excluding the 
entry of radioactively contaminated material. This requirement must be 
met if dosage control is to be possible. It is from these centers that 
personnel would then emerge into the surrounding contaminated area 
to discharge their assigned duties. The decontamination of "work 
areas" will be a requirement in many cases. The provision of clean 
food and water supplies will be the predominant practical problem en¬ 
countered in this phase. The acceptance of a single, sub-lethal 
acute dose as the criterion of acceptable dosage is not sound practice* 
Tt is required that a new set of medical standards be developed to 
cover continued exposures at various rates, comparable to the 30 
roentgens per two years which is the current standard for peacetime 
operations permitted by the Navy. Criteria for properly prorating 
the revised acceptance doses are needed. 

Unlike Industrial practices, fall-out influences the environment. 
Although this fact does not affect the currently practiced principles 
of radiological procedure, there are drastic practical influences on 
the conditions to which these principles are applied. Table I 
illustrates some of the more obvious differences brought about by tho 
radiological effect on environmental conditions. Notice that only 
items 9 and 10 stay the same under the two situations. 
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These differences are further illustrated by Figures 1 and 2. Note 
particularly the relationships between areas and contamination density. 
Figure 1 shows the lay-out of a building designed to conduct research 
work with multicurie quantities of radioactive materials. The highest 
activity area marked by the heavy stippling is the interior of the cells 
which are surrounded with 3-foot concrete walls for shielding, maze 
entrances to further reduce radiation levels when the cell doors are 
open, doors to confine activity to the interior, and air flow (shown by 
the arrows) for regions of low contamination to those of higher con¬ 
tamination. In Figure 2, the principles employed are identical, but 
the pressure is in the opposite direction. The situation in Figure 2 
is that which one encounters when one is required first to survive and 
subsequently live with the effects of a major fall-out event without 
evacuation. The same principles of ventilation control would lead to 
positive pressure ventilation in the shelter* 

This completes the summary of the qualitative picture. Let us now 
consider the more quantitative aspects which determine the character 
and ultimate cost of countermeasures systems. Quantitative knowledge 
is essential because any protective system can be overwhelmed by 
increasing either the gamma radiation levels or the contamination 
densities. Also, any system can safely handle any radiological con¬ 
dition lower than that for which it was designed. However, when 
applied unnecessarily, this usually involves very high costs and ex¬ 
cessive impediments to normal practice, since remote or semi-remote 
control operations are usually involved in coping with severe 
radiological conditions. 

Three types of protection are required: from whole-body gamma 
radiationj from skin burns arising predcminatly through direct con¬ 
tamination of the person; and internal deposition of radionuclides 
from all methods of entry into the body. 

It is recommended, with limitations to be noted presently, that the 
absolute maximum dosages considered acceptable for countermeasures 
systems design purposes be: 


For gamma exposure. 

Not over 5>00 r accumulated over 5 years. 

Not over 200 r in ary one year, 

Not over 20 r in any event under maximum protection. 
For beta exposure, 

Not over 5,000 rep in 5 years. 

Not over 2,000 rep in any one year. 

Not over 200 rep in any single exposure. 
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TABLE II ^ 

Reference Contours Related to Delivery of 1000 Rep 
or 10,000 Rep in One Week 


Time After Burst 
of Start 

of Initial Exposure 

Reference Contour ^J 
Leading to 1000 rep in 

1 week (r/hr at 1 hr) 

Reference Contour 
Leading to 10,000 rep 
in 1 week (r/hr at 1 hr) 

1 hr 

7.8 

78 

1 wk 

107 

1,070 

1 mo 

U55 

U,55o 

1 yr 

7,5oo 

75 ,ooo 


(a) Calculations by Broido, A,, and Teresi, J, 

(b) The calculations assume the applicability of the t” 1 *^ decay lav, 
and uniform contamination over the whole body* Whole-body beta 
exposure UV Q assumed to be about £000 rep. The 10,003 rep 
column applies to small particles only. 


DOE ARCHIVES 




For internal deposition, 

89-90 

The accumulation of 5 Me Sr 7 as the body burden should be 
considered the absolute upper limit* 

It must be stressed that these criteria have a limited biological 
safety factor and those personnel getting the maximum of any one of 
these might be expected to show a high incidence of delayed effects. 
Further, a quantitative basis for reduction in allowable whole-body 
exposure in the presence of radionuclides deposited in the bone is "not 
available. 

The problem of controlling gamma exposure has been dealt with many 
times* Suffice it to say that exposure time control is the essential 
ingredient and that design appears to call for a shielding factor of 
the order of 1000 or more. 

The probable scope of the beta burn problem is shown in Table II. 
Here one sees an estimate of the reference contour lines expressed in 
r/hr at 1 hr which would be expected to give 1000 rep and 10,000 rep 
in one week. The sele.ction of these two levels was based on the 
assumptions that: 

The burns would come largely from weapon debris in direct- 
contact with the skin. 

No iimediate burn would be developed with less than 1000 rep. 

A local burn would be certain above 10,000 rep from a snail 
particle (point source). The current effects table (AFSWP) states 
that incapacitation will, occur in 3 to 6 hours from 3000 to 5000 rep 
over a relatively small area (a few square centimeters) of the body. 

Health habits will insure that no contamination stays in one 
particular locality on the boc^r for more than 1 week. 

No subsequent contaminating material will expose exactly the 
same area of the skin by depositing in exactly the same place. 


Note that the reference value shifts upward as the initial time of 
exposure increases from T . 

° DOE ARCHIVES 

The internal deposition problem can be assessed from the Sunshine 
project which reached the conclusion that the Sr' u tolerance level in 
soil would be reached at about the 120 r/hr at 1 hour reference line. 
During the first 18 months, the activity of Ba li: - La 111 and Sr ^ is 
about 10 times this amount or more. Allowing for build-up from de¬ 
tonations occurring over a period of time and the possible addition 
of inhalation, as a source of internal deposition, reference contour 
of 1 r/hr at 1 hour can be taken as defining the point at which one 
must start becoming concerned. 
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The general character of the applicable countermeasures systems e 

as related to the reference contour lines from a single detonation are 
summarized In Table HI, 



Biological: 


Maximum "acceptable" dosages prorated over ? years. Particular 
attention should be paid to late biological effects occurring between 
*>-10 years. 

Medical Effects Table for skin effects as functions of time 
and area. Particular emphasis should be paid to burns from small "hot" 
particles. 

Instruments: 


Contamination meters for. clothing, water, air, food, etc. 
must be able to cancel out surrounding radiation field between — 0.1 

mr/hr to 1 r/hr. 

Dosimetry should distinguish between deep dose and surface 

dose. 


Calibrations with a point source appear to be inadequate. 


Calibrations should reflect the distributed geometry which is actually 


encountered. 


Characteristics of radionuclides: 


Beta energy and mount as function of time up to 2 years. 
Gamma energy and amount as function of time up to 2 years. 
Fission Yields. , < 

Solability in the range of 10”" 5 to 10 • 

Exchange properties with soil and other surfaces, particular¬ 
ly Sr and Pu. 


Uptake of radicelemerrts into growing foods. 


d-oe archives 


Radiation Fields characteristics: 


Radiation fields from contamination distribution other than 
the infinite flat plane. Knowledge of important variations in either 
the spectrum or air ionization are needed. 

A "simulated" or "synthetic" gamma spectrum as a function of 
time is needed for both dosimetry and shielding development work. 

Typical radiation fields from basic geometric models should 
be developed schematically so that actual situations can be analysed 
by a summation of parts. 
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Shielding: 

Basic types of construction geometry should be schematized. 

Build-up factors should be measured for energies of about 
O.b Mev, 0.7 Mev, and 1.5 to 2,0 Mev. 

Build-up factors should be based on data for slab, edge, and 
corner orientations. 

Fall-out: 

Mathematical model should cover region to the 1 r/hr at 1 hr. 
reference contour line. 

Meteorological conditions leading to deposition of "hot 1 * 
spots in lower activity areas should be identified possibly through 
refinement of the model. Probable number and size but not location 
should be statistically predictable. 

Chemical composition of the weapon debris is extremely im¬ 
portant in relation to contamination density and solubility. 

Countermeasures; 

Reduction in the radive field without shielding by removal 
of the weapon debris should aim at 90 to 99 per cent reduction. 

Water treatment should aim at a reduction of~lC)5 in total 
activity and~lCP for Sr. 

Contaminability of clothing and skin from the environment 
must be quantitated at least roughly. 


£>ge archives 




TABLE III 

Comparison of Countermeasures Requirements as Related to Initial Contamination Level 
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TABLE IIT (continued) 

Comparison of Countermeasures Requirements as Related to Initial Contamination Level 
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PREDICTION OF DOSE RATE AND DOSAGE CONTOURS AS FUNCTIONS 
OF HELD AND METEOROLOGICAL CONDITIONSj COMPARISON OF 

PROBLEM SOLUTIONS 

GDR Roger ¥ Paine, Jr., USN 
Headquarters, Armed Forces Special Weapons Project 


Early in the planning phase of the fall-out symposium, it 
was decided that a ready means of comparing the various methods 
of generating fall-out contours would be a desirable addition 
to the symposium summary report, and to the symposium itself If 
such a means of comparison could be made ready in time. In the 
last analysis, the contours generated by a model for a specific 
weapon yield and for a specific we at ter situation probably form 
the best means of comparison with other models. Accordingly, 
speakers who had consented to present their methods at the sym¬ 
posium under Topic F were sent six actual weather soundings 
taken at various localities in the United States in recent years 
and asked, time permitting, to work out a series of contour 
solutions for two or more of these weather situations. In order 


to make & distinct break with the CASTLE BRAVO event, against 
which most of the methods had already been calibrated, yields 
of 1 MT and 50 HT were chosen. The participants were asked to 
generate by thsir methods for surface bursts of each of these 
yields, contours of 100, 500, end 1500 r/hr referred back to 


H-l hour, and contours of 100, 500, and 1500 
la ted up to U8 hours after burst time. 


r total dose accumu« 
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This "homework", as the participants chose to call it, rep¬ 
resents a considerable amount of hard work by the agencies involved, 
some of whose methods are not easily adaptable to generate the 
particular contours asked for. It is therefore very gratifying 
that couplets problem solutions by each method were turned in for 
one of the weather situations, and at least partial solutions by 
each method for a so cord. Particular thanks are due the Arssy 
Signal Corps, the Air Weather Service, and Technical Operations 
Incorporated for completing work cm all six weather situations. 


Owing to the great volume of material involved, only prob- 
Isa solutions for the two situations for which almost complete 
returns are available will be compared in detail in this report. 
Solutions for these situations are sufficient to point out the 
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more important conclusions which can be drawn from the conparisons, 
which are; 

For situations which do not involve abrupt wind shears, the 
contours generated by the several methois, although differing in de¬ 
tail, do not differ greatly in the direction of fall-out nor in the 
order of magnitude of the area covered. 

Because of the basically different assumptions as to the loca¬ 
tion of radioactive materials in the stabilized cloud, the presence of 
an abrupt wind shear between 15,000 and 50,000 feet is apt to cause 
considerable differences in the predicted direction of fall-out for 
certain contours,although predictions of the magnitudes of the areas 
covered do not differ greatly. 

None of the methods predict very high intensity dose or dose 
rate contours, even for the hypothetical 50 MT yield. A maximum dose 
rate of the order of 10,000 r/hr or less at H-t-1 hour appears to be 
indicated, and this is likely to occur only in a very small area at or 
near surface zero. 

In comparing the problem solutions, it is well to keep in 
mind the basic assumptions and limitations of the methods being com¬ 
pared. Some of these account directly for the differences observed. 
Examples are: 

AFSWP Method 

This method is based on ground surveys at JANGLE-S and CASTLE 
BRAV0, and consists of "idealized 1 ’ contours, which follow a single 
scaling wind direction. Abrupt wind shears and unusual weather con¬ 
ditions are not easily handled. The method is suitable for planning 
purposes only; not for post-shot analysis. An order-of-magnitude 
contour can be drawn for any weapon yield between 0.1 KT and 100 MT 
in two or three minutes by untrained personnel. 

AKDC Method 

The model assumes that 90 % of the bomb debris activity is in the 
stem of the stabilized cloud, and 10^ is in the mushroom. Mean 
effective particle sizes are assumed for the cloud and parts of the 
stem, and Stokes Law fall rates for spherical particles are used. Wind 
and weather conditions are allowed for. The method is calibrated to 
CASTLE BRAVO, but is adaptable over a wide range of yields. A prob¬ 
lem solution requires several hours by trained personnel. 










Amy Signal Corps Method 


The model divides a hypothetical stabilized bomb cloud consist ing 
of super-posed cylinders into disc or cylindrical wafers or compart¬ 
ments, each associated with a particular particle size category and fall 
rate. Each disc is then brought to the ground according to the winds 
acting on it, and ground values are then summed over all the discs to 
obtain contours. A different model must be generated for each weapon 
yield, and for different localities (tropical, polar, etc). Wind and 
weather conditions are allowed for. A particle size distribution 
hypothesized by the Rand Corporation is used. Aerodynamic particle 
fall rates are used. The method is calibrated to the AFSWP-507 report¬ 
ed survey of CASTIE BRAVO. Because of the finite number of discs in 
the model, contours are often lumpy and may involve artifact "hot spots". 
A problem solution requires several hours by trained personnel. (The 
dose contours presented are for 18 hours after arrival time, rather than 
for H+18 hours). 

Air Weather Service Method 

The method is essentially a radex plot resulting in rough sectors 
within which fall-out may be expected, together with estimated times 
of arrival along vectors from ground zero for particular altitudes from 
which it is assumed particles start their fall. The method does not 
generate contours, and thus is not directly comparable to methods that 
do* The plot can be drawn in about 15 minutes by trained personnel* 

Los Alamos Method 

The method was devised for use during CASTLE in forecasting the 
gross results to be expected from surface bursts of about 10 MT fired 
in the northern Marshalls. The plots submitted are modified only by 
direct yield dependence. The effect of yield on initial cloud geometry 
is not included, nor is the effect of the different tropopause height 
of the middle latitudes from that in the Marshalls. The method gives • 
very little detail anywhere and none at all near the origin. It is not 
intended for use in detailed post-shot analyses, but rather for radsafe 
planning. A problem solution probably requires about an hour for an 
individual familiar with the method. 

Navy Method DOE ARCHIVES 

This method, which is based largely upon a radex wind plot, is a 
system of weighting incremental square areas according to the degree 
of fall-out expected, relative one to another. It results in wind- 
sensitive contours which have only relative values* Actual values may 
be fitted to the contours as a result of one or more post-shot survey 
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readings taken in the contaminated area. For purposes of this com¬ 
parison, the contours submitted have been roughly normalized in area 
to the corresponding contours of another method or methods (e.g., LASL) 
in order to show agreement in shape. A problem solution requires some¬ 
thing more than an hour for personnel familiar with the method. 

NRDL Method 

This is a wind and weather dependent system which assumes the bulk 
of the radioactive material originates in the lower portion of the mush¬ 
room, and utilizes an aerodynamic particle fall rate which varies con¬ 
siderably with the altitude of the particle. A problem solution 
requires several hours by personnel familiar with the method. 

Rand Corporation Method 

The method gives wind and weather-sensitive contours, based on an 
assumed particle size distribution and the hypothesis that 90$ of the 
fission product radioactivity falls out from the mushroom cloud and 10$ 
from the stem. An aerodynamic rate of fall is used which is somewhat 
different from that used by NHDL, but which also varies markedly with 
altitude. A problem solution requires a great amount of hand calcu¬ 
lation. The method has also been programmed for machine solution, and 
the machine method was used for the solutions presented here. 

Technical Operations, Incorporated, Method 

The method utilizes an inverted cone cloud model, the NRDL assump¬ 
tion of particle size distribution,_and fall rates which increase with 
incr ease in .altitude. 

A problem solution requires some w^TTra'cS'e’’-'^ 
Hrcah ajTTfour^'utilizing^trained personnel. 

The problem situations presented for comparative solutions were 
as follows; 
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CONDITION "A" 


Wintertime situation of an abrupt, approximately 90° shear at a height of 
approximately 1 * 0,000 feet. 


Dodge City, Kansas 

37° 1*6'N, 99° 58' W. 



1500 GMT 

Elevation: 

28 Dec 1953 

2625 feet 

Height 

ft, msl 

Wind Biree- Wind Speed 
tion Degrees Knots 

Pressure 

Millibars 

Temperature 
Degrees C 

Relative 

Surface 

260 

7 

925 

- 7 »li 

8 U 

5,000 

21*7 

12 

81*2 

- 3.0 

36 

10,000 

273 

19 

695 

- 8.7 

30 

15,000 

307 

13 

570 

-13.2 

MB 

20,000 

008 

16 

1*63 

• 2 lj .0 


25,000 

0U5 

1*1 

373 

-314.2 


30,000 

036 

52 

300 

-1*5.5 


35,000 

357 

29 

238 

-53.5 


1 ( 0,000 

21(3 

1*7 

187 

-55.5 


h5,ooo 

255 

1*3 

11(8 

-56.0 


50,000 

255 

55 

115 

-60.2 


55,000 

260 

1*7 

90 

— 60.6 


60,000 

272 

51* 

71 

— 6 l *0 


65,000 

269 

ho 

55 

-58.8 


70,000 

265 

36 

1* 3 

-59.3 


75,000 

285 

38 

3U 

-56*5 


80,000 

285 

U5 

27 

-56.0 
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CONDITION "B" 


Gradual Shear of Approximately 90° 


Washington, D. C, (5ilver Hill) 

38° 50« N, ? 6 ° 57* W. 

0300 GMT 28 Sept 1952 

Elevation: 289 feet 

Height 

Wind Direc¬ 

Wind Speed 

Pressure 

Temperature Relative 

Ft, msl 

tion Degrees 

Knots 

Millibars 

Degrees C Humidity 

Surface 

Calm 

0 

1010 

9.7 85 

5,000 

356 

11 

853 

11.8 16 

10,000 

009 

20 

708 

5.8 39. 

15,000 

325 

li» 

688 

- 0.2 MB 

20,000 

282 

19 

U85 

-10.9 

25,000 

263 

3U 

391; 

- 21 . 1 ; 

30,000 

263 

U7 

320 

-31;.0 

35,000 

273 

37 

255 

-Hi.o 

lj 0,000 

308 

27 

203 

- 51.5 

1j5,ooo 

292 

23 

160 

-60.2 

50,000 

290 

26 

123 

-66.5 

55,000 

290 

20 

97 

-65*1; 

60,000 

268 

11 

76 

-6U.0 

65,000 

276 

21 

58 

-61.5 

70,000 

293 

7 

he 

-59.0 

75,000 

293 

8 

36 

-55.5 

80,000 

285 

10 

26 . 

-51;.3 

65,000 

250 

11 

23 

- 52.0 

The 

salient features to note about Condition "A 11 are that the 

winds are 

generally from the west, 

except for a band between 15,000 


and 35,000 feet, wherein the wind3 are from the north and the north¬ 
east, The mean wind is about I;0 knots, _ _ _ 

doe archives 

For Condition "B", the winds are again generally from the west* > 

except for the region between the surface and about 15,000 feet, within 
which the winds are northerly. The mean wind is about 20 knots# 

The problem solutions are compared by superposition on the charts 
which follow, using coded overlays as indicated on each chart. Miss- -■*> 

ing solutions either had not been received in Headquarters, AFSWP, at 
time of press make-xip (1 February 1955), or the material received was 
not readily adaptable for comparison. Tables giving numerical compari¬ 
sons of contour areas and downwind and crosswind extents follow the 
charts. 
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Comparison of Certain Contour Parameters 
CONDITION "A" 

1 KT Surface Burst 
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* Indicates contours not received for this condition at time of press make-up. 

•sBi- Indicates method does not predict a significant contour for the conditions specified. 

( ) Indicates smallest value predicted. 

_ _ Underlining indicates largest value predicted. (Navy and Tech Ops values not Included 
in this comparison, since the contours had to be normalized from data furnished.) 
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CLOSING REMARKS 


v 

Dr. Herbert ScoviUe, Jr, 

/- Armed Forces Special Weapons Project 
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Before opening the general discussion of the conference, I should 
like to summarize briefly my cwn impressions of what has been accomp¬ 
lished, First, we should remember that the prime objectives of this 
symposium were to make all the data on fallout available to everyone 
who has been making studies in the field. The second objective was 
to put forth the various theories and fallout models so that their 
assumptions and methods could be compared. Obviously, it is not pos¬ 
sible to reach any decision as to a best method at a meeting of this 
sort. Actually, there is probably no such thing as a best method 
since different groups have different requirements for information. 

I believe this symposium has rather graphically shown, in case 
anyone had previously doubted it, that the radiological fallout in¬ 
volves very complex phenomena. The initial stages in the process- 
depend upon reactions occurring in the fireball and the properties of 
the atomic cloui, Later stages include atmospheric transport problems, 
particularly those related to the settling velocity of particulate 
matter. Finally, one has the radiological phenomena which determine 
the radiation fields in actual locations and which are effected by the 
radiochemical properties of the bomb debris. The entire fallout pro¬ 
cess is moreover strongly dependent on such external variables as the 
nature of the soil at the point of detonation and the meteorological 
conditions. It is, therefore, no wonder that all the data are not 
self-consistent and that a variety of models can be developed to ex¬ 
plain the test observations. 

With respect to the data, I believe that it is important to 
remember that to a large extent it derives from one or two actual test 
shots. For example, the data from the BRAVO shot has been used for 
calibration of all theories. Other shots have produced bits and 
pieces of data which are very useful in checking parts of the theories 
but with the exception of the JANGLE surface shot, all others have 
certain basic differences from the situation of greatest operational 
interest; i.e«, surface land detonations of high yield weapons. Fur¬ 
thermore, as Lt Colonel Lulejian has pointed out, the weather conditions 
which actually existed at the ERAVO shot were not good for differenti¬ 
ating between soma of the proposed models. In particular it was not 
satisfactory for determining what fractions of tie active material are 
in the cloud and stem* I think it is also important to remember that 
obtaining data on fallout at the Pacific Proving Grounds is an extreme¬ 
ly difficult process, A great deal of credit is due to those who have 
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struggled to get this information, where most of the fallout was occur¬ 
ring over the open ocean. In this respect, there is a new bright ray 
of light for the future as a result of the very fine work which was 
done on the spur of the moment, using the ocean itself as a collecting 
medium. It is indeed a rare event when a project produces important 
results within a week of its conception, and all who participated deserve 
a great deal of credit for the success of the oceanographic sampling 
project. In tire future it would seem that oceanographic sampling com¬ 
bined with aerial surveys would provide a very practical solution to 
obtaining good fallout data at the Pacific Proving Grounds. 

Despite all the difficulties in obtaining good data and the com¬ 
plexity of the phenomena, I think it would be a mistake to get the 
impression that we are completely ignorant of the magnitude of the fall¬ 
out problem in operational situations. Actually, I believe that the 
general scale of the pnenorasna, as derived from the BRAVO fallout data, 
can be believed with a great deal of confidence. The available data 
from both the JANGLE surface shot which took place over a silicate sand 
and the ocean data from the lagoon barge shots all support the thesis 
that an appreciable fraction of the radioactive material produced in 
the explosion will fall out in the neighborhood of the detonation. We 
must conclude that a surface burst of a yield coup arable to that of 
BRAVO will produce large areas, approximately 5,000 square miles, where 
the fallout will be lethal to exposed personnel. In certain specific 
situations this area might be only half as great and in others the area 
might be even greater, but I believe it can be guaranteed that this area 
will not be limited to a few hundred square miles. Likewise, I believe $ 

it is important to note that none of the speakers indicated that the 
contamination level outside the area of blast and fire damage would be 
likely to exceed ten or twenty tines this lethal level. This is extreme¬ 
ly important from a defensive point of view since it means that adequate 
radiation shielding ban be obtained anywhere without construction of 
massive, expensive shelters. 


Coining to the models, one has seen during the past two days a 
large number with widely varying degrees of complexity. Some will solve 
problems in a few minutes; others will require hours or fancy machines 
to produce the specific fallout patterns. Each of these models have 
their specific operational usefulness. While the complicated medals 
can not be used directly to give operational information immediately 
following an attack, they are particularly useful in arriving at a 
detailed understanding of the phenomena. Furthermore, if a wide variety 
of typical meteorological situations are studied by means of these 
machine calculations it may be possible to get a rough picture which 
would be satisfactory for operational purposes by merely comparing the 
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actual wind conditions with those which had been previously used in a 
machine solution. This approach would be similar to that described by 
Dr« Felt, who indicated that at the Nevada Test Site it is possible to 
predict with considerable confidence, based on previous experience, the 
fallout levels which will occur. 

In sum!nary, I believe that we can conclude that: (a) the general 
scale of the fallout problem under various detonation conditions can be 
reasonably well defined, (b) a number of models are now available which 
will provide reasonable estimates of the specific fallout patterns which 
are likely to occur under different meteorological conditions, and fin¬ 
ally, (c) there is a real need for further experimental data to check 
the various theories. In this connection, the oceanographic techniques 
offer considerable promise in the Pacific, but truly satisfactory re¬ 
sults win be obtained whan it is possible to detonate large yield 
weapons on the surface of a large land mass. 
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GENERAL DISCUSSION 


Question: Dr, Kellogg, Rand Corp. 

Is there any initial dosage rate at -which it would be 
advisable to evacuate a contaminated area? It would seem 
that there must be a limit somewhere. 

Answer: Dr, P, C. Tompkins, USNEDL 

No, In ny opinion, it is completely feasible to reduce the 
radiological effects at any location down to a safe value so 
that the limiting factor will be blast or fire. I am equally 
sure that one can design shelters to provide blast protection 
quite close to the limits of the crater, but here the problem 
may become economic. 

Question: (Unknown) 

Does this hold for civilian populations too? 

Answer: Dr, P. C. Tompkins, U5NRDL 

In principle yesj however, practically there may be situations 
where evacuation is the best. 


Question: Dr, Larger, FCDA 


When you said you rule out evacuation, does this mean you 
rule out permanent evacuation also? 


Answer: Dr, P, C. Tompkins, U5KRDL 


DOE 


In general, yes, I do not believe that we are going to run 
off and leave 5,000, 10,000, or 15,000 square miles deserted. 
If a single vreapon were going to be dropped, then evacuation 
might be more practical but in a situation where many weapons 
may be released then this measure does not appear too useful. 


Comment: LTCOL Lulejian, ARDC 


I should like to show this slide of 111 bombs exploded over 
this country, 106 on cities of 100,000 or more. This points 
out that in such a situation there are very few locations to 
which it is safe to evacuate. Furthermore, it points out 
that when you have large numbers of bombs, meteorology becomes 
less important since the over-lap tend3 to blanket such large 
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Comment; Dr. Knapp, OEG - 

This last chart certainly makes evacuation not look very 
good, up in the north, anyway. However, I would like to 
point out that people living within 5 or 6 miles of the 
center of a city might be killed fcy blast. The only 
solution for these people if thqp have an hour or two 
warning time might be to go upwind and take shelter. 

Comment; Dr. Felt, LASL 

In the event that we find ourselves forced to live in a 
contaminated area following an attack, I might point out 
that about twenty-five of us lived almost two months in 
such an area and after the first week people adjusted to it 
quite readily. After a couple of mistakes, it never occur¬ 
red to us to go outside without taking necessary precautions. 

Comment; LCDR Miller, BuAer 

If one is planning to measure fall-out by surveying wide 
areas with aircraft, I believe that the accuracy of the 
method should be checked more thoroughly, experience 
at CASTLE indicated considerable variation in the correction 
factors to reduce aerial readings to surface values. 


Comment; 


Dr. Werner, USNRDL 
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I believe that the aerial techniques for obtaining the 
intensity over the water are reasonably well established. 
However, care should be exercised to avoid the impression 
that these can give you the total amount of material in 
the water without additional information on the distri¬ 
bution of activity with depth. Perhaps the apparent cor¬ 
relation which was observed at CASTLE resulted from a con¬ 
stant distribution with depth which can not be relied on 
in all cases. 


Comment; Mr. LeVine, NYOO, AEG 


We too had difficulty initially in correlating aerial data 
with ground surveys but by the last shot the techniques had 
been sufficiently well developed so that we obtained a con¬ 
stant correction factor in almost every case. It should be 
noted that the aerial surveys picked up activity which went 
off to the northeast and which was not detected by other 
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means. This activity disappeared on the second day indicat¬ 
ing that particles had settled below the surface. There is, 
of course, considerable question as to how much activity had 
'disappeared before the first survey was made. This observa¬ 
tion may give some clue as to the presence of activity in the 
stem since the fall-out pattern from the stem was to the 
northeast in contrast to the pattern from the cloud which 
was to the west. 

Comment: Dr. Scoville, AFSWP 

One of the reasons for the discrepancies in aerial survey 
measurements results from the non-uniformity of the source. 
When a large area is uniformly contaminated, the altitude 
correction factors should be constant. However, for the 
Ihiwetok situation, where the contamination is heavy on a 
small land mass and low over the water, it is quite natural 
to expect big variations from one location to the other. I 
do believe, however, that it will always be necessary to 
calibrate the aerial surveys over the ocean with oceano¬ 
graphic measurements in order to obtain depth distribution. 

Comment: LTCOL Lulejian, ARDC 


Question: 


We did lots of aerial measurements at Nevada and these show 
the importance of the geometric factors. If the contamina¬ 
tion is not uniform or if there is shielding, then the 
altitude relationships will vary. 


MAJ Hord, AFSWP 


DOE ARCHIVES 


What is there about the thermocline which traps the fission 
products? 


Answer: 


Mr. Isaacs, Scripps 


In the ocean there is a layer, about 100 meters thick, in 
which there is considerable mixing; but there is very little 
mixing between this layer with the water below it. There¬ 
fore, diffusion of fission products through this boundary, 
i.e., the thermocline, is very slow. 


Question: Dr. Scoville, AFSWP 

Both the shots at which oceanographic measurements were made 
were lagoon shots where the fall-out was primarily associated 
with very small particles. In a surface land shot where the 
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particles are ICO micron or larger, would they be expected 
to have an appreciable settling velocity through the 
thermocline? 

Answer: Dr. Folsom, Scripps 

We do not know. In order to check this at an actual shot it 
will be necessary to collect water samples below the thermo¬ 
cline and, where possible, even on the bottom. 

Comment: Dr. Kellogg, Rand Corp. 

With respect to fall-out model construction, it is interest¬ 
ing to note that a great many agencies have started by using 
* a simple model, finding that it does not satisfy the experi¬ 

mental results and finally going to more and more complex 
models, including eventually machine calculations. It should 
be pointed out that the machine calculations, once they are 
developed, will again reverse the trend and decrease the time 
required to solve a specific fall-out problem. 

Comment: Dr, Graves, IASL 

I think what Dr. Kellogg said is true. It is possible that 
complicated calculations may now be made by machine. On the 
other hand, we have seen examples of various computations 
here and they do not disagree enough to be able to evaluate 
them with the available data. I winder whether we ars yet 
at the stage where these more complicated calculations are 
justified until we have more input data. 


Question: Mr. Isaacs, Scripps 


DOE ARCHIVES 


In this connection, one matter worries me. This is the type 
of material over which the bomb is exploded. We have talked 
about particles produced from soil but soil does not enter 
into it if one has a surface burst in New York City. In such 
a case, what would be the particle size? I also question how 
realistic a very careful model would be. 


Answer: Mr, Tompkins, CRL 


As far as cities are concerned the Chemical Corps is looking 
into this now. I think this problem can be approached by 
trying to develop a better understanding of the mechanism 
of the formation of the original particles in the cooling 





Comment: 


Comment: 


Question: 


Answer: 


fireball. If we understand this mechanism, we should be 
"better able to predict what will happen for different types 
of surfaces. 

Dr. Felt, LASL 

I think that Dr. Scovilla mistakenly gave the impression that 
it is not worthwhile for people to consider data resulting 
from other types of shots, particularly those at Nevada. 

These produced lots of information, some good, some poor. 

The models should be able to say something about the Nevada 
shots if they are good models, and these shots provide a lot 
more data for checking the various systems. 

Dr, Scoville, AFSWP 

I agree. I only meant to imply that the only test which 
approaches the situation which one is most interested in 
operationally was the Bravo shot, but certainly other shots 
can be useful in providing checks on various steps in the 
different fall-out models. 


Dr. Tompkins, USNRDL 


I wonder if Mr. Young of NOL would be willing to discuss some 
of his work on mass subsidence which might explain some of 
the other fall-out observations. 


Mr. Young, NOL 


DOE ARCHIVES 


We have observed mass subsidence from HE tests and underwater 
shots. A few pictures we have seen at CASTLE indicated that 
this may have occurred during the first five minutes but un¬ 
fortunately there was insufficient light to get good pictures. 
This effect might have contributed to the early spreading out 
of the active material along the surface. 
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